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Abstract 



Monopoles are solutions of an SU{2) gauge theory in satisfying a lower bound for en- 
ergy and certain asymptotic conditions, which translate as topological properties encoded 
in their charge. Using methods from integrable systems, monopoles can be described in 
algebraic-geometric terms via their spectral curve, i.e. an algebraic curve, given as a poly- 
nomial P in two complex variables, satisfying certain constraints. In this thesis we focus on 
the Ercolani-Sinha formulation, where the coefficients of P have to satisfy the Ercolani-Sinha 
constraints, given as relations amongst periods. 

In this thesis a particular class of such monopoles is studied, namely charge 3 monopoles 
with a symmetry by C3, the cyclic group of order 3. This class of cyclic 3-monopoles is de- 
scribed by the genus 4 spectral curve X, subject to the Ercolani-Sinha constraints: the aim of 
the present work is to establish the existence of such monopoles, which translates into solv- 
ing the Ercolani-Sinha constraints for X. 

Exploiting the symmetry of the system, we manage to recast the problem entirely in terms 
of a genus 2 hyperelliptic curve X, the (unbranched) quotient of X by C3 . A crucial step to 
this aim involves finding a basis for ifi (X , Z) , with particular symmetry properties according 
to a theorem of Fay. This gives a simple form for the period matrix of X ; moreover, results by 
Fay and Accola are used to reduce the Ercolani-Sinha constraints to hyperelliptic ones on X. 
We solve these constraints on X numerically, by iteration using the tetrahedral monopole so- 
lution as starting point in the moduli space. We use the Arithmetic- Geometric Mean method 
to find the periods on X: this method is well understood for a genus 2 curve with real branch- 
points; in this work we propose an extension to the situation where the branchpoints appear 
in complex conjugate pairs, which is the case for X. 

We are hence able to establish the existence of a curve of solutions corresponding to cyclic 
3-monopoles. 
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Introduction 



The concept of a magnetic monopole as a pointlike magnetic charge was introduced by 
Dirac almost 80 years ago in [DirSl], in the context of electromagnetism. It was a ground- 
breaking idea, as for the first time it put electricity and magnetism truly on the same footing 
in Maxwell's equations, and more importantly provided an argument for the quantisation of 
electric charge. 

With the study of nonabelian gauge theories, where the U{1) gauge group of Maxwell's theory 
is enlarged to a nonabelian one, it was soon realised that objects like the Dirac monopole are 
not just peculiar of electromagnetism. In the case of SU (2), 't Hooft [tH74] and Polyakov 
[Pol74] numerically found that the Yang-Mills equations admit soliton solutions, namely 
smooth fields configurations localised in a region of space, which at large distances be- 
have like Dirac monopoles. Shortly afterwards, in [PS75] Prasad and Sommerfield exhib- 
ited analytic solutions in a simplified case, with vanishing Higgs potential; Bogomolny stud- 
ied this case further [Bog76], introducing the Bogomolny equations (1). Hence this class of 
monopoles is referred to as BPS monopoles: these are the monopoles we are concerned with 
in the present thesis. 

Interestingly enough, monopoles have never been found experimentally. Nevertheless, 
the physical interest in the consequences of the existence of such objects, and, more rele- 
vant to us, the beauty and richness of their mathematical description make monopoles a 
rich and dynamic area of research still nowadays. 

Consider an SU{2) gauge theory on the Euclidean M?, with connection A, and Higgs field 
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(j). Monopoles are then solutions of the Bogomolny equations 

F = (1) 

(with certain boundary conditions); here F is the field strength and X> the covariant deriva- 
tive associated to a gauge field A, while ★ is the Hodge star. 

The Bogomolny equations can be obtained via a dimensional reduction of the selfduality 
equations in 4 dimensions, keeping one of the coordinates constant: several remarkable 
properties can hence be understood using tools coming from instanton theory. In particular, 
in [Hit82] Hitchin, adapting the twistor approach, showed that monopoles may be identified 
with certain bundles over P\ subject to some nonsingularity conditions; the solution is then 
determined by an algebraic curve C in 

Moreover, Nahm managed to adapt the ADHM instanton construction to the case of 
monopoles in [Nah82] . The resulting Nahm's equations can be presented in Lax form: this al- 
lows to use methods from integrable systems in this context. In particular an algebraic curve 
can be constructed, the spectral curve, which is again C. 

These two methods, which Hitchin showed to be equivalent in [Hit83], both yield the spec- 
tral curve as fundamental object. The main advantage of this description is that the spectral 
curve is a gauge invariant object which determines the solutions completely, in principle. In 
practice, however, very little is known about explicit solutions. Most of the work in this di- 
rection concerns spectral curves which are, or can be reduced to, elliptic curves. This is due 
to the fact that elliptic curves are indeed well understood, and are subject of a vast literature 
spanning centuries now. Symmetric monopoles constitute a class where some explicit so- 
lutions are known: in general, symmetry simplifies the problem, reducing the solutions to 
ones written again in terms of elliptic functions. 

Ercolani and Sinha provided a powerful insight in the investigation of exact monopole 
solutions [ES89] . They showed how one can solve the Nahm equations in terms of Baker- 
Akhiezer functions, introducing methods from the algebraic-geometric theory of integrable 
systems in the monopole setting (in particular [Kri77]). Their procedure, although concep- 
tually simple, proves rather challenging to implement explicitly, especially for higher charge 
monopoles. This is due to the fact that it requires a rather explicit knowledge of several 
objects on the spectral curve, which is in general a difficult task for higher genus curves. 
Braden and Enolskii [BE06] manage to implement Ercolani-Sinha construction for a charge 
3 monopole with a particular symmetry, hence finding explicit solutions for monopoles in 
this class. In the present work of thesis, we aim to extend the results of [BE06] to a more gen- 
eral class of monopoles, namely charge 3 cyclic monopoles. 
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A charge 3 cyclic monopole has a (genus 4) spectral curve X of equation 



+ awz^ + Pz^ + - /3 = 0, (2) 

with a, /3, 7 s M, subject to some constraints. We focus in particular on solving one set of 
Ercolani-Sinha constraints, involving relations amongst the periods of X: the number of 
cases when these constraints can be solved is fairly limited. 

The cyclic group C3 acts on X, and we have an unbranched covering tt : X ^ X ofa. genus 2 
curve X. Results by Fay [Fay73] and Accola [Acc71] allow us to express several objects on the 
genus 4 curve in terms of those on the genus 2 curve. For this simplification to take place, the 
first step is the explicit construction of a basis of Hi {X, Z) with particular symmetry proper- 
ties: this proves to be nontrivial, and constitutes one of the main results of this thesis. 
Using the theory of unbranched covers, we manage to reduce the Ercolani-Sinha constraints 
to hyperelliptic ones, hence making them more manageable. Moreover, the genus 4 theta 
function is expressed as a product of hyperelliptic ones, which means that the solutions for 
this class of monopoles are expressible in terms of hyperelliptic functions. 



In order to solve the Ercolani-Sinha constraints, we make use of the Arithmetic-Geometric 
Mean (AGM) method. This method was developed by Gauss [Gau99] to find elliptic integrals 
in a more efficient way, and modified by Richelot and Humbert for the case of genus 2 sur- 
faces with real branchpoints in [Ric36, HumOl]. We manage to extend it to genus 2 curves 
with complex conjugate branchpoints: we give all the details for the subclass of such curves 
to which the quotient monopole belongs; we remark that the methods we use make possible 
such an extension also to generic genus 2 curves with complex conjugate branchpoints. 



Using this extended AGM method, we solve the Ercolani-Sinha constraints numerically, 
via an iteration starting from the known solutions of the tetrahedral monopoles. We find a 
one parameter family of charge 3 cyclic monopoles, by giving explicit expressions for those 
parameters a, /3, 7 in (2) satisfying the Ercolani-Sinha constrants: this provides a quantita- 
tive confirmation of results obtained by Hitchin, Manton and Murray [HMM95] , and Sutcliffe 
[Sut97] obtained by geodesic monopole scattering arguments. 

Finally, a study of the Igusa invariants for the curve (2) allows us to conclude that 4 of the 
above solutions correspond to curves admitting elliptic subcovers: hence, in the class of 
cyclic charge 3 monopoles there are 4 elliptic ones. 



Plan of the work 

Chapter 1 gives some background about monopoles. After a short overview of the origi- 
nal differential geometric formulation, we briefly present the Nahm and Hitchin descrip- 
tions; within this setting we introduce symmetric monopoles, which are a main theme of 
this work. We then describe the Ercolani-Sinha formulation, which, using tools from inte- 
grable systems, provides a suitable framework for exact monopole solutions: in this context, 
we examine in some detail a certain symmetric charge 3 monopole studied by Braden and 
Enolskii. 

In Chapter 2 we recall some elements from the algebraic-geometric theory of integrable 
systems. In particular, we focus on the Lax formulation and, within this framework, we dis- 
cuss the spectral curve associated to a system in Lax form and how the solution can be ex- 
pressed in terms of theta functions on this curve. This introduces some of the techniques 
that are at the basis of the study of monopoles introduced in Chapter 1, upon which the 
present work of thesis is based. 

We study the spectral curve for the cyclic charge 3 monopole in detail in Chapter 3. We 
examine its general properties, and in particular the C3 symmetry, and its quotient with re- 
spect to C3, making use of Fay's theory of unbranched covers. One of the main results of 
this work is indeed the explicit construction of a particularly symmetric basis which satisfies 
Fay's theorem. This allows to express several objects on X in terms of objects on X: in par- 
ticular, the Ercolani-Sinha constraints reduce to hyperelliptic constraints, for which we are 
able to find a solution. 

Chapter 4 is dedicated to the study of the invariants of the curve X, to investigate for 
which values of the parameters the curve X admits a reduction to an elliptic curve. 

Finally, solutions for the Ercolani-Sinha constraints reduced to X are calculated explicitly 
in Chapter 5, using the Arithmetic-Geometric Mean method. After examining the original 
method, due to Gauss, and its adaptation by Richelot and Humbert to genus 2 surfaces with 
real branchpoints, we present a further extension to curves with complex conjugate branch- 
points. Using this extended AGM method, we manage to solve the Ercolani-Sinha constraints 
iteratively, finding a curve of solutions in the moduli space. From the analysis of the curve in- 
variants performed in Chapter 4, we also deduce the existence of 4 elliptic monopoles within 
the class of those with cyclic symmetry and charge 3. 
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CHAPTER 

1 

Monopoles 



We give here some background about monopoles. After a short overview of the 
original differential geometric formulation, we briefly present the Nahm and 
Hitchin descriptions, which are more algebraic-geometric in flavour; within 
this setting we introduce symmetric monopoles, which are a main theme of 
this work. We then describe the Ercolani-Sinha formulation, which, using tools 
from integrable systems, provides a suitable framework for exact monopole so- 
lutions: in this context, we examine in some detail a certain symmetric charge 
3 monopole studied byBraden and Enolskii. 



1.1 A very short introduction to monopoles 

The setting is a gauge theory, the prototype of which is electromagnetism. 
In the geometric formulation of electromagnetism, the fundamental objects are a 1 -forms 
on M^, the electric field E, and a 2-form on R^, the magnetic field B. Out of these, one can 
build the Maxwell field tensor F = B + cdiAE e n'^{Mi), where AI4 = x M is the 
four dimensional Minkowski space-time. In these geometric terms. Maxwell equations (in 
absence of sources) can be written as follows: 

dF = 0, d*F = 0, (1.1.1) 

where * is the Hodge star. The first of eq. (1.1.1) implies that F is closed on M4, hence can 
be expressed as F = dA, where A is a 1-form, the vector potential. 

From its very definition, the vector potential is defined up to an exact form, or in other words 
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it is invariant under a gauge transformation A — > A + da. Geometrically this means that A 
is a connection form of a (trivial) U{1) bundle over M4, and F is its curvature. Note that 
the setting where the necessity of the bundle description really arises is, in fact, that of the 
so-called Dirac monopole [DirSl]. A Dirac monopole is a solution of the Maxwell equation 
with a magnetic charge at the origin of M'^. The magnetic field B corresponding to such a 
magnetic monopole is given by 

B — - — k{xii1x2 a dx3 + cyclic). 

In terms of the "usual" vector B this would read 

B(x) 



where x — (a;i, X2, 2:3), r — |x|. 

With this magnetic field, there exists no globally defined form A such that F = dA. However, 
it is possible to find such a form locally, for instance 

1 1-0086* 

Atv = -; . „ 60 on Un 

47rr sm 

1 I + COS6I 
As = -- ^-^60 on Us, 

where Un = Sn and Us ^ Ss ^ K, where Sn (resp. Ss) is the sphere minus the north 
(resp. south) pole. Thus F = dA holds locally. Note that the following relation holds: 

Ajv — As + da, where a — — — 4>- 

which is precisely the Maurer-Cartan equation for a connection on a U{1) bundle. 
In general, the connection viewpoint is extremely powerful when the space has a nontrivial 
topology, as in this case, where there is a singularity at the origin. However, in this work we 
are mostly interested in monopoles without singularities, defined in the whole of M^: this 
means that for a connection there is always a globally defined gauge potential. Nevertheless, 
these monopoles share important similarities with the Dirac monopole, including some sort 
of topological nature. 



1.1.1 SU{2) monopoles 

Pure Yang-MiUs theory is a gauge theory directly modelled on electromagnetism, but with 
gauge group a non-abelian Lie group instead of U{1). The connection A is interpreted (in the 
quantum theory) as giving rise to a vector particle for every generator of the Lie algebra of G, 
with zero mass; in order to incorporate mass in a gauge invariant way, one has to include in 
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the theory a new field 0, the Higgs field. 
The action is then (proportional to) 

S = [ ({F,F) + {V(b,T>(j)) + Viq^)dvoi) (1.1.2) 

where F is a gauge invariant potential, 2? is the covariant derivative associated to the con- 
nection A, and (F, F) = Tr(F A ★F). 

In the following, we only consider a theory on equipped with the Euclidean metric, with 
gauge group SU{2), a Higgs fields taking values in the adjoint bundle to the principal SU{2) 
bundle on and zero potential. 
The action then reads 

S = [ ((F,F) + (I?0,I?0)); 

since S is interpreted as the energy of the system, one has to ensure that it is finite. This is 
achieved by imposing the following field behaviour at infinity^ (for details see [JT80]): 

l + 0(r-i) 

= 0(r-2), (1.1.3) 



\\m\ = o{r-'), 

where ^ denotes the angular derivative (along a sphere). 

Integrating the action density over a ball of radius R centered at the origin, one has: 

/ ((F,F) + (I>0,I>0)) = / ((F-*I>0,F-*I?0) + 2(*I>0,F)). (1.1.4) 
Upon using the Bianchi identity for F: 

d(0,F) = (*X>0,F), 

which implies that 

/ (*I?0,F) - / (0,F), 

where S\ = Let the connection A be defined on a rank 2 vector bundle E. For large 

i?, as 1 1^1 1 -> 1 by (1.1.3), its eigenspaces define complex line bundles L and L* over 5*^, with 
Chern class ±k. With the decay conditions | | = 0(r^^), the projection of the curvature 



' Note that this ensures the finiteness of the energy only when the fields satisfy the Bogomolny equations 1.1.5. 
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F on these line bundles approaches the curvature of the bundles themselves: 

lim / ((/),F) = ±47rfc, 



hence 



S= (F - *I>0, F - ★l?^) ±47r/c. 
Therefore, if fc > 0, the action S is minimised if (A, 0) satisfies the Bogomolny equations: 

F = (1.1.5) 

A Higgs field solving these equations can in fact be shown to have (in a suitable gauge) the 
following asymptotic expansion: 



I \ k li 



0(r 



-2\ 



This adds a further constraint to the asymptotic conditions (1.1.3); hence we can now give 
the following 

Definition 1.1. An SU(2) monopole of charge k is a solution of the above Bogomolny equa- 
tions with (minimal) energy Ank and boundary conditions 

^=0(,-2) (1.1.6) 

Remark 1.1.1. We notice that monopoles also arise as dimensional reduction of a Yang-Mills 
theory in the four dimensional Euclidean space. More specifically, consider an SU{2) gauge 
theory on the Euclidean W^, with Higgs field taking values in the adjoint bundle and action 
as in eq. (1.1.2), with a non zero potential: 

S= f ((F, F) + (!></., + A(l - l^n^d^x) (1.1.7) 

If we consider only the solutions which are invariant by translations along, for instance, the 
fourth coordinate, these can be dimensionally reduced to R^, where the first three compo- 
nents of A can be interpreted as components of a connection on M^, and the fourth as the 
Higgs field. 

Taking A = in the potential, and keeping the boundary conditions given by the finiteness of 
the action (1.1.7) results in the dimensional reduction of the self-duality equations F = ★F 
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to the Bogomolny equations (1.1.5). 

Thus, SU{2) monopoles can be interpreted as static self-dual solutions of Yang-Mills equa- 
tions on M^. 

Remark 1.1.2. Note that setting A = in (1.1.7) is known (mainly in the physics literature) as 
taking the BPS limit of the Yang-Mills action. Indeed, 't Hooft and Polyakov found the first 
monopole solution in a non abelian gauge theory, namely a static solution of the fields equa- 
tion coming from eq. (1.1.7) with A 7^ 0; their analysis was mainly numerical. Shortly after 
this, Prasad and Sommerfield found an analytic form for the solution in the case A = 0, which 
is a spherically symmetric monopole of charge 1. Then Bogomolny studied further this A = 
limit, also introducing the Bogomolny equations (1.1.5). Hence this class of monopoles is of- 
ten referred to as BPS monopoles: in this work, we refer to them simply as monopoles. 

Remark 1.1.3. Finally, we point out that in the case of a gauge group U{1), Bogomolny equa- 
tions reduce to 

B = grad (f> 

where (j)—-^- This corresponds indeed to the Dirac monopole. 

1.2 Nahm construction 

The above description of monopoles, whilst very beautiful, does not prove very useful 
when finding explicit monopoles solutions. Hence, several different approaches have been 
developed to this aim, both approximate and exact. In this work we are mainly con- 
cerned with Nahm's approach [Nah82] , and the subsequent formulations given by Hitchin 
([Hit83, Hit82]), and Ercolani and Sinha ([ES89]); in this section, we describe Nahm's formu- 
lation. 

In [Nah82] Nahm managed to adapt the Atiyah-Hitchin-Drinfeld-Manin (ADHM, see 
[ADHM78]) construction for instantons to the case of monopoles. The equivalence between 
the geometric description of monopoles and Nahm's algebraic description was proven by 
Hitchin in [Hit83] , where he also showed the equivalence of yet another formulation, that we 
refer to as Hitchin's formulation (explained in section 1.3). 

The crucial construction needed to understand this alternative formulation is the Nahm 
transform, which is a two way transform: it takes solutions of the Bogomolny equations 
(with boundary conditions (1.1.6)) to solutions of the Nahm's equations, a first order differ- 
ential equation between matrices satisfying certain conditions (see end of this section), and 
vice versa. 
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In analogy to the ADHM construction, Nahm considers a Dirac operator on R^, given by: 

\ / a- ■ V ~ {(t) + iz) 

Dl j y cr ■ -D + {(t) + iz) 

where cr = (cti, (T2, 0-3) are the Pauli matrices, p acts on £^-normalizable x valued'^ 
functions on M^, and D'^ is the adjoint of D. 
After using Bogomolny equations, one finds: 

D,Dl = I? • I? - ((^ + + iz), 

which means D^Dl is a positive operator, and hence has no (£^-normalizable) zero modes; 
this in particular implies that the £^-kernel of Dl is empty. One can show* that has k- 
dimensional £^ -kernel A'^z only if z e (-1, 1) • 

Consider then an orthonormal basis for A^^, namely k functions Tp^ such that 



and satisfying 

Jr3 oz 

These functions can be used to define the three Nahm matrices as follows: 

T'J{z)= [ {i^\XaV). (1.2.1) 



These matrices satisfy a number of properties, given below. 
Nl Nahm equations: 

dt 2 



N2 Ti is analytic for z e (-1, 1) and has simple poles at 2; = -1, z ^ I, where the residues 
form irreducible fc- dimensional representations of su(2). 

N3 T,{z) = -T^iz) r,(z)=rf(-z). 

A proof of these properties can be found either in Nahm's paper [Nah82] , or in the review 
[WY07]. 



A solution of eqs. (1.2.2) satisfying conditions Nl, N2, N3 can be mapped back to a con- 
nection and a Higgs field via the inverse Nahm transform, which has a similar structure to the 

^ Both C^s are representations of su(2): the first copy is the spin space where the 0-^ act, the other is the space 
where the connection and the Higgs live. 

^Using an index theorem it follows that Dz has an index fc if 2 e (—1, 1) , and zero otherwise. 
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map introduced earlier, and is given schematically as follows (for details see again [Nah82] 
or [WY07]). 

Consider the differential operator 

-Da; "T" - T^cr' - Xjaj. 
dz 

It can be shown that its £^-kernel is 2-dimensional (this makes use of N2). Choosing an 
orthonormal basis* (wi, ^2) for Er^, define a connection and a Higgs field by 

(t>ab^i / zvl{z)vb{z)dz, 

I (1.2.3) 

/ d 
zvl{z) — Vb{z)dz. 

A and cj) so defined are indeed smooth solutions of the Bogomolny equations with the ap- 
propriate boundary conditions for a charge k monopole. 
We can then state the following theorem. 

Theorem 1.2 {Nahm [Nah82], Hitchin [Hit83J). An SU{2) monopole of charge k, namely 
a solution of the Bogomolny equations with (minimal) energy Ank and boundary conditions 
(1.1.6), is equivalent to a solution of the Nahm equations Nl, subject to the conditions N2, N3. 



Remark 1.2.1. The Nahm transform and its inverse are isometrics, and the Nahm transform 
followed by its inverse gives back the original monopole. 

Remark 1.2.2. Nahm's equations can be cast in Lax form, which allows to use methods from 

integrable systems; we explain this in some more detail. 

Define 

A_i = (Ti + iT2), Ao = Ai^Ti- iT2, 

and hence, introducing a spectral parameter C, eP^: 



A ^ A^iC^ + Ao + AiC, M^-Ao + AiC. 



ThenNl becomes 

dz 

i.e. a Lax form for Nahm's equations. This permits to use integrable system methods to find 
a solution to Nahm's equations. Indeed, an algebraic curve appears, namely the spectral 



^Note that the ambiguity in the choice of this basis is reflected in an ambiguity in the definition (1.2.3) of the 
fields up to a gauge transformation 
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curve, given by the equation (see Chapter 2) 

P(C, il) = det(7yl -A)= dct (r; + (Ti + 1X2) - 2iTzC + {Ti - iT2)C') = 0, (1.2.4) 

and Nahm's equations describe a linear flow on the Jacobian on this curve. As described in 
Chapter 2, there are (almost) algorithmic approaches to deal with this class of systems: in 
particular, in this work we make use of the finite-gap integration method, due to Krichever 
(see section 2.3.1 for more details), in a framework due to Ercolani and Sinha (see [ES89] and 
section 1.5). 

Remark 1.2.3. Finally, it has been possible to solve Nahm's equations in a number of non- 
trivial cases, and then to carry out the inverse Nahm transform, at least numerically: in this 
way a number of monopole solutions of different charges have been constructed explicitly; 
this has proven useful especially to describe qualitative features of monopoles, such as the 
energy density distribution. 



1.3 Hitchindata 

Another approach to monopoles in terms of a certain class of holomorphic bundles on TP^ 
is given by Hitchin; here we only sketch very briefly this construction, referring to [Hit82] for 
details. 

The idea is to consider the set of oriented geodesies on with the Euclidean metric, 
which is in fact isomorphic to TP^; a solution to the Bogomolny equations gives rise to a 
holomorphic bundle on this surface. 

To set notation, introduce homogeneous coordinates [Co : Ci] on P^ = [/o U C/i, with Ui 
standard open sets. Let C = Co/Ci be a coordinate on U^; hence standard coordinates (C, 77) 
on7r^^([/o) c TP^ are defined by (C, ??) — > (and similarly for t/i). 

With respect to these coordinates and the open covering of TP^ given by Ui — TT^^{Ui), the 
line bundle U is defined by the transition function exp ( — j | ) in the intersection C/q n f/i . 
Recall that in P^ the transition function goi = C" 011 f^o n Ui defines the line bundle 0[n) 
(which is the unique line bundle of degree n on P^); denote the pullback of this line bundle 
to TP^ also by ©(n). Take then the tensor product of this line bundle 0[n) with the line bun- 
dle U defined earlier, to obtain the bundle L^{n) U ® 0{n). 
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Consider the differential operator: 

where 7?^ is the covariant derivative along the oriented straight line 7. 
Define a bundle E on the space of geodesies by associating to each 7 the (2-dimensional) 
kernel E-y of D~^. If A and satisfy the Bogomolny equations, one can endow E with a holo- 
morphic structure (with some additional properties, see Theorem 4.2 in [Hit82]). 
The bundle E has two holomorphic sub-bundles, E^, whose fibers E^ are solutions of 
D^f — which decay at ±00 along the line 7. It can be seen (Theorem 6.3 in [Hit82]) that 
E^ ~ L{-k) and ~ fc). The set of those curves 7 for which i?+ i?:^ forms a curve 
in TP^, called the spectral curve S of the monopole. Since a decaying solution decays expo- 
nentially, the spectral curve is also the set of lines along which there is an £^ -solution. 
One can give a more precise characterisation of the spectral curve. That the quotients of E 
by E^ satisfy 

E/E+~L{k), E/E-~L*{k). 

The curve 5* is defined by the vanishing of the map i?+ — > E/E^ and hence by a section of 
{E+Y (g)E/E- ~ 0{2k). In terms of the standard coordinates (C, v) on TP^, S is then defined 
by an equation of the form: 

P{V: = + ai(C)^'"' + ■ • • + ak-i{C)v + «fe(C), (1-3.1) 

where {C),0 < j < k are polynomial of degree at most 2j in 
This curve satisfies the following properties ([Hit82, Hit83]) 

HO S is compact and has no multiple components. 

HI Reality condition: the curve S is invariant under the standard real structure t on TP^ 
The involution r is defined by reversing the orientation of the lines in R^, and in coor- 
dinates takes the following form 

r:(C,^)^ (-i,-^)- (1.3.2) 

This invariance condition translates to the following relations amongst the coefficients 

at 

«,(C) = (-i)^C'^%(-i). 

H2 is holomorphically trivial on S and L{k — 1) is real. 

The first statement follows from the fact that, on S, E^ and E~ are isomorphic, which 
is equivalent to say that i?+ ® E^ (~ L^). Both these properties are consequences of 
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the boundary conditions (1.1.6). 

H3 H"{S,L''{k - 2)) = for s e (0,2). This is equivalent to the non- singularity of the 
monopole determined by S. 

The above findings can be summarised in the following theorem. 

Theoreml.3 {Hitchin [Hit83]). AnSU{2) monopole of charge k, namely a solution of the Bo- 
gomolny equations with (minimal) energy Auk and boundary conditions (1.1.6), is equivalent 
to a spectral curve of the form (1.3.1), subject to the Hitchin conditions HI, H2, H3. 

Incidentally, note that from the above properties one can deduce, using Riemann- 
Hurwitz formula, that the genus of S is gs = {k - 1)'^. 

If fc = 1, the spectral curve (1.3.1) satisfying HI, H2, H3 takes the form : 

77 = (xi + ix2) - 2x3( - (xi - ix2)C (1.3.3) 

where (xi, X2, xa) = x e and can be interpreted as the centre of the monopole; indeed, 
such a curve is the set of all oriented lines through the point x. This curve is also called a real 
section as it defines a real section of the bundle TP^ — > . 

A fc-monopole also has a well defined centre (as well as a total phase), whose coordinates 
(ci , C2 , C3) enter in the spectral curve as coefficients of the polynomial ai (C) as follows: 

ai(C) = (ci + ^C2) - 2C3C - (ci - ic2)C, (1-3.4) 

Note: Nahm's and Hitchin's approaches are indeed equivalent: this is shown 
in the fundamental paper [Hit83], where a link among these two formula- 
tions, and the original geometric description is given. We remark that both 
in Hitchin's and in Nahm's data an algebraic curve appears, directly as in HI, 
or indirectly as in eq. (1.2.4): this is indeed the same curve in two different 
manifestations. 

1.4 Symmetric Monopoles 

Although the methods outlined above provide a great deal of simplification, in general it 
proves quite involved to find an exact solution for a monopole, and not many solutions have 
been found thus far: a notable exception is the case of symmetric monopoles. Symmetric 
monopoles are fc-monopole solutions which are invariant under various symmetries: this 
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constrains their spectral curves, and in some cases it is quite immediate to infer tiie non- 
existence of monopoles with certain symmetries. It is more difficult to obtain results on the 
existence and on the exact form of the solutions, although some solutions have been found 
already in the early paper [HMM95] . The fundamental paper [HMM95] is the basis for the 
exposition below: we have chosen to present this topic in Hitchin formalism, rather than in 
Nahm formalism because, despite using the latter in most of this thesis, we believe that sym- 
metry is more natural and transparent in the former. 

A general monopole curve has the form given in eq. (1.3.1), subject to conditions HI, 
H2, H3; here we investigate the form of these curves when the monopole is required to be 
invariant with respect to certain symmetry groups of of rotations around the origin. To 
implement this symmetry, the monopole is taken to be centered at the origin, which means, 
from eq. (1.3.4), that ai{Q = 0. But before doing this, a special type of symmetry, namely 
inversion, is discussed. 



1.4.1 Inversion 

The inversion map in is the reflection in the (xi, 2:2) plane, namely: 

As this map reverses orientation, it induces an anti-holomorphic map on the space TP\ 
which in standard coordinates reads 

Notice that 4> is very similar to the real structure t of eq. (1.3.2); in fact 

(j)OT{C,T]) = (-C,f?). 

It follows then that if P(C, 77) = is a monopole curve, and hence invariant under t, then the 
inverted curve is P(-C, v) = 0- 

It is interesting to study the monopoles which are invariant under inversion; in particular, 
their moduli space presents an interesting structure (see [HMM95], section 4). Moreover, 
they form a large class, and are also easy to characterise: their spectral curve is given by 
polynomials P{(, i]) — which are even in C- 
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1.4.2 Rotations 

Let us now recall that an 30(3) rotation in corresponds to a Mobius transformation in 
and hence TP^ In particular, a rotation by an angle 6 around the unit vector (ni, n2, ^s) is 
represented in TP^ as : 

-{b + ia)C + {d-ic)' (^-(^b + ia)( + {d-ic)f 

where 

. . e , . e 6 , 

rii sin 2~ sm 2 ~ ' '^^ 2 ~ ^' 2~ 

Thus, a monopole is invariant under a rotation if its spectral curve P{(, 77) = is invariant 
under the associated Mobius transformation, namely P{C\v') = is the same curve. One 
can then examine the consequences of invariance under a number of uncomplicated sub- 
groups of 50(3). 



Example 1.4.1 [Cyclic and dihedral monopoles). The first simple example is the cyclicgroup 
ofrotations Cn around the 0:3 axis. It is generated by 

27rz 27rz 
C=enC) T] = e n 7]. 

Hence, for P(C, 77) = to be invariant under C„, all terms must have the same total degree, 
mod n; since the leading term is rj'', all terms must have degree k mod n. In particular, for 
-P(C) 7?) = to be invariant under Ck, all terms must have degree 0, mod k. 



The groups C„ are extended to dihedral groups _D„ by adding to the generators a rotation 
by TT around the xi -axis, which corresponds to the transformation: 

C = \, = (1.4.1) 

Let us consider some examples of such symmetric curves, namely charge 3 monopoles 
with either C3 or D3 sjmimetry. A general charge 3 monopole curve takes the form 

+ v{a4C^ + asC^ + a^C^ + aiC^ + ao)+ 

+ + /34C' + ^3C' + 02^ + AC' + /3o) = 0, 

where the coefficients are subject to the following conditions: 

a4 = ao, 03 = 0:1, 0:2 = 02; 

P6 = -P0, 135 = Pi, P4 = -P2, P3 = Pa- 
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Imposing the C3 symmetry, the curve reduces to 

+ a77C' + /3C^ + 7C^ - ^ = 0, (1.4.2) 

where a, 7 e M, and also /3 can be made real by a change of coordinates (rotation about the 
2:3 -axis). 

Imposing on (1.4.2) the extra symmetry (1.4.1), one obtains that for D3 symmetric curve 
7 = 0, hence: 

rf + ar](^ + /3(C^ - 1) = 0. (1.4.3) 

Example 1.4.2 [Platonic monopoles). An interesting class is given by those monopoles 
which are invariant under the Platonic groups; in these cases, Klein's theory of invariant 
polynomials proves very useful in implementing the symmetry requirements on the spectral 
curve. For instance, by symmetry considerations alone one can conclude that no octahedri- 
cally symmetric monopole exists for fc = 2 and fc = 3; for fc = 4, it would have the following 
form: 

7/^ + a(C* + 14C^ + 1) = 0. (1.4.4) 

Analogously, the simplest monopole curve with tetrahedral symmetry is of the form (for 

/c = 3) 

+ iaC(C^ - 1) = 0. (1.4.5) 
Notice that after a rotation, this curve becomes 

77^ + a(C^ + 5V2C^ - 1) = 0, (1.4.6) 

thus exhibiting the C3 symmetry around the 2:3 axis (cf. eq. (1.4.2) ). The simplest curve with 
icosahedral symmetry is (for fc = 6) 

7/6+aC(C'° + llC^-l)=0. (1.4.7) 

We remark that the curves (1.4.5), (1.4.4), (1.4.7) only satisfy the reality conditions HI: for 
them to describe monopoles, they also have to satisfy conditions H2 and H3. One can prove 
(see [HMM95]), that the curves (1.4.5), (1.4.4), are indeed spectral curves for monopoles of 
charge 3 and 4, for suitable values of a; on the other hand, there is no charge 6 monopole 
with icosahedral symmetry. 

Note that the quotients of the curves (1.4.5), (1.4.4), (1.4.7) by the respective symmetry 
groups are elliptic curves: this proved a great advantage in obtaining the above results. In- 
deed, the explicit solutions for the tetrahedral and icosahedral monopoles are expressed in 
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terms of elliptic functions, and the constant a in terms of elliptic periods. 

These last examples are in fact paradigmatic of the general situation for symmetric 
monopoles: the symmetry, and reality conditions impose (sometimes stringent) constraints 
on the spectral curve, but this is not enough to conclude that the constrained curve does 
indeed correspond to a monopole. For this to be the case, conditions H2 and H3 (or N2 
and N3 in Nahm's formulation) need to be satisfied as well: this proves to be a very diffi- 
cult step in general, as it involves relations between integrals on the spectral curve, which 
become harder to find explicitly the higher the genus of the curve. In the above cases, the 
spectral curve admits an elliptic quotient, so the integral reduce to hyperelliptic ones, which 
are manageable; we shall see that the same happens for a certain symmetric 3 -monopole 
(cf. [BE06] and section 1.6). In the case of the cyclic 3-monopole, which is the object of the 
present work, the general quotient is hyperelliptic (cf. section 3.1.2), which is a slightly more 
involved case. 

1.5 Ercolani-Sinha formulation 

Based on Nahm's description of monopoles, Ercolani and Sinha developed a method for find- 
ing solutions of Nahm's equations making use of Krichever's method for solving integrable 
systems in terms of Baker-Akhiezer functions on the spectral curve S. They provided an ex- 
plicit expression for these functions, and hence translated Hitchin's constraints H2, H3 to 
constraints on these functions, and hence to constraint on periods on the curve S. 
We describe the Ercolani-Sinha method in some detail here, as this is the method that we 
apply in this thesis to the case of the cyclic 3-monopole: we follow the expositions given in 
the original paper [ES89], and in the more recent [BE06]. 

Ercolani and Sinha consider the differential operator 



does not take the usual form of an eigenvalue problem, since Ai is dependent on z. One can 
obtain the standard eigenvalue equation 



d 
dz 



+ M{z) = — + -Ao{z)+A,{z)C 



related to the Lax equations; studying its spectral theory, one sees that 





(1.5.1) 
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(1.5.2) 



after performing the following gauge transformation 

Q,(z) - C-\z)A,{z)C{z), if = C(z)*, 

with 

Hence the reduced Nahm equations become: 

^ = [g,g+], (1.5.3) 

where Q = C"^Q-i + Qo + CQi and Q+ = Qq + CQi- 

Note that Qi{z) = Qi{0) = Ai{0), and since Ai is symmetric, we can assume it is diagonal: 

Qi(0) = Ai(0) = diag(pi,...,pfe). 

Comparing with the expression (1.2.4) for P{(, 77), we see that the pj are indeed the roots of 
P(C, rj) near C = 00: 

-^2^ - [[[^ ' Pj j ■ (1-5-4) 

One can apply Krichever's method to the reduced Nahm equations (1.5.3) to find a solu- 
tion for Qq in terms of Baker-Akhiezer functions. This relies on the following fundamental 
theorem by Krichever (more details can be found in section 2.3.1). 

Theorem 1.4 {Krichever, 1977). Let S be a smooth algebraic curve of genus gs with n ^ I 
punctures Pj, j = 1, . . . , n. Then for each set of gs + n — 1 points . . . , (5gs+„-i in gen- 
eral position, there exists a unique function {t, P) and local coordinates Wj{P) for which 
Wj{Pj) = 0, such that 

1. the function ofP e 5* is meromorphic outside the punctures and has at most simple 
poles at 5 s (if all of them are distinct); 

2. in a neighbourhood of the puncture Pi the function takes the form 

(z, P) - (^ji + oc% (2) , wi^wi (P) , m e N+ . 

We can use this result to express $ in eq. (1.5.1) as follows. Let be the columns of $ in 
eq. (1.5.1), normalised so that 

exp(CAi(0)z)$|oo =M„ (1.5.5) 

Consider the k points Pj above C = 00 as punctures; we denote their rj coordinate by oo^ . 
Then one has the following result: 
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Theorem 1.5 {Ercolani-Sinha [ES89] ). The j-th component of^i, normalised by (1.5.5) is 
given by Theorem 1.4 on S, where the punctures are the k points Pj above ( = oo. 
From this, the matrix Qq can be reconstructed from 

{Qo)ji ^ -{Pj - Pi)aji = -{pj - pi) lim Cexp(Cpiz)*j(2;,P). 

To summarise the complete Ercolani-Sinha method, one can express the components of 
<I> in terms of these Baker-Akhiezer functions, and hence reconstruct Qo! then, using 1.5.2, 
one can determine C(z), and therefore A^i, Ai^, Ai, from which reconstruct the matrices Tj. 

We explain in detail how to obtain $ and Qp in the next subsections. 

1.5.1 Solving Nahm's equations via the spectral transform 
Baker-Akhiezer functions 

The functions 5"; of section 1.4 can be expressed explicitly in terms of theta functions on S; 
in order to do so, we introduce the following basic ingredients: 

1. The differential "fao- From eq. (1.5.4) we see that the rational function 77/^ has the 
following asymptotic behaviour: 

^ PjQ near 00 j . 

Hence its differential satisfies 

d - (^-^ + 0(1)) dt near ooj , (1.5.6) 

where t = 1/C is a choice of local coordinate. Using the above as motivation, one can 
define a second kind differential 700 such that 



7oo (-^ +0(1) nearooj; (1.5.7) 



Pi 

7oo = 0; (1.5.8) 
7oo holomorphic away from 00. (1.5.9) 



2. The vector U. After Theorem 1.4 and the subsequent observation, one can introduce, 
for a particular flow, the vector with components: 




(1.5.10) 
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This differential encodes ttie flow on the Jacobian corresponding to the time evolution 
of the system. 

3. Normalisation factors Vj. To obtain the normalisation required in Theorem 1.4 we 
introduce the following normalisation factors 



Vj — lim 



700 +^ 



Note that these in fact depend on the basepoint Pq, and are finite quantities, because 
of the asymptotic property (1.5.7). 

We can now give an explicit expression for the non normalised function in terms of theta 
functions; we introduce first the following functions: 

e{AiP)-z, + zU) U,^.9iAP)~R^) ^ ^^ ^^^^ 



where A denotes the Abel map, /C the vector of Riemann constants based at Pq (see Appendix 
A for some properties of these objects). Here we have set 

gs+n-l n 

Z^=Zt + A[P,) = A{^j) + JC, A{5,)-Y,A{P,) + K,, 

s=l 3=1 
gs-l Ss-1 

R,= J2^ (Ss) + A (Pj) + /C, Sj=J2-^ i^s) + A (<5g,-i+,) + /C. 

s=l s=l 

Also, by Abel's theorem, Zj is equivalent to an effective divisor A ^ of degree gs- 
From eq. (1.5.11), we obtain the following expression for the normalised functions 



.J - -F^TT^e V^o / . (1.5.12) 



The functions in eq. (1.5.12) have all the required properties of Theorem 1.4: for a fur- 
ther discussion on this point we refer to section 2.3.1, where we examine in more detail how 
to build Baker-Akhiezer functions on a spectral curve in a more general context. 

Using the expression (1.5.12) for the <i>j in Theorem 1.5, we obtain the following expres- 
sion for Qo in terms of Baker-Akhiezer functions: 

(Qo) - . .i.,-.A 0{A{Pi)-Z3 + ^U)e{A{P,)-Z,) 

Wo),i- pi)cjie e{A{P3)-A3 + zU)o{A{Pi)-Zjy ^^-^-^^^ 



where 



c,i = lim CgAP)^ P - (C, v) e S. (1.5.14) 

P-i-OOi 
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1.5.2 The Ercolani-Sinha constraints 

As we have now expressed Qo concretely in terms of Baker-Akhiezer functions, we are able to 
impose Hitchin's constraints H2, H3 explicitly on these Baker-Akhiezer functions: these are 
implemented as constraints on the periods of the curve S. We examine this in some detail in 
this section, since in the applications of the Ercolani-Sinha method we are most interested 
in, namely the 3-monopole of [BE06], and the cyclic 3 -monopole studied later in this thesis, a 
large part of the analysis has been devoted in fact to the implementation of these constraints. 



Let us consider then Hitchin's constraints in detail. 



i L{k - 1) is real. [Har78]. 



ii Triviality ofL?. is trivial if and only if it admits a (nowhere vanishing) holomor- 
phic section /, namely, in local coordinates, if and only if one can find holomorphic 
functions /o on C/q n 5* and /i on Ui n S, such that, on C/o n C/i n 5 one has 

/o-expf-2^Vi- 



c 



Taking the logarithmic derivative we get 



dlog/o = -d i^Yj +dlog/i. (1.5.15) 

In order to avoid any essential singularity of /o on Uo, then d log /i has to cancel those 
of -d {2r]/C) at oo; this means, using eq. (1.5.6), that 

dlog/i - 2 (-^ + 0(1)) ne&vooj. (1.5.16) 

Then, using eq. (1.5.15) we can define 

m,= £ dlog/o ^-^iii/l^S/i, 
n,= 2^ £ dlog/o^ 2iii/^°S/i. 

Comparing eqs. (1.5.6) and (1.5.16), we can write: 

1 " 

7oo = -dlog/i +i7r^mj, 

1=1 

where the Uj are the canonically a-normalised holomorphic differentials on S, and the 
second addend ensures that the normalisation (1.5.8) is achieved. 
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Integrating 700 around b- cycles yields the Ercolani-Sinha constraints 

f Joo ^ iT^Tij + in'^Tjimi, (1.5.17) 
•'^j 1=1 

where t is the period matrix of S; these constraints are equivalent to the triviality of 
on S. This condition implies that the vector U of eq. (1.5.10) (which appears in the 
Baker-Akhiezer functions (1.5.12)) takes the following form: 

U^^n+Km (1.5.18) 

or, in other words, J7 is a half-period. 

Note that Ercolani and Sinha in [ES89] take n to be zero: Braden and Enolskii in [BE06] 
show that this need not be the case. 



iii H3: H'^ [S , {k — 2)) = 0. This condition can be restated in terms of the theta divisor as 

L\k-2)e JS^-\Q for (-1,1). (1.5.19) 
This constraint must be checked using explicit knowledge of the 6-divisor. 
Summarising, we have the following theorem: 
Theorem 1.6 {Ercolani-Sinha Constraints [ES89J). The following are equivalent: 
1. is trivial on C. 



We mention here another equivalent characterisation given by Houghton, Manton and 
Ramao, namely: 

Theorem 1.7 {Houghton, Manton and Romao [HMROOJ). The following condition is equiv- 
alent to the conditions of Theorem 1.6: there exists a 1 -cycle c = n • a + m- fa such that for every 
holomorphic differential ri, 

c 

Moreover, Braden and Enolskii show a further property of this cycle in [BE06] : 
Corollary 1 .8. The cycle c satisfies 

T*C = -c. 
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1.5.3 Braden-Enolskii extensions to the Ercolani-Sinha theory 

Braden and Enolskii in [BE06] propose an extension to the Ercolani-Sinha formulation, also 
giving simpler expressions for and Qo (see section 3 of [BE06]]. 

We begin by providing an expression for in terms of theta functions, alternative to eq. 
(1.5.12): 

(a{P) - A{(x,,) + zU - 6^.^ (~JC) z fy^-zu, 

^,{z,P)^g,{P) \ L_ Vy ^ P„ (1_5_21) 

(AiP)-Ai^,)-IC)0^^^ (z[/-/Cj 

Here 9w,^^ are theta functions with characteristics, A is the Abel map, z e ( - 1 , 1 ) , and P e S. 
The vector IC is defined by 

}C = IC + A ^(n - 2) ^ oofc^ , (1.5.22) 

where IC is, as above, the vector of Riemann constants. 
In [BE06] it is established that 

Lemma 1.9. The following equivalence holds: 

H°{S, L'{k -2))^0-^ e{sU -JC)^Q (1.5.23) 

fors G (0, 2), withJC as in (1.5.22). 

Thus the problem is to determine when the (real) line sU — K. intersects the theta divisor 

e. 

We remark that IC satisfies the following properties: 

1 . /C is independent of the choice of base point of the Abel map; 

2. e{jC) = 0; 

3. 2JC G A; 

4. for n > 3 we have IC e Osinguiar- 

The point IC is the distinguished point Hitchin uses to identify degree 5 — 1 line bundles with 
J{C). The proof of these properties together with the following lemma further constraining 
the Ercolani-Sinha vector may be found in [BE06] : 

Lemma 1.10. U ± K. isa non-singular even theta characteristic. 
Moreover, Braden and Enolskii give another expression for Qo{z): 
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Theorem 1.11 [BradenandEnolskii [BE06]). The matrix Qo{z) (which has poles of first order 
atz — ±1) may be written 

O (z) , = e > ~ ^' MA(^i)-A(^.)) ^(Aooi) - Aooj) + {z + 1]^ - ^.fa-.,) 
^'f(oo„oo,) 9{[z + l]U-^) 

(1.5.24) 

Here E{P, Q) = £{P, Q) / dx{P)dx{Q) is the Schottky-Klein prime form, U - K. = \p+ ^rq 
(p, q ^ "L^) is a non-singular even theta characteristic, ande^i — eij — ±1 is determined (for 
j < I) bye.ji = <^jj+i^j+ij+2 ■ ■ ■ Then — 1 signse^j+i = ±1 are arbitrary. 

For the proof of this theorem, we refer again to the paper [BE06] . 

1.6 A certain symmetric charge 3 monopole 

We explore here in some detail a class of charge three monopoles with a certain symmetry, 
introduced by Braden and Enolskii in [BE06] : this constitutes a particularly important ex- 
ample, because it is one of the few cases (and the only one of charge greater than 2) where 
the Ercolani-Sinha method has been fully implemented. Moreover, it is particularly rele- 
vant for this thesis as we aim to generalise this to a slightly less symmetric case, the cyclic 
3-monopole; the work we present in the later chapters is, in fact, partly based on this. 

1.6.1 The spectral curve, a homology basis and the Riemann period ma- 
trix. 

Braden and Enolskii in [BE06] consider the following spectral curve of a charge 3 monopole 
(cf. eq. (1.4.3)) 

ry^ + ^C^ + 7C^ - /3 = 0. (1.6.1) 

The corresponding Riemann surface, which we denote'^ by Xq, has genus 4; viewing it as 
a branched three sheeted cover of the Riemann sphere P^, Xq has 6 branchpoints, \, i ^ 
1...6: 

A2 = pAi, A3 = p^Ai, A5 = /C1A4, Ae = p^A4, (1.6.2) 

where A = 27(12^^ + 7^). Hence the curve can be parametrised in the following way, which 
we use later 

6 

= - Ai), where w = -/J^^t?, z = C- (1.6.3) 

i=l 

^This choice of notation is clarified in sections 3.1, 3.2. 
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We also mention here the basis of the holomorphic differentials used in [BE06] : 

dz dz z dz z^dz 

w 

This curve (1.6.1) admits the following automorphisms: 

s : {z,w) — !• (z, pw), (1.6.4) 
TZ : {z,w) — > [pz, w), (1.6.5) 

r: {z,w) ^ (^~^,-^^ , (1.6.6) 

0: {z,w)^ {-l'-^) ' 

where p — exp(2i7r/3); s corresponds to shifting a point up one sheet, TZ rotates by 27r/3 with- 
out shifting sheet, t is the real involution of eq. (1.3.2), and </> is an inversion. 

The first task of [BE06] is providing a canonical basis for Hi{Xo, Z): this is given here in 
Figure 1.1. Also, we give an expansion for these cycles in terms of "basic arcs" as follows. 
Denote by 7^ (i, j) the arc going from branchpoint to branchpoint Xj on sheet k: 

7fe(i,j) = arcfc(Ai, A-,), iy^j^l,...,6 on fc-th sheet 

Then: 

ai =71(1, 2) +72(2,1), bi -71(2, 1)+ 73(1, 2), 
a2 -72(3, 4) +73(4, 3), b2 =72 (4, 3) +71 (3, 4), 

as = 73(5, 6) +71(6, 5), b3 = 73(6, 5) +72(5, 6), (1.6.8) 

ao = 73(1, 2) + 71 (2, 6) + 73(6, 5) + 72(6, 1), 

bo = 73(1, 2) + 71 (2, 5) + 72(5, 6) + 73(6, 3) + 71 (3, 4) + 72(4, 1) 

We remark that here the monodromy (calculated with respect to infinity) is [1, 2, 3] (more 
details on monodromy and how to calculate it are given in section 3.3). 
This basis has the following property under the symmetries s and TZ: 

TZtia^) - a,+fe, 7^^(b,) = b,+fe, ; (1.6.9) 

s*(bo)=ao, s*(bfc) = afc,n = 1,2,3 

where the last map has the effect of just shifting up (in our conventions) one sheet. 



The choice of such a symmetric homology basis allows us to relate the matrices of a and 
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Figure 1.1: Symmetric iiomology basis of [BE06] 

b-periods, and hence simplifies greatly the Riemann period matrix. 
Introducing vectors 

b= (61,62,63,64)^ = 



Vl, . . . , ® Vl 

ai J a4 



V2, • • • , $ V2 

ai J a4 



C = (C1,C2,C3,C4) = <?> V3,...,* V3 



ai Ja4 



d= (di,d2,rf3,d4) = [ f V4,. 



ai J an 



T 

V4 I , 



one then finds the following result: 

Proposition 1.12 [Welhtein [Wel99]; Matsumoto [MatOl]; Braden-Enolskii, [BE06]). For 

the curve Xq of equation (1.6.1), one has 



A = \ (t)Vi \ = (a;,b, c, d) 

afc / i,fe=l,...,4 



(1.6.10) 



B 



V,; I = {pHx, p^Hb, p^Hc, p^Hd) = HAA, 

i,/c=l,...,4 



whereH = diag(l, 1, 1, -1) andk = diag(p, p^, p^, p^). Then the Riemann period matrix is of 
the form 



n = AB-' = p[H-{l-p) 



XX 



x^Hx } 



(1.6.11) 



andTb is positive definite if and only if x^ Hx < 0. 
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1.6.2 Solving the Ercolani-Sinha constraints 

For convenience we rewrite tlie curve of eq. (1.6.1) as follows: 

= + bz^-l = {z^ - a^) (^z^ - (1.6.12) 

where 5 = anda'^ = + yp~T4). 
The Ercolani-Sinha constraints here read: 

A + B = 1^(1,0,0,0). (1.6.13) 

where^ v = Qj3i (in general v depends on normalizations). 
These constraints may be solved using the following result. 

Proposition 1.13. The Ercolani-Sinha constraints (1.6.13) are satisfied for the curve (1.6.1) if 
and only if 

X = C{Hn + p^m), (1.6.14) 



where ^ 

V 6/33 

^ ^"1^ g 15) 

[n.Hn — m.n + m.Hm] [n.Hn — m.n + m.Hm] 

The proof, which follows from Proposition 1.12, can be found in [BE06]. 
Thus, solving the Ercolani-Sinha constraints amounts to imposing the four constraints 
(1.6.14) on the periods Xk. For the curve (1.6.1) one can actually calculate these integrals 
in terms of special functions (Gauss hypergeometric functions) and the constraints in fact 
reduce to a number theoretic one. 



We remark that in this case all the integrals reduce to just 8, namely: 

dui^Ii, Jdu, = Ji, i^l,...,4, 



where Ai = (d,0),A4 = {f3,0)and(3 = 4, and it is intended that these integrals are computed 
on the first sheet. This follows from the symmetries of the curve, for details see [BE06] , p. 50 



^ To see this observe that (1.7) requires that -25ifc = /„.„+^.u fi'*"' for the differentials r2(i) = = — 

nam. _ nrj 



= — — In the parametrisation (1.6.3) we are using we have that 



ap 

drj 



: I — = I — = _3/33 / qW, 
Jai w Jui _/3- 3 n Jdi 



I at —13 3 r) 

Imposing — 2 = f^^^' = — (n.x + pm.H.x) and so n.a; + pm.H.x = v, with the value of v stated. 
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(also cf. section 3.10). In particular, usingli and J\, which in terms of hyperelliptic function 
read 



dz 1'K\fZa ^ / 1 1 



Xi(5)= / - = —,FA-^,-^^,-a 







(1.6.16) 







the x-periods are 



xi =-(2J'i+Ii)p-2Ii- J-i, X2 = (Ji + +2Ji, 



(1.6.17) 



where we remark that X2 = pxi and 2:3 — p^xi by symmetry. 

Hence, using (1.6.14) and (1.6.17) the Ercolani-Sinha constraints can be rewritten as 

X, = ^{eiu, + p^rui) = (<5,Ii + PiJi) + (7,Ii + 5ai)p. (1.6.18) 

where e, = 1 for z — 1, 2, 3 and 64 = —1. Now this expression can be solved for n^, to,, as in the 
following 

Proposition 1.14. For each pair of relatively prime integers m, n such that 

(m + n){m — 2n) < 

we obtain a solution n = (^n m - n -to 2n-nij,m = (^m -n n-m Snj to the 
Ercolani-Sinha constraints for the curve (1.6.12) as follows. First we solve the equation 



2n-m _ 2-^1(5, |; 1,0 
m + n " 2^1(5, 

fort; we have then 

^ _ 1 - 2i ^ _ -b+VW+A 



(1.6.19) 



and we obtain jS from 



=-{n + m)^ « ^F,{\,\-l,t) (1.6.20) 
3V3 (1 + a6)3 3 3 



witha^ = t/{l - t). 



We remark that the constraint gcd(TO, n) = 1 ensures that we have indeed a primitive 
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vector in the lattice Z^. We also point out that eq. (1.6.19) can be rewritten as follows: 

n = = (1.6.21) 

m + n Ji 

The quantity TZ plays an important role in simplifying the period matrix of a curve whose pa- 
rameters satisfy eq. (1.6.19) (see eq. (1.6.26)). 



Hence, the problem of finding solutions to the Ercolani-Sinha constraints has been re- 
duced to that of finding solutions to the transcendental equation (1.6.19): remarkably, this 
can be done. One finds in Ramanujan's second notebook several results of generalised mod- 
ular equations, and various theta functions identity, some of them recently proven (see 
[BBG95], [Ber98]). We report here some of the solutions to the Ercolani-Sinha constraints 
found using these methods by Braden and Enolskii, referring to [BE06] for details. 



n 


m 


(2n — m)/{m + n) 


t 


b 


2 


1 


1 


1 

2 





1 





2 


1 1 5\/3 

2 ^ 18 


-5\/2 


1 


1 


1 

2 


1 5^3 

2 18 


5%/2 


4 


-1 


3 


(63+ 171v^- 18v^)/250 


(44 + 38v^ + 26^)/3 


5 


-2 


4 


1 1 153^3-99^2 

2 250 


9^458 + 187^6 



(1.6.22) 



Hence, we can see that only a discrete (countable) number of solutions appear, meaning that 
there is at most a countable number of monopoles with spectral curve of the form (1.6.1). In 
particular, we point out from the table above that the tetrahedral monopole of eq. (1.4.5) is 
one of them. 



1.6.3 Covers and reduction 

We briefly describe how the curve under consideration covers elliptic curves, which provides 
a noteworthy simplification of the period matrbc (as well as a better understanding of some 
of the results based on Ramanujan's identities, see [BE06]). This is particularly relevant in 
view of our subsequent study of the cyclic 3 -monopole, where again a quotient appears. 

The results about the covers are summarised in the following 

Lemma 1.15. The curve with equation = with arbitrary values of the parameter 

b, is a simultaneous covering of the four elliptic curves £±, £1,2 as indicated in the diagram, 
where C* is an intermediate genus 2 curve 
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^0 = {x, y) 




£+ = {z+,W+) £_={Z_,W_) £i^{zi,Wi) £2^{Z2,W2) 



The equation for Xf) is 

= (/,2 +5)2+4 

The equations of the elliptic curves are 



£±: wi= z±{l- z±){l~k\z±), (1.6.23) 
f 1 : zl + wi^ + 3zi + 6 = 0, (1.6.24) 
£2 ■■ W2^ + Z2^ + fez2 - 1 = 0, (1.6.25) 



where the Jacobi moduli k± are given by 

2 p{pM±l){pMTl? 
± (M±1)(M + 1)3 

withM = f, K = (2i- = (52 + 4)^. 



When a curve covers one of lower genus, the period matrbc admits a reduction, i.e. can be 
expressed in terms of a lower genus period matrix, and similarly, the associated theta func- 
tions can be expressed in terms of lower dimensional theta functions. 
As the curve under consideration has many covers, it admits many reductions too; the case 
where a true simplification occurs is when the Ercolani-Sinha vector reduces as well: this is 
the object of the rest of this section. 



Before doing this, we give an expression for the period matrix of a curve of the form 
(1.6.1), satisfying Ercolani-Sinha constraints. Combining equations (1.6.10), (1.6.11) and 
(1.6.17), one obtains: 
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A = 



B = 



/ -1 - 2p- (2 + p)7^ l + 2p 1 + 2/9 + (1 - p)7e -1 + /9 

2 + p + (l-p)7^ -2-p i-p-{2 + p)n -1 + p 

-1 + P + (1 + 2p)7^ 1-p -2-p + {i + 2p)n -1 + p 

y 3 + 3p-3p7^ -3p-3{i + p)n o 

/ 2 + p + {i-p)n 1-p i-p-{2 + p)n 2 + p^ 

-l+p+{l + 2p)n l + 2p -2 - p + {1 + 2p)TZ 2 + p 

-l-2p-{2 + p)n -2-p l + 2p+{l- p)n 2 + p 

y 3-3(l + p)7^ 3-3p7^ y 



Jl 



I2 

(1.6.26) 



where 7^ is a defined in eq. (1.6.21). 

The Riemann matrix n = [A, BY admits a reduction with respect to each one of its columns: 
we consider here a reduction with respect to the first column. 
Using eq. (1.6.14) it follows that 



HA = 



A 
B 



i l\ 





where M is the 2(7 x 2 integral matrix 



i{iin + p^m) 
i{pn + Hm) 



(1.6.27) 



„ , ni — mi n2 — m2 — —n^ — nii 7712 TO3 —1114^ 

Ar = 



—mi 



-7712 



-TO3 



~TO4 Til 712 1^3 ^4: 



For every two Ercolani-Sinha vectors n, m one has 



M^JM = d 



1 

-1 



d = n.Hn — m.n + m.Hm — ^^(£j"-j ^ '^j'^j + ^j'l^'j)- 

7 = 1 



Theorem 1.16. (Braden and Enolskii [BE06]) For the symmetric monopole one can reduce 
by the first column using the vector (1.6.27) luhose elements are related as in Proposition 1.14, 
with gcd{ni, mi) = 1. Then 

d ~ 2{ni + mi)(mi — 27ii) 
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and ford ^ there exists an element a of the symplectic group Sp2g{'^) such that 

( {p + 2)/d a/d ... \ 



rj, = cr o n, 



a/d 


V 



Letting pmi + qni = 1 then 



a = gcd(mi + 4ni — q [mi — 2ni], ni — 2mi — p [mi — 2ni]). 

When a — 1 a further symplectic transformation allows the simplification t'h — p/d. 
Under a the Ercolani-Sinha vector becomes 



aoU = ao {m^ + n'^n,) = (1/2, 0, 0, 0). 



The proof of Theorem 1.16 is constructive: in particular there is an algorithm by Martens 
[Mar92a, Mar92b] for constructing a which has been implemented in M aple by Braden and 
Enolskii [BE06]. 

The theory of Weiestrass reduction implies that the theta functions reduce correspond- 
ingly, taking the form^: 



m=0 
d-1 

E' 

m=0 



ma p + 2 
d 



f 




{Dw;Dt*D) 



{dz-d{p+2))e{{wi + —,W2,W^)]T*), 



where D = diag((i, 1, 1). 



Using the above expression for the theta function in eq. (1.5.24), we conclude that one 
can express Qq in terms of genus 1 theta functions, or 

Proposition 1.17 {Braden and Enolskii [BE06]). For symmetric monopoles the theta function 
z -dependence ofQo{z) is expressible in terms of elliptic functions. 

This means in particular that the last constraint H3 in Ercolani-Sinha formulation (1.5.19) 
can be reduced to the problem of finding zeros of elliptic functions. Using this, Braden and 



^When gcd(o, d) =^ 1 a smaller multiple than dj = d would suffice here with correspondingly fewer terms in the 
sums < m < di — 1. 
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Enolskii manage to limit the countable number of possible solutions of table (1.6.22), estab- 
lishing a stronger result: 

Theorem 1.18 {Braden and Enolskii [BE09J). The only curves of the family (1.6. 1) that yield 
BPS monopoles correspond to tetrahedrally symmetric monopoles. These have b = ±5\/2 and 
(m,n) = (l,l) or (0,1). 
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CHAPTER 

2 

Algebraic-geometric aspects of integrable 
systems: some basics results 



Integrable systems are a vast area of research, and we do not attempt to give 
even just a short summary here. We only recall some elements on the Lax for- 
mulation, which leads to a description of such systems in terms ofRiemann 
surfaces; this is the point of view we take in this work of thesis. We illustrate 
how to associate a spectral curve to a Lax pair with a parameter, and how this 
allows to realise the dynamical evolution of the system as a linear flow on the 
Jacobian of this curve. 

**Integrable systems, what are they? It's not easy to answer precisely. The question can oc- 
cupy a whole book, or be dismissed as Louis Armstrong is reputed to have done once when 
asked what jazz was - 'If you gotta ask, you'll never know!'" 

N. Hitchin [HSW99] 

The study of integrable systems is a very old subject, which has been of interest to both math- 
ematician and physicists for centuries. Despite this, or perhaps because of this, it is difficult 
to give a general definition that encompasses all the different ways in which integrability 
manifests itself: the main feature of such systems is considered to be, depending on the sit- 
uation, the existence of "many" conserved quantities, or the ability to give explicit solutions, 
or the appearance of algebraic geometry, or the presence of certain topological or algebraic 
structures, etc. . 

Hence, much work has been done in the last century on the subject on integrable systems, in 
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both finite and infinite dimensions, using and indeed developing many sophisticated math- 
ematical techniques, from PDEs to symplectic geometry, from differential topology to rep- 
resentation theory. Here we are mainly concerned with algebraic-geometric methods in in- 
tegrable systems, involving, for instance, algebraic curves, theta functions, Abelian varieties 
etc. . 

We already see the appearance of such objects in several classical dynamical systems: it 
was known since the nineteenth century that solutions for many classical tops, or already 
the simple pendulum, were expressible in terms of elliptic functions. A posteriori, we see 
that there is an algebraic curve in the background. Indeed, modern techniques make this ex- 
plicit: we can build this curve from the beginning, use it to obtain constants of motion, and 
eventually to represent the solutions. This is achieved through the Lax formulation, which 
is the main subject of this chapter. 

Several classical integrable systems admit a Lax description: for instance, in finite dimen- 
sions, all the completely integrable Hamiltonian systems, which constitute a large and im- 
portant class; moreover, several nonlinear equations, e.g. the Korteveg-de Vries equation, the 
Kadomtsev-Petviashvili equation, the nonlinear Schrodinger equation, etc. . In particular, 
the theory of the so called finite gap solutions of KdV type equations led to the more general 
results for solutions to Lax equations in terms of Baker-Akhiezer functions. The main result 
that we discuss in this chapter is in fact how to express solutions of a Lax equation in terms 
of Baker-Akhiezer functions, using a result due to Krichever [Kri77] . We remark, though, that 
the introduction of algebraic-geometric methods in integrable systems has many contrib- 
utors: we refer to [BBE+94] and references therein for more details, and to [BBT03] for an 
introduction. 

This chapter is based on [BBT03, AHH90]. 

2. 1 The Lax formulation 

Within the realm of integrable systems, what we are concerned with in this thesis are systems 
admitting a Lax formulation: 

Definition 2.1. A dynamical system admits a Lax formulation if its time evolution can be 
described by an equation of the form 

^^A = A=[A,B], (2.1.1) 

where A and B are N x N matrices, and [A, B] = AB - BA. The pair {A, B) is called a Lax 
pair. 
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It can be seen that the entries of B are functions of those of A, hence the equation of the 
motion (2.1.1) is nonlinear. 

In general, the Lax formulation for a given system is not unique: from the form of eq. 
(2.1.1) it is immediate to see that if A, B axe a Lax pair for a certain system, then so are 

Ai = and Bi = MBAJ-^ + MM, for any TV x iV matrix M. 

An important feature of such systems is the immediate appearance of conserved quanti- 
ties. Indeed, the solution of eq. (2.1.1) is of the form 

A{t) ^ g{t)AiO)g{t)-' with s.t. ^ = ^5"' 

It follows that all the quantities of the form 

/(fe) = XrA*^ (2.1.2) 

are constants of motion. The presence of "many" constants of motion is in general a fun- 
damental feature of integrable systems; note, though, that from eq. (2.1.2) nothing can be 
inferred about their independence. 

Since the components of the characteristic polynomial for L are functions of these traces, 
the whole spectrum of L is conserved: the time evolution of a system in Lax form is isospec- 
tral. 

2.2 Lax pairs with a parameter 

In many cases, the Lax pair depends on a complex parameter A: 

iA = [AA,BA]. (2.2.1) 

Lax pairs of this form arise quite naturally for many systems, and the study of the analytical 
properties of the Lax equation depending on a parameter gives considerable insight into the 
understanding of the dynamical evolution of the system. A fundamental object for this aim 
is the spectral curve. 

Definition 2.2. Given a Lax pair Ax, Bx depending on a parameter , the characteristic equa- 
tion for Ax 

A) = det(.4A -/il) = (2.2.2) 
defines an algebraic curve C, called the spectral curve. 

37 



2.2. Lax PAIRS WITH A PARAMETER 



In the following we assume it is irreducible, and that it can be completed to obtain a Rie- 
mann surface. 

Notice that, since the coefficients of the characteristic polynomials are functions of the 

of eq.(2.1.2), the spectral curve is preserved by the flow. 
To understand equations of the type (2.2.1), one has to study the spectral curve and its struc- 
ture. In particular, the dynamical information about the system is encoded in a line bundle 
on C, which describes a linear flow on the Jacobian J{C) of C. Under these circumstances, 
together with a restriction on the form of Bx, the solution can be written in terms of theta 
functions. The rest of this chapter explains this in some more detail. 

Firstly, we recall that the genus of the spectral curve can be calculated to be 

N{N - 1) 

k 

where rik is the order of Ax at a pole A^., and the sum runs over all poles. 

The first important object on C is the eigenvector bundle, a line bundle on C built as 
follows. For every values of the pair (/i, A) = P e C we have an eigenvector v{P, t) e V, 
hence an eigenspace: assuming the eigenvalues all have multiplicity 1, this eigenspace has 
dimension 1 at every point on C; this defines the bundle of eigenvectors L^. 
Introducing the eigenvector map 

/t : C ^ Py ; P ^ Cv{P, t) C V, 

we define the eigenvector bundle^ Lt as 

where Pic'*(C), is the group of line bundles of degree d. This degree d of Lt is computed to be 
( e.g. in [BBT03]) 

d = g + N-l. 

Note that the eigenvector bundle does depend on t, since, although the eigenvalues are 
constant with time, the corresponding eigenvectors are not. Recalling that Pic''(C) ~ J{C), 
we have that the time evolution of the dynamical system in Lax form (2.2.1) can be translated 



'Note that this is the dual of the bundle of eigenvectors. 



38 



2.2. Lax PAIRS WITH A PARAMETER 



into time evolution on the Jacobian of its spectral curve: 



(2.2.3) 



It is possible to show the following result (see e.g. [Gri85]). 

Theorem 2.3, Under certain conditions on Bx, the flow (2.2.3) is linear on J{C). 

We recall that the Jacobian of a Riemann surface is a complex torus. Hence, the above re- 
sult means that, for an integrable system in the form (2.2.1), it is possible to realise the linear 
flow on the torus J(C). Note the relation with completely integrable Hamiltonian systems 
here: for such systems, the Arnold-Liouville theorem states that the phase space admits a 
(local) fibration by tori, built such that the dynamical flow is linear on them.The interesting 
thing is the appearance of a torus in both descriptions: the advantage is that the Jacobian 
has a given canonical linear structure, while the Arnold-Liouville tori have a linear structure 
defined by the flow itself. 

The conditions on B\ in Theorem 2.3, for the flow defined by (2.2.1) to be linear, are that 



where P{z, w) is an arbitrary complex polynomial of two variables, and ( , )+ denotes taking 
the regular part of P, with respect to its Laurent series expansion in A. 

Note that here, given a matrbc equation of the form (2.2.1), we have built a linear flow of 
line bundles Lt. It is possible to do the converse, namely, given a linear flow of line bundles 
onC, of degree 5 + — 1 to construct a solution (t) to eq. (2.2.1) explicitly in terms of theta 
functions on J{C). This is done as follows. 

Take a linear flow of line bundles of degree g + A^ — 1 on C. Take an "initial value" line 
bundle Lq. represented by a positive divisor 5, or = L{5). 

Denote by • ■ ■ , V'^ a basis for i/°(C, Lq). We remark that Lq has precisely A^ sections when^ 



Bx = {P{Ax,\-^))+ 



(2.2.4) 



<limH^{C,0{Lo)) = 0, 



(2.2.5) 



i.e. if the divisor 5 is nonspecial: we assume that this is the case throughout. 



^This follows from Riemann-Roch theorem, upon using that degio = 9 + ^ — 1: 

dimH"(C,C'(Lo)) = degLo + 1 - g + dimH^{C,OiLo)) = N + dimH'^{C,0{Lo)). 
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Choose a Ao e C over which there are N distinct points on C, ordered as Pi = 
(Ao, ^i), . . . , Pjv = (Ao, fJ.N), and define the N x N matrbc * as follows: 

(*(Ao))., - V^'(P,)- 

Proposition 2.4. The Lax matrix Ax can be expressed as 

Aa„ -*(Ao)-/x-(*(Ao))-\ (2.2.6] 

where fi = diag(/ii, . . . ,hn)- 

Note that eq. (2.2.6) can be extended by continuity to give Ax in a neighbourhood of Aq, 
and hence globally, since Ax is a polynomial. 

In order to have an explicit expression for Ax, we need an explicit representation for the 
basis Tp'^ of Lq ; we choose one as folios. As the"initial value" line bundle Lq is represented by 
the divisor S, sections of Lq are then meromorphic function of C with poles at S, i.e. 

div(V') > -S. 

Finding a basis for such sections amounts to a gauge choice for Ax. We can fix the gauge by 
imposing the Ax is diagonal for a fixed value of A: we choose the value A = oo. This results in 
the eigenvectors being proportional to the canonical vectors at A = oo. 
Let D°° ^ + ...P^, where the P^ are the N points above A = oo, and let A* ^ 6 - D°° + 
P°° , = - Pf° . It follows from the assumption (2.2.5) that^ 

dimff"(C,0(LAj) = 1- 

This implies that sections of La. also represents sections of Lq = L{5); in particular, for each 
i there is one meromorphic section ip^ such that 

diYiil)"-) > -5 + D\ (2.2.7) 

The V^* may be normalised by^ 

,/,*(P°°) = 1. (2.2.8) 
These last two conditions uniquely determine a basis (V'i) for -ff"(C, 0(Lo))- 

^For a thorough discussion of this point we refer to [BE06], sect. 2.3.2. 

""This is possible since {P^ ) 7^ 0, as the eigenvectors are proportional to the canonical vectors at A = cxd. 
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2.3 Theta functions 

The final step is writing these sections in terms of theta functions; the Riemann theorem 
on the zeros of the theta function is fundamental in this construction and is recalled in Ap- 
pendix A, together with the basic properties of some objects used below, such as the Abel 
map A, the vector of Riemann constants JC, and the theta function 9. 

Write the divisors A'= and 5 as A'^ = X;f=i and S = Ef=f Si. 
Onecanfinde e esuchthat0(e + yl(A,f)) 7^ 0, e{e + AiP°°)) ^ 0, e{e + A^S,)) ^ 0. 
Then a function satisfying (2.2.7) is given by 

First, we notice that the theta functions in the expression (2.3.1) have been assembled so 
that their automorphy factors cancel, i.e. the function is indeed well defined. 

By the Riemann theorem, the zeros of 9 (a{P) 

- J2j=i ^i^j) - ^) are the g points 

5i, . . . , (5g, hence has poles at these points. 



As for the rest, we notice that the second fraction in (2.3.1) can be expressed as a product 
of iV- 1 quotients of the form 9{A{P)-A{P^)+e) /9{A{P)- A{5i) + e), for appropriate k and 
/. Let us analyse the properties of 0{A{P) — A{Q) + e) first. From the Riemann theorem, we 
know it has g zeros: one of them is at Pi = Q, and the others, say P2, ■ ■ ■ ,Pg, are independent 
of it. Indeed, it follows from the Riemann theorem thatA{Pi) + . . . + A{Pg) = A{Q) + e — /C. 
Since Pi = Q, the remaining g -1 zeros are determined by an equation that is independent 
of Q. In particular, this implies that the function 9{A{P) - A{P;^) + e) /0iA{P) - A{di) + e) 
has only one zero at P^ and one pole at Sf, the other 5 — 1 zeros of the theta functions at nu- 
merator and denominator cancel. So overall the product contributes - 1 poles at 5i,. .. ,Sg 
andA^ — 1 zeros at P°°, with j ^ i. 



Hence, a section normalised as in eq. (2.2.8) is given by 

^UP) ^ S^-J- (2.3.2) 

2.3. 1 Time evolution and Baker- Akhiezer functions 

In the previous section we have reconstructed the Lax matrbc ^1^(0) for a specific "initial 
value" line bundle at i = 0. To describe the flow of matricial polynomials A\{t) corre- 
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spending to a linear flow of line bundles Lt we can repeat the above construction for every 
time t to obtain a time dependant basis for H^{C, 0{Lt)), and hence a time dependent Lax 
matrix. Instead of doing this, though, the time dependance can be made more explicit using 
Baker-Akhiezer functions rather than meromorphic functions. 

The sections i^l of Lt may be represented as functions with zeros at A^, poles at 6 and ex- 
ponential singularities on a divisor D. 



Theorem 2.5 {Krichever, 1977). Given a nonspecial divisor 5 of degree g + ?i - 1, there exists a 
unique function ipj (i, P) and local coordinates (P) for which ^ij (Dj) = 0, such that 

1. the function ipj ofP e S is meromorphic outside D and has at most simple poles at 5 (if 
all of them are distinct); 

2. in a neighbourhood of the puncture D I the function takes the form 

{z, P) = e^""'"" (^ji + a% [z) w;f ^ , wi=wi{P), m G N+. 

The integer m is arbitrary, but in applications it is determined by the given flow; in the case 
of monopoles it is 1. 



Introduce the second kind differential fiJo such that^ 



-d(/i7'" + 0(/ij)) ne&vDj- (2.3.3) 
= 0; (2.3.4) 

XI'q"' is holomorphic away from D. (2.3.5) 



t^!™' = 7^A ^n'- (2.3.6) 



Introduce then the vector 

The function ijjj of the above theorem is then given, for our case, by defining 

e{A{P) + tul-^-A{A^)-JC) ,oiH _ 

e{A{P)-Am-K.) ' ^^-^"^^ 

^Note that in the monopole case there we have used a different second kind differential 700 1 with different nor- 
malisation properties (1.5.7). The only difference lies in a different normalization for the exponential exp(/j^ td}j^^ ) 
in (2.3.7). 
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hence the expression for a normalised section V'J of Lt reads 

^liP)^^^^m^. (2.3.8) 

From the considerations after eq. (2.3.1), we conclude that this function is well defined, 
and has the required poles and zeros; moreover, it has exponential singularities at D. 
We point out that the multiplication by h\ does not add any extraneous poles, as its denomi- 
nator cancels with one of the factors at the numerator of eq. (2.3.1). Thus, the poles of ^t(F) 
only come from those of (2.3.1), hence are independent of time. 

Denote by DQ{t) the divisor of the zeros of h\, which hence contains all the time evolution 
information: by the Riemann theorem it satisfies 

A{Do{t)) - ^(A^) = (2.3.9) 

So we recover the fact that the flow is linear onJ{C). 

To conclude, using this theta function representation in eq. (2.2.6), we obtain a solution 
for A\{t) in terms of theta functions, realised as a linear flow on J'(C). 



2.4 General remarks 

This chapter has been a quick overview of algebraic-geometric methods in solving Lax sys- 
tems, to give a general idea of why these methods are so powerful in this context. Firstly, 
algebraic geometric techniques have a wide range of applications, since many known inte- 
grable systems admit Lax pairs: from classical completely integrable Hamiltonian systems, 
to many nonlinear PDEs, to quantum systems as well. 

Moreover, in this context, one manages to make precise sense of the idea that an inte- 
grable system admits a "simple" or closed solution, either "by quadratures", or expressible in 
terms of known functions, or linear in terms, for instance, of angle action variables. We see 
that for Lax systems the flow is indeed linear with respect to the given linear structure on the 
Jacobian of the spectral curve. 

The most important outcome of these algebraic-geometric techniques is probably the 
possibility of obtaining an explicit expression for the dynamical flow of a Lax system in terms 
of Baker-Akhiezer functions. The beauty and value of this result lie in the generality and 
compactness of this expression. In practice, though, the concrete realisation of this formula 
requires a deep understanding of the spectral curve, to obtain an explicit expression for the 



43 



2.4. General REMARKS 



objects needed, i.e. the vector of Riemann constants, the theta function, hence the period 
matrix, which in turn depends on a specific choice of homology basis. This is easier to do for 
the case of elliptic or hyperelliptic curves, but already for trigonal curves it becomes a more 
complicated task: this is, for instance, the case of the monopole curve (1.4.2), which is the 
object of the present thesis. 

It is often the case, however, that in the presence of symmetry, even for more complicated 
spectral curves it is possible to simplify eq. (2.3.8) by reducing the solution. Introducing a 
particularly symmetric homology basis, the period matrix takes a simpler form and it is then 
possible to express the theta functions in terms of lower genus ones, often hyperelliptic or 
elliptic. Again, this is the approach we take for the cyclic charge 3 monopole. 
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CHAPTER 

3 

The spectral curve 



This chapter is dedicated to the study of a (family of) curve(s) corresponding to 
a charge three monopole symmetric under the cyclic group C3, introduced by 
Hitchin, Manton and Murray in [HMM95] (see also section 1.4). Some general 
properties of the curve are examined, in particular its C3 symmetry and the 
corresponding quotient, in view of Fay's theory of unbranched covers [Fay73]. 
The main result is the explicit construction of a particular symmetric basis, 
which satisfies Fay's theorem. This allows to express several objects on this 
genus 4 curve in terms of objects on the genus 2 quotient: in particular, the 
Ercolani-Sinha constraints can be reduced to relations between periods of a 
genus 2 curve, hence becoming more manageable from a computational point 
of view. Throughout, a constant comparison is made with the class of symmet- 
ric monopoles studied byBraden and Enolskii in [BE06], corresponding to the 
limit a ^ of the present case, which serves both as a motivation and a con- 
sistency check. 



3. 1 The Spectral Curve 

Some general properties of the curve introduced in eq. (1.4.2) are presented here; we remark 
incidentally that we should more correctly refer to this as a class of curves, depending on the 
parameters a, /3, 7: for each value of these parameters we have a curve, which corresponds 
to a monopole only if it satisfies the Ercolani-Sinha constraints of Theorem 1.6. 
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3.1.1 Properties 

Let us consider the Riemann surface X defined by tfie following equation: 



+ awz^ + I3z^ + - /3 = 0, 



(3.1.1) 



where^ a, /3, 7 € M. 

The surface X is not hyperelliptic, and has genus 4, hence a basis for the holomorphic dif- 
ferentials consists of 4 elements, for which a standard choice is 



zdz 



Ul 



Us 



3w7^ + az"^ 
z^dz 
+ az^ 



U2 



U4 



Sw"^ + az^ ' 
wdz 
+ az^ 



(3.1.2) 



The above curve can be seen as a cover of the complex sphere with 3 sheets, and 12 ramifi- 
cation points, whose z-coordinates are given by 



B^ 2 - 



L (-97 - 2zV3a^/^ ± A^f 



18^ 



1/3 



^3,4 - 



18^ 



B 



i+4k 



(3.1.3) 

i = 1,...,4, k = 0,l,2, 



where 



A± = 324/3^ - 3(2a3/2 ± SiV^jY' 



p = 63'". 



In Figure 3.1 we give a qualitative sketch of the branchpoints; we remark that, despite hav- 
ing chosen specific values for a, /3, 7 for Figure 3.1, this is the picture of the curve X that we 
keep in mind for the rest of the present work: indeed, there is no lack of generality in doing 
so, as the general properties of X do not change with the parameters (unless a — 0, which is 
a degenerate case, examined in section 3.2); this is due to the fact that the relative positions 
of the z-coordinates of the Bi do not change varying a, {3, 7. The values for the parameters in 
Figure 3.1 are chosen for numerical convenience. 

Moreover, we are also able to find explicitly the monodromy around each branchpoint: 



Bi, B2, Bj, Bs 

B3, B4, Bg, Bio 

B5, Be, Bii, B12 



[1,3] 
[1,2] 
[2,3] 



'The most general form for this surface would be + awz'^ + fiz^ + -yz^ 



= 0, with a, /3, 7 e C. But for 
it to be the spectral curve of a monopole, it has to be invariant under the real involution r (condition HI in section 
1.3); this implies that a, 7 e M; moreover it can be easily seen that (3 can be made real with a change of coordinate 
(a rotation of z of an opportune angle) 
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[1,2] 

n ,3] 

[2,3] 



Figure 3.1: Branchpoints and monodromy for X 

This is represented as well in Figure 3.1: as indicated in the figure, a continuous line means 
that around branchpoint Bi there is monodromy [1,2], and so on. More details about mon- 
odromy are given in section 3.3, here we only remark that it is calculated numerically, but it 
doesn't change varying the parameters (as long asa^O). 



The curve admits the following automorphisms 

inversion : 
real involution : 
cyclic involution : 



1 w 

^ 2 



(/) : (z,w) - 
T : (z, w) — 
a : (z, w) — > {pz, pw) 



z z- 
1 w 

z 



(3.1.4) 
(3.1.5) 
(3.1.6) 



Comparing the above with the symmetries of the curve (1.6.1) studied in [BE06], we note that 
the inversion and the real involution appear in that case too, the former being an inversion 
with respect to the origin, and the latter corresponding to the reality condition of constraint 
HI in Hitchin formulation (see section 1.3). 

The cyclic symmetry is the symmetry by the action of C3: this is again a symmetry of the 
curve (1.6.1), and in the notation of section 1.6.1 it reads a = TZo s. This last automorphism 
indeed characterises this family of monopoles, as explained in section 1.4.2; in fact, this is 
the main focus of the present investigation: we are hence going to exploit this symmetry 
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further in what follows. 

3.1.2 The quotient with respect to C3 

First, we remark that a sends the branchpoints one into the other as follows (cf. eq. (3.1.3)) 



Bi 


a 


■B5 


17 


■ ^9, 


B2 


a 


■Be 


<7 


■ ^10 


B3 


a 


■B7 


<7 


■^11 


B4 


a 


■Bs 


<7 


■ B12 



The quotient under the cyclic involution is described by the map^ 



TT : (z, «;) ^ (x, 2/) = I -, + ^) ) . (3.1.7) 



Z Z' 



So the quotient curve, which we denote X, is 

y2 (a;^ + ax + 7)^ + 4^^ (3.1.8) 

It is a genus 2, hence hyperelliptic, Riemann surface, in standard form. Seen as a 2-sheeted 
cover of the Riemann sphere, it has six branchpoints whose x-coordinates are 

6 J3 6^3 





2p^a 

2 




2a 




2" 5 

51 



- -p d_ 2- , Bq - -p 2- , 

6 ^3 6 J3 



(3.1.9) 



where 

(5± = -108 7 - 216/3 i + 12 ^^12^3 + 81 {-/±2(3if. 
As (5- = these branchpoints can be split in complex conjugate pairs 

Be^Bl, B5^B^, Bi^B^. (3.1.10) 

We remark in particular that the quotient map tt : X ^ X is an unbranched cover of X 
by X. This can be seen using the Riemann-Hurwitz formula, i.e. 

2 • (5 - 1) = 2 • deg(7r) • (g- 1) +^mult(7r)p =^ ^ niult(7r)p = 



2 Note that in fact it is the composition of the quotient map with a map that puts the quotient curve in hyperel- 
liptic standard form. 



48 



3.1. The Spectral Curve 




2 




4 





Figure 3.2: Branchpoints and monodromy for X 



This fact is very relevant, in that it allows us to use the theory of unbranched covers (see sec- 
tion 3.4) with interesting results. 

Note that the branchpoints B.^ are not the images of the branchpoints of X under tt: this 
is due to the form (3.1.7) of the map tt, which mixes the coordinates z and w. 
Figure 3.2 shows again a qualitative sketch of the branchpoints for the curve X, and their 
monodromy, with the same choice of parameters of Figure 3.1, and the same caveats as in 
the previous case; more details about monodromy are given in section 3.3. 

A standard basis for the holomorphic differentials on X is given by 

da; xdx rQ 1 n A 

Ui = , U2 = . (3.1.11) 

y y 

We notice that the differentials U2 and U4 on X are invariant under a and hence descend to 
differentials on X; or better, in other terms: 



This observation allows us to considerably simplify some integrals, and hence the period 
matrix: this is examined in section 3.9. 



TT* Ui = -3U2, 



TT* U2 — — 3U4. 



(3.1.12) 
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3 n 




^ Branchpoints X ^ Branchpoints B . varying a 



Figure 3.3: Branchpoints for a — > 

3.2 The case a = 

This section is dedicated to the analysis of the spectral curve in the case q = 0, namely the 
case studied by Braden and Enolskii in [BE06], and reviewed in section 1.6. This is indeed 
a genuinely different case, in that we have truly different curves, with different properties, 
e.g. more symmetry. We are particularly interested in making contact with the case a = 0, as 
the present work aims to generalise the results in [BE06] , and therefore, it is useful to have 
that as both a comparison and an inspiration in what follows. 

In the case a = we recover the curve of eq. (1.6.1), which we report here again for 
convenience 

+ + - /? = 0, (3.2. 1) 

together with some of its basic properties. The corresponding Riemann surface, which we 
denote by Xq, also has genus 4, but only 6 branchpoints, Ai, i — \ . . given in (1.6.2). In- 
deed, letting a — > in eq. (3.1.3), we see that the branchpoints collide pairwise to give the 
\i (see Figure 3.3). 

In Figure 3.4 we also give the monodromy, which is the same for every branchpoint, 
namely [1,2,3]; this can be seen from an explicit calculation, but also by appropriately taking 
the limit of the monodromies of X. This is explained in detail in section 3.3. 
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The basis of the holomorphic differentials used in the following is 

dz , Q zdz ,Q z'^dz „Q w dz 

= U2 — U3 = U4 — 

w w w w 

Note that the ordering is different from that used in [BE06], as this basis corresponds, up to 
a factor of 1/3 to the one for X given in eq. (3.1.2) 

lim u, = -u?. 

a^o ' 3 ' 

Taking the limit for a — > in the quotient curve X, we obtain the curve Xq. We notice that 
the limiting operation does not change the curve dramatically as it changes the curve X: in- 
deed, no degeneration appears, unlike the X case, and the branchpoints just move slightly; 
the one remarkable effect is that this curve still has the extra symmetry observed above. 



Notation We recall that all the objects with a hat are objects on the genus 4 
curve X (3.1.1), i.e. the spectral curve, while objects without a hat belong to 
the quotient curve X, obtained quotientingX of (3.1.8) with respect to the C3 
action. Also, objects with a hat and a superscript are intended to be on the 
limit for a — > of the curve X (3.2.1). Summarizing this in a diagram: 



7r 

>■ 

X 
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3.3 Remarks 

In this section we examine in more detail some technical aspects of the previous work. 
Throughout this work, we view a Riemann surface S, given by an equation of the form 
V{z,w) = 0, as a branched n-sheeted cover of the Riemann sphere, with covering map 
TT : S . In practice, this means that we treat it as n copies of the complex plane, each cor- 
responding to a sheet, that we glue appropriately along the branchcuts to obtain the surface. 
First we examine in detail how to explicitly calculate monodromy and how to order sheets 
in practice given a Riemann surface, emphasising the fact that these are two distinct oper- 
ations. We then explain how they actually tie together, giving information on how to glue 
sheets to a branchcut. All this is needed to understand the picture of the Riemann surfaces 
X, Xo and X of the last section, and their homology bases, which are introduced later in this 
chapter: therefore we conclude this section applying these considerations to our specific 
case. 

Monodromy 

Let us remark a very important point (which is often overlooked in the literature) : the mon- 
odromy around a branchpoint depends on the basepoint, the point where the loop is based. 
Indeed, we recall that the monodromy representation of a holomorphic map tt : X ^ Y is 
defined as the representation of the fundamental group based at a point of the base, tti {Y, P) 
in terms of permutations of n elements, n being the number of preimages of P under tt (cf. 
e.g. [Mir91],p. 84). 

So, in our explicit calculation of the monodromy around a branchpoint B^, we first choose 
a reference point on the base P^, which is the same for every branchpoint, and we or- 
der its preimages on X, Pi = {z^,wf),. . ., P„ = (z^, w^); then we consider, on the base 
space P^, a loop based at zq around Bi. We lift this closed path to the Riemann surface 
(using the Acontinuation command in Maple); if the lifted path is an open path, its ini- 
tial point is going to be Pi and its final point Pj: then we say that the monodromy around 
the branchpoint has cycle structure . . .]. In other words, we start from an ordering of 
the w-values above z^, i.e. (wf , , . . . , w^), and analytically continuing these values along 
path enrcircling the ramification point, we return to the same set of w-values, but per- 
muted (^f^(i), w'^(2)i ■ • ■ : ^^(n))- ^^^^ defines a permutation o-^, the monodromy data for 
the branchpoint S,; the permutations ai generate the monodromy group based at z^. A dif- 
ferent choice of ordering the ui -values at z^ leads to a conjugate set of permutations and 
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Figure 3.5: Cuts and paths for monodromy calculations 



isomorphic monodromy group. 

In this picture, the concept of branchcut is introduced, namely, using the classical ter- 
minology, we "cut" the complex sphere by eliminating, for every branchpoint Bi, a curve 
(usually a straight line) connecting Bi to the basepoint Pq, hence obtaining the sheet, which 
is a plane. In practice, to calculate the monodromy the closed path we use is the one as in 
Figure 3.5: qualitatively, the path follows the cut on one side, then encircles the branchpoint 
along a small circle of radius e, then follows the cut on the other side. Note that conven- 
tionally this path is oriented such that it goes around the branchpoint anticlockwise. Notice 
also that the monodromy depends on how we choose the cut: for instance, if in Figure 3.5 
we consider the dashed cut and follow the path around it as just described, the monodromy 
around Bi is different, and in fact conjugate to that obtained following the other cut. Thus, 
when we talk about monodromy around a branchpoint we always mean monodromy calcu- 
lated with respect to a specific choice both of basepoint, and of cut; in particular, in section 
3.1, the monodromy is always understood to be calculated along the path explained above, 
using the point as infinity as basepoint, and the cuts as in Figures 3.1, 3.4. Details about this 
are given at the end of this section. 

Sheet ordering 

In the picture above, we wish to think of the n cut planes as sheets. Firstly, we need a specific 
way of identifying, or labelling, sheets. The notion of sheet ordering is distinct from our or- 
dering at zq. In practice we achieve this as follows: we fix another (regular) point on P\ the 
base of the covering, say z^. Corresponding to it we have n points on the Riemann surface 
Si = (z-^, wf ), where wf are solutions of V{z^,w) = 0; now, we order the wf, and we say 
that wf labels sheet 1, etc. . The points Si :— {z^ , wf) are then our reference points for each 
sheet: to understand on which sheet a given point Q = (z*^, w'^) is, we have to compare its 
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Figure 3.6: Sheet labeling 



w-value with that of one of the reference points Si. 

In order to do so, choose a path 7 on starting from z'^ and ending on z^, with the only 
constraint that it does not have to cross any cuts. Then lift it to a path 76!] (for instance, us- 
ing the Acontinuation command in Maple) with the condition that 7(0) — Q: the endpoint 
of 7 on S is one of the Si, which tells us then on which sheet Q is. This is pictured in Figure 3.6. 



In this work we have always chosen z = as reference point for sheet ordering (or sheet 
base, in short), and have ordered the sheets as follows: 



sheet base z = 



sheet 1 w = a e E 
sheet 2 w = a 
sheet 3 w = p a. 



(3.3.1) 



This ordering in the limit a — > yields the sheets ordering in [BE06], thus making the com- 
parison with the Braden-Enolskii case much easier. 



Monodromy and gluing sheets 

Before explaining in detail how to calculate monodromy for the curve X (where we choose 
the basepoint to be at infinity), we clarify one further very important point concerning the 
very meaning of monodromy. 

The notion of sheet ordering at is distinct from the ordering at z^. Indeed, depending 
on the path chosen to link z^ with we could reach a different wf, so the 7«-values at z^ do 
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Figure 3.7: Example I 

not distinguish sheets alone. In practice which wf one reaches depends on the "sector" be- 
tween cuts where the path approaches z^. Sheets and monodromy then appear as distinct 
constructs-^. 

In particular one cannot say just from the monodromy information at a branchpoint which 
sheets come together, or vice versa. In fact, the monodromy data around a branchpoint, to- 
gether with the sheet ordering, gives information about how to attach the n copies of C to a 
given cut, rather than about which sheets are connected at a certain cut. Thus, if a branch- 
point has monodromy [i, j, . . .] , it does not necessarily mean that crossing the corresponding 
cut one goes from sheet i to sheet j. We make this more clear from an example, as it is a 
fundamental point in understanding the homology basis of section 3.7. 

Example 3.3.1. Consider the curve of equation 

+ z^w + + 1 = 0, 

with 6 branchpoints 

hl^d-K hl^pdrj, h^^pd+K hl^p'd-J, bi=Ai, hl^dlK 
with d± = 1 ± 2 and p — exp(2i7r/3); they are represented in Figure 3.7 as grey spheres. 

Take the point — l/2exp(i7r/3) as basepoint for the sheet ordering, which lies in the 
^ We thank Tim Northover for many useful discussions concerning this point. 
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b2 



Figure 3.8: Example II 



"sector" between bi and 1)2, and choose an ordering for the sheets (wf , wf > "'f ) (the black 
dots in Figure 3.7). This means, looking at Figure 3.7, that sheet 1 corresponds to the colour 
red, sheets 2 to blue and sheet 3 to green. 

Take the origin to be z^, the reference point for monodromy, and impose an ordering 
of the w[ out of that at , by analytic continuation along a straight line from z^ to z^: 
{w^, w[) = (-1, -p, -p^). In the figure this is represented by the height. Then the mon- 
odromy around each branchpoint is computed to be 



bi, b4 ^ [-1, -p^] = [1, 3], b2, bs ^ [-P, -p^] = [2, 3], b3, be ^ [-1, -p] = [1, 2], 



and in the figure is represented as a change in height. 

This pictorial representation clarifies the fundamental differences between sheets and 
monodromy, together with their link: the sheet ordering at imposes an ordering of the 
preimages wf of z^; the monodromy at a branchpoint is a permutation of this w[ ordering, 
after going around a branchpoint (represented in Figure 3.7 by a change in height); this tells 
us how to attach a sheet at each cut. 

For instance, starting in the sector between bi and 1)2, the preimages of z^ are ordered 
(wf , wf, ); going around bi, the monodromy is cri = [1,3], hence the order is now 
("'ffi(i)' ^,^1(2)' ""^^1(3)) = (^f' ^D- t'^i^ means that sheets 1 and 3 come together at this 
cut. Going from there around bg, the monodromy is o-g = [1:2], hence the order is now 
(^ffeCaj'^^eCaj'^^lli)) ^ (^fj^fj^f)- comparing with the ordering we had in the previous 
"sector", this means that sheets 2 and 3 come together at this cut. This is sketched in Figure 
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Figure 3.9: Branchpoints and cuts under inversion 



3.8. 

From this we see that starting at the point (z^, w[) on sheet 1, the red one, after crossing the 
cut at bi we end up on the green sheet, namely sheet 3; after crossing one more cut at we 
end up on the blue sheet, namely sheet 2; and so on. 

One can continue in this fashion for every branchpoint: the case of the cuts at bi and bg is 
already explicative of the general situation. 

Monodromy calculation for X and Xq 

We can now clarify the technical details about how we calculate monodromy in section 3.1, 
where we consider, for both curves, the point at infinity as the basepoint for monodromy, 
and the cuts are straight lines from the branchpoints to infinity. The way to deal with this, 
numerically, is to use the inversion (3.1.4) to bring the point at infinity to the origin, hence 
making all the paths finite. Note also that the branchpoints get mapped one into the other 
as follows: 

HBi) = Be, =55, 0(^3) =B4, 

(and the rest are obtained recalling that </)^ =1): so, for instance, the monodromy around 
Bi with cuts going to infinity is obtained from the monodromy around Br with cuts going to 
zero. This is represented in Figure 3.9. 

There is just a small subtlety, namely that the sheet ordering of (3.3.1), i.e. using the ori- 
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Figure 3.10: Cuts for a -j> 



gin z = as the sheet base^, is not preserved by (p: hence, to obtain the sheet ordering of eq. 
(3.3.1) in the case of outwards cuts to infinity, we have to choose a different ordering in the 
case with cuts to zero, induced from the previous one via 0. 

This last remark allows us to clarify a final point, namely what happens in the limit when 
a — )■ 0. In this limit the branchpoints in X collide pairwise, and so do the corresponding 
cuts: from this one can deduce the monodromy around the resulting cut out of that of the 
starting two cuts. Let us focus for simplicity on the colliding branchpoints B12 and Bi : as 
said earlier, in practice we calculate the monodromy after inversion, so the cuts we consider 
are those of Figure 3.10, with basepoint . Lifting a closed loop based at , if it starts, say, 
at Pi, going around Bi it ends at P3, then going around B12 it ends at Pi, now, in the limit, 
the cuts collide, and hence only the initial and final part of the path remain, meaning that, in 
the limit a — > 0, if the lifted path starts at Pi, going around Ai, it ends at P2 (as in the upper 
p art of Figure 3.10). 

This can be seen in an equivalent way, considering paths crossing the cuts, as shown in the 
bottom part of Figure 3.10. 

■* Since the origin is already used as monodromy base, we cannot take it as sheet base as well; this is solved by 
slightly moving the sheet base away from the origin (we choose the point ^ = 0.1 + 0.1i), this also has the advantage 
that after inversion it corresponds to a finite point, which makes it easier to compare the respective ordering of the 
«)-values by analytic continuation to the Si. 
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Doing this again for the other cases one obtains that the monodromy around Ai is [1, 2, 3] 
along the cut going inwards; as explained earlier, the ordering of the sheets has been cho- 
sen to match before and after inversion, so we can conclude that the monodromy around Ai 
along the cut going outwards is indeed [1, 2, 3]. 

One proceeds in the same fashion the other branchpoints (and corresponding cuts), to ob- 
tain the monodromy in Figure 3.4. 

3.4 Interlude: unbranched quotients 

In this section we describe some results in the theory of unbranched covers, which plays an 
important role in the study of the curve X, and simplifies the problem significantly. 
The general setting is that of a Riemann surface with a cyclic group acting on it, such that the 
quotient induced by this action is unbranched. 

In the first subsection we present a theorem of Fay [Fay73] about the existence, under these 
circumstances, of a basis for the first homology group with certain symmetry properties. 
This proves to be crucial in expressing several objects on our curve X in a much simpler way 
[Fay73, Acc7 1] ; moreover, we remark that the theorem is non constructive, and there are few 
examples in the literature of such a basis: here we provide a new one. 

3.4. 1 Fay's symmetric homology basis 

Consider a Riemann surface S of genus g which admits an (effective holomorphic) action of 
a finite cyclic group G of order p 

$ : E — >t. 

This action induces a quotient 

TT : S — > S, 

where E is another Riemann surface of genus g. We consider the case in which the covering 
TT has no branchpoints. 

It is then possible to find a canonical basis (see Appendix A) of Hi (E, Z), whose 2g cycles 
we denote by^ 

So, oi, . . . , dg_i; bo, bi, . . . , bg_i; (3.4.1) 

^Throughout the following discussion quantities with a hat are understood to be on the "original" curve, quanti- 
ties without a hat are on the quotient curve. 
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and a canonical basis of Hi (S, Z), whose 2g cycles we denote by 



with the following properties with respect to the lift of the covering tt* 



(3.4.2) 



7''*(afc + „(g_l)) — Ofe, 

'^*ibk+n{g-l)) = bk, 



7r*(oo) = ao, 
7r*(bo) ^P, bo; 

n = 0, 1; k = 1, 



(3.4.3) 



In particular, this implies that the cycles on E have the following behaviour under 



$"(afe) =afe+„(g_i), 

$"(&fc) =bfe+n(g-l)- 



(3.4.4) 



<i':^(bo) =bo. 



Here, and throughout in the following, the equality signs mean equals in Hi{T., Z), so when 
representing the cycles as closed paths on the surfaces we have to replace "equals" with "ho- 
mologous". 

The above properties can be stated at the level of the holomorphic differentials, i.e. for 
iJi(S,M). If c2>o,--,'^g-i is a basis for i/i(E,M) dual to that in eq. (3.4.1), then eqs. (3.4.4) 
become 

<i>*"(u>o) = 'I'o, ^*''iu:k) = ^fe+„(3-i). (3.4.5) 



3.4.2 Simplification in the period matrix, theta functions and vector of 
Riemann Constants 

Using a basis satisfying the symmetry properties above, one manages to obtain a particularly 
symmetric expressions for the Riemann period matrbc, and a consequent simplification in 
the associated theta function. 

Fay-Accola Theorem 3.1. With respectto the ordered canonical homology bases ho, . . . , ag_i, 

bo, bi, . . . , b-g-i and ao, • • • Og_i, bo, bi, • • • bg_i satisfying (3.4.4), the a-normalised Riemann 
period matrices of S and S take the respective forms 



I 



Too tq To 



To 



fj" M M' ■■■ M(p-i) 



V ^0 



fj M' M" 



M 



1 



-Too To 



N 



(3.4.6) 
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where 



'4>"bo Jit>"bj 
Nij = Tij + Tij+g_l + . . . + T'i,j(p-1)(5-1)' *7 J = I7 • ■ -3 - 1 

This particularly symmetric form for the period matrices allows to simplify both the theta 
function and the vector of Riemann constants for S, that can be expressed in terms their 
correspondents on E. 

Recall that the map tt : E ^ S can be lifted to a map tt* ; J'(S) -> Ji^)> and hence to a map 
from to C^, which we still denote tt*: with the above choices (3.4.1), (3.4.2) of homology 
bases, we have, for z e 

TT*Z = 7r*(zo,Zl, . . .,%-!) = (pZo.Zi, . . . . . . ,Zi, . . .,Zg_i), 

or in characteristic notation 



TT 



ao ai ... ag_i \ ^ / ao ai ... ttg-i ... ai ... ag_i 
/3o /3i ■ . ■ /Sg-i / \pl3o Pi ... ... /3i ... 



We can then state the following theorem, proven by Accola^ [Acc71] and Fay [Fay73] in- 
dependently. 



Fay-Accola Theorem 3.2. For z e C'-' one has 

0[e]{7r*z;T) 



where cq{t) is a non-zero constant independent of z ancf 



= co(f). (3.4.7) 



„ ,0, ...,0 , e = 7r*e. (3.4.8) 
2p J 

Moreover, the following relation between the two vectors of Riemann constants holds 

p-i 

tt*Kq = /Cq + (5 - 1) ^ -4($"(g)) + e. (3.4.9) 

n=0 



^Note that in Accola's proof the covering need not be induced by the action of a cyclic group, but here we give the 
result in this particular case as it is the case under consideration. Also, in this more general setting, the form of e is 
not known explicitly. 

^Note that this differs from proposition 4.1 in [Fay73] by a factor of 2 in the denominator in the first entry: this is 
due to a different choice of conventions in the basis we use for the Jacobian, which differs by a factor 2ix from what 
Fay considers. 
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3.5 A symmetric homology basis: preliminaries 
3.5. 1 Fay's symmetric basis for the curve X 

As the curve X satisfies tlie Jiypotiieses above, we can now restate tlie properties that Fay's 
basis has to satisfy in this particuiar case. Using the same notation as above (cf. eqs. (3.4.1), 
(3.4.2)), eqs. (3.4.3) become 



7r*(bj) = bi, 
and eqs. (3.4.4) become 



7r*(ao)==ao, (3.5.1) 
7r*(bo) = 3bo, i = 1,2,3, 



(3.5.2) 

<(bo) = &o, n= 1,2,3. 



It turns out that in this case it is most effective to find a basis on X satisfying (3.5.2), and then 
project it to obtain a basis on X satisfying (3.4.3). In 77^(E,Z), eqs. (3.5.2) mean that one 
cycle stays invariant under the group action, and three pairs are mapped one into the other 
by this action. 



3.5.2 do is invariant under a 

We remaric that Fay's conditions in the case of a finite cyclic group of order 3 require oq to be 
invariant as well. 

Indeed, from equations (3.4.3) it follows that the matrrK representing the involution cr on a 
basis with such symmetry properties is of the following form 



/ 



Ql 0-2 0,3 0-4 15 16 fly 



Only the first row is arbitrary, as there are no conditions on how ao behaves under <t* . But 
imposing that that ct is a symplectic transformation and that al is the identity is enough to 
fix all the ai to vanish except for ai = 1. Hence the cycle do is indeed invariant. 
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Therefore, we can rewrite Fay's conditions (3.5.2) as follows 



11 J 



o-"(ao) = Qo, 



(3.5.3) 



- bo. 



n = 1,2,3. 



These are the conditions we focus on in the following. 

As remarked above, here the equalities are to be intended in Hi{X); so, when we repre- 
sent the cycles by closed paths in X, the equal signs should actually mean "homologous". 
Indeed, it turns out that even at the level of paths "true" equalities hold: the only exception 
is do, whose images under a'^ axe just homologous. 

3.6 The relation to the case a = 

It turns out that the above requirements (3.5.3) do not determine a homology basis uniquely: 
indeed, we found many. An obvious ambiguity can be seen already at this stage, for in- 
stance, given a basis satisfying (3.5.3), cyclically interchanging the cycles Oj, hi,i = 1,2,3 
yields a basis with the same symmetry; there is in fact even more freedom in the choice of 
the dj, bi, i = 1, 2, 3. Therefore we can impose some extra conditions on this basis, in partic- 
ular that in the limit a ^ some of the cycles are mapped to those of the homology basis 
in [BE06] . This would be a desirable feature as it allows us to relate the present work with 
[BE06] , whose results we aim to generalise. 

We will be more precise on this in section 3.7, and we refer to section 1.6.1 for more details 
on the basis of [BE06] : here we only give again this basis in Figure 3.11. We remark again that 
here the monodromy (calculated again with respect to infinity) is the same for every branch- 
points, namely [1, 2, 3]. The notation for the cycles is explained in the figure, and is used 
throughout this work. We also recall that this basis satisfies different symmetry properties 
than those we are interested in, see eqs. (1.6.9). 

3.7 A symmetric homology basis 

In this section we find explicitly a basis for Hi {X) satisfying Fay's symmetry conditions, as in 
eqs. (3.5.3). 

We remark again that, as Fay's result is non-constructive, providing an example is not a trivial 
task. Here we do not provide an algorithm to implement this; we think that indeed this would 
be possible, in principle. In fact, there are already algorithmic ways to construct a homology 
basis: for instance, Tretkoff and Tretkoff algorithm (see [TT84] and Appendix B), which is al- 
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Figure 3.11: Symmetric homology basis of [BE06] 
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gebraic in nature. Once the Tretkoff and Tretkoff 's basis is found, one could impose on it the 
required symmetry conditions, and get a symmetric basis in terms of the original Tretkoff 
and Tretkoff's basis. The problem with this is that usually one does not get a "simple" result: 
for instance, the Tretkoff and Tretkoff's basis is already quite complex (in the sense that the 
cycles are very involved, cf. Appendbc B), and a linear combination of its elements is even 
more so. This defeats the whole purpose of having a basis with a certain symmetry, which is 
indeed to simplify the problem. In particular, in the calculation of the periods, having "sim- 
pler" cycles leads to a much simpler and compact result, see section 3.9. 

Hence, the way we follow is not algorithmic; nevertheless, the methods used provide 
some insight about how to tackle this sort of problems, and the same logic can be used in 
similar cases. 

Moreover, we emphasise that many results described below were greatly simplified by Tim 
Northover's package for Maple described in Appendbc C.2. 

3.7.1 The cycles dj, bi 

First, we focus on the search of the three pairs dj, b^, with the condition that they are mapped 
to one another by a. 

As mentioned earlier, we want to relate out constructions as much as possible with those 
a = as in [BE06] , hence initially we work in this limit, and then we try to extend the cycles 
to the case a ^ 0. We have remarked earlier that the cycles a^, b^, z = 1 ... 3 in Figure 3.11 
do not have Fay's symmetry properties. Nevertheless, we can choose one pair of them, say 
ai, hi, to be a", b°, and then obtain the others applying a: 

a; = ai, b° = bi, (3.7.1) 

a^+i=a^ai, b°+i=a,^bi, fc = l,2. (3.7.2) 

They are shown in the upper part of Figure 3.12. 

Note that the cycles thus chosen still satisfy the second equality of (1.6.9), namely 

s,{bl)=nloaM)^cil 1,2,3. 

To obtain the corresponding cycles on X, we exploit some properties of the limit a ^ 0. 
We recall that in this limit the branchpoints in X and the corresponding cuts collide pair- 
wise; this means that if a cycle crosses both cuts, its arc between the cuts becomes degener- 
ate. This is illustrated in the example in the bottom part of Figure 3.10, where we examine 
the cuts corresponding to B12 and Bi. In that picture consider, for instance, the path on X 
starting on sheet 2, crossing the cut and hence continuing on sheet 3: we see immediately 
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that there is a unique path on X that gives the previous path; the same happens for the other 
possible paths crossing this cut. 

From this example it is then clear how to carry out the inverse process, i.e. starting from a 
cycle in the case a = how to reconstruct its "missing" arc between the collided cuts, unam- 
biguously. Hence, we obtain the cycles a^, bi, as in Figure 3.12 with the required properties as 
ineqs. (3.5.3): 

a:?(bo) = bo, n-1,2,3., 
and also 

lim di = ai, lim bi = bi. (3.7.3) 



3.7.2 The cycles do, bo 

As for the remaining two cycles, we focus first on bo. 

It can be found using similar considerations: bo has to be invariant under a and is a three- 
fold covering under tt of a cycle on X, and needs to have intersection zero with all the other 
cycles. Implementing first the symmetry requirements, we split the z and w planes in three 
sectors, and as above we look for a cycle which has a component of similar shape in each 
sector. From these considerations, we choose a cycle bo as in the lower part of Figure 3.13. 
One can check that its intersection with the a^, bj is zero, either by inspection, or using 
Northover's extcurves (see Appendix C). 

As for oo, we focus first on the condition of its intersection with bo to be -1, and used 
the invariance under a (as shown in section 3.5.2) as a consistency check. Since we want the 
maximal possible contact with [BE06] in the case a = 0, one can check whether one of the 
cycles ao, bo has the correct intersections with the other cycles in the limit: it turns out that 
ao does. As explained in the previous section, we can lift this cycle to the case a nonzero, to 
obtain the cycle in Figure 3.13. 

One can then check, for instance with Northover's code, that the basis so obtained is indeed 
canonical, and moreover that it satisfies the required symmetry properties. In particular, we 
notice that while the cycle bo is manifestly invariant (just looking at the figure), the cycle do 
is invariant up to homology. 
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3.7.3 Summary and arc expansion 

We have found a canonical basis for Hi{X) with the symmetry properties as in (3.5.2), 
namely 

(3.7.4) 

'T:(bo) = bo; n = l,2,3. 
Moreover, its limit for a — > satisfies 

Ti-lo cr* ( lim bi ) = lini di, i = 1, 2, 3, 

and also 

d^ = lim di — ai, = lim bi = bi, dp = lim do = a4. 

a— a— >0 Q— ^0 

These cycles can be expanded in terms of "basic arcs" as follows. 

Denote by {h j) the arc going from branchpoint Bi to branchpoint on sheet k: 

lk(i,j) = aiCk{B,,Bj), i^j = l,...,12 on sheet fc. 

Then we have the following 

di =7i(l, 2) +72(2,1), bi =7i(3,l)+72(l,12)+73(12,3), 

d2 =72(5, 6) +73(6, 5), b2 =72 (7, 5) +73 (5, 4) +71 (4, 7), 

£3 =73(9, 10) +71(10, 9), b3 = 73(11, 9) + 71(9, 8) +72(8, 11), (3.7.5) 

do = 71 (3, 10) + 73(10, 9) + 71(9, 8) + 72(8, 12) + 73(12, 3), 

bo = 71 (2, 8) + 72(8, 11) + 73(11, 4) + 71(4, 7) + 72(7, 12) + 73(12, 2). 

and, as a 0, 

7i(l, 2) + 72(2, 1), b? = 71(2, 1) + 73(1, 2), 
72 (3, 4) + 73 (4, 3) , b^ = 72 (4, 3) + 71 (3,4), 

73(5,6) +71(6,5), b§ =73(6,5) +72(5,6), (3.7.6) 
73(1, 2) +71(2, 6) +73(6, 5) +72(5,1), 
73(1, 2) + 71(2, 5) + 72(5, 6) + 73(6, 3) + 71 (3, 4) + 72(4, 1). 

We point out explicitly that both these basis are indeed canonical, i.e. that the intersection 
matrix for both is** — J, where ^ = ( 04 ) • ^^^^ checked by manually, computing 

^We have chosen tlie minus sign to agree with the conventions of [BE06]. 
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b" 
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the intersections of each pair of cycles in Figures 3.12, 3.13; here it has also been verified us- 
ing Northover's code ext curves, which provides a function to find the intersection matrix of 
a given homology basis (see C.2). 

For completeness we also report here, for the case a = 0, the matrix of change of basis 
between the cyclic basis, described above (cf. Figure 3.13,3.12), and the symmetric basis of 
[BE06] (cf. Figure 3.11): 



M 



1 
1 






-1 









-10 



10 



1 
-1 -1 
1 



(3.7.7) 



This matrix has been obtained numerically using Northover's code extcurves for Maple (see 
Appendix C and the Maple file change_basis_alphaO .mw). 



3.8 The projection of the symmetric basis on the quotient 
curve 

Fay's basis projects to a basis for Hi (X) as follows (see [Fay73]) 

7r*(ai) = ai, tt^Oo) = ao (3.8.1) 

^4b,) = bi, 7r,(bo)-3fao; i = 1,2,3. (3.8.2) 

Thus to obtain such a basis one just projects the cycles of Figures 3.12,3.13: the fact that the 
branchpoints of X do not get mapped by n to branchpoints of X (cf. section 3.1.2) makes 
the projection less straightforward. This has been implemented in Maple with Northover's 
extcurves (see Appendix C and Maple file projected_basis .mw), and the result is given in 
Figure 3.14. 



70 



3.8. The projection of the symmetric basis on the quotient curve 





Figure 3.14: Cyclic homology basis for the quotient curve X 
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With the same notation as in the previous section, the arc expansion for these cycles 
reads: 

01 =7i(2,6) + 72(6,2), 
bi =71 (6, 4) +72(4,6), 

(3.8.3) 

ao = 71 (3, 4) + 72 (4, 6) + 71 (6, 1) + 72 ( 1 , 6) + 71 (6, 4) + 72 (4, 3) , 
bo = 71 (3, 4) + 72 (4, 5) + 71 (5, 4) + 72 (4, 3) . 

We therefore have a homology basis for the hyperelliptic curve X which is different from the 
usual one for an hyperelliptic curve in standard form. Since it is the projection of the basis 
for the homology of the genus 4 curve X, it turns out that many objects on X which depend 
on the basis used [e.g. the period matrbc), can be "reduced", namely expressed in terms of 
objects on X. As the curve X is hyperelliptic, and in standard form, calculations are much 
easier in this case, thus the problem simplifies greatly. 



3.9 Symmetries and the period matrix 

We are now in a position to use the basis found in section 3.7 to obtain a very simple form 
for the period matrbc of X. Consider all the symmetries of the monopole curve, the real 
involutionr (3.1.5), the inversion^ (3.1.4), and the cyclic symmetry a- (3.1.6). These maps lift 
as maps from Hi {X, Z) to Hi {X, Z), and their matrix form in the basis of Figures 3.12, 3.13 
is 
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3.9. Symmetries and the period matrix 



The first two are found using Northover's ext curves for Maple , while the last follows from 
the very defining properties of the symmetric basis (cf. section 3.5.2) and has indeed been 
checked with Northover's code. 



These symmetries hugely simplify the period matriK as follows. 
Firstly, we exploit the cyclic the symmetry a; it acts on the differentials as follows 



f7*Ui = p^Ui, 



(7 U3 = p U3, 



Cr*U2 = U2, 



CT U4 = U4. 



(3.9.2) 
(3.9.3) 



Hence, the matrices of periods take the following form 



A 



^ yo ^ 

zi xi wi yi 
p^zi xi pwi yi 



B 



y pzi xi p^wi j/i J 

For instance, we get 



^ Xo Yo\ 
Zi Xi Wi Yi 

\ pZi Xi p'Wi Fi j 



(3.9.4) 



Zq — f Ui = (j) Ui = il) (7 Ui = 



Zq = 0, 



Z2= <P Ui 



Ui = cp a'Ui = (p p^Ui = p^ zi, 



and similarly with the others, to obtain (3.9.4). 



We now use the other symmetry ip, acting on the holomorphic differentials as follows 



(fi Ui^ U3, ip U3= Ui. 



(3.9.5) 



With considerations similar to those reported above for a, we have the following identities 



W2 = -Zi, 



Wi 

Wi 



-Z2, 
-Z2, 



W3 
W3 



'Z3. 



(3.9.6) 
(3.9.7) 



Hence the matrices of periods simplify to 



A 







^2 . 



Xi 



yo 

„2 



p Zl yi 

p-zi xi -Zl yi 
y pzi xi -pzi yi j 



B 







Xo 

\ pZx Xi 





-p^Zi 
-Zl 



Yo\ 

Yi 
Yi 



-pZi Yi j 



(3.9.8) 
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3.9. Symmetries and the period matrix 



As for the real involution, one gets: 

T*ui = U3, T*U3 = ui, (3.9.9) 

T*U2 = — U2, r*U4 = — U4. (3.9.10) 

These relations, together with the above expression (3.9. 1) for M^. , allow us to simplify further 
the structure of the period matrbc. We have then that 

Zi = -{wi + zi) = pzi — zi, Wi — —{zi + wi) = p'^zi — zi, (3.9.11) 

3X1 = 2x0 + xo, 3Yi = 2yo + yo, (3.9.12) 

Xo = 2xi+xo + xi+2Xi, Yo^2yi + yo + yi + 2Yi. (3.9.13) 

We prove for instance the first one: 



W - / U3 = / r*ui = / ui by (3.9.9) 



ao Ja-2 J 0.3 J bo J bi J b-- 



ui by (3.9.1) 

Zl(p^ + p) + + p"^) ^ = -(Wl + Zl) = pZl - Zl. 



Note that the expressions in (3.9.11) automatically satisfy (3.9.6), (3.9.7). Also, the Rie- 
mann bilinear relations are satisfied as well. 

Moreover, they can be rearranged with the identities above to get some new relations 

Zi = -p^Zi, Wx^Zx. (3.9.14) 



Finally, the Riemann period matrix takes the following form 

( a h h h\ 

b c d d 

b d c d 

b d d c 



(3.9.15) 



in accordance with Fay-Accola Theorem 3.2; here we are also able to find explicitly the quan- 
tities a,b, c, d to be 



xiYq - Xpyi 
xm - xoyi ' 

2 Zi ^ 1 xqYi - Xiyg 

3 zi 3 xiyo - xoyi 



XiYi - Xiyi 

xiyo - Xoyi ' 
1 Zi 1 xqYi - Xiyo 
3 zi 3 xiyo - Xoyi ' 



(3.9.16) 
(3.9.17) 
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3.10. An ALTERNATIVE FORM FOR THE PERIOD MATRIX WHEN a = 

3.10 An alternative form for the period matrix when a = 

We know the expression of Fay's cyclic basis in the case a = (see the upper parts of Figures 
3.12, 3.13 ); in fact, this particular case was indeed our starting point in section 3.7. To this 
homology basis corresponds a period matrix with a particularly simple form, different than 
those found in the literature thus far (see [BE06] and references therein). In this section we 
give this alternative simple form. 



Firstly, we remark that in this case several simplifications occur: in particular, as pointed 
out in section 1.6.2^, all the integrals reduce to just eight, namely 

a /3 
Jdui=Ii, J dui ^ Ji, 1 = 1,..., 4. 



where Ai = (a, 0), A4 = (/3, 0) (we recall that the branchpoints are defined in (1.6.2)) and 
^ = 4, and it is understood that these integrals are computed on the first sheet. 
This is because, making use of the symmetries of the curve, we see that 

■yiiO.W'Bi) -fiiQ.W'Bi) 

U2 = P 2:2(0;), / U2 = p^ J2(a), 

J u!] = p%(5), J uO = p^J3(5), 

where again the integrals I^, Ji are computed on the iirst sheet. 

These integrals can be expressed in terms of hypergeometric functions as follows (cf. [BE06], 
p. 50i°): 

Mia 



Jiia) 



gr{iy (1 + 56)5 

47r^ a 



37f ; " 3^ (1 + 56)1 ^^^^3'3' 



, , 27r(52 /2 2 A 2n dP' ,12, 
12(a) - -^2^1 -,-;l;-a6 = — = 2Fi{-, --,1,1), 



^Due to a different choice of the ordering for the differentials, the formulae presented here agree with those in 
[BE06] up to a permutation of the subscripts of I, J. 
^"See the caveat in the previous footnote. 
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3.11. The Ercolani-Sinha conditions 



27r /2 2 1 



3^3 (1 + a 



~3' 3' 



J3(5) 



1 



2F1 



3' '3' 56 



27r a / 1 2 



3\/3 (1 



a" 3 



^3' 3' 



27r /II 1\ 27r 5 ,12 

3V3a V3 3 a'^ J 3v3 {l + a^)^ 3 3 

wltht = <5V(1 + a^)- 

One can check that the following relations hold 



~ Ji{a) May- 



Ma) + J2(5) = 0, 14(5) - J4(5) - 12(5). 



(3.10.1) 



Using the above expressions for the integrals, and relations (3.10.1), in the arc expansion 
(3.7.6) of the elements of the cyclic basis in the case a = 0, the matrices of periods simplify 
as follows 



/ 



^0 = 



{-ZR-i) p-iR 
l + 2p {2-R)p+l + R 
l-p {2-R)p+l + R 







-3i? + 3)p + 3 



-l + p {-2-R)p-l-2R 
-l-2p {-2-R)p-l-2R 



( 



J2 



y-p-2 {2-R)p+l + R p + 2 {-2- R) P-1-2R j 



\ 



Ji j 



I 



{-6R-3)p-3R + 3 



{-6R + 3) P-3R + 6 



\ 



l-p {-l-R)p+l-2R p + 2 {l-R)p + 2 + R 

-p-2 {-1-R)p+1-2R -l + p {l-R)p + 2 + R 
\^ l + 2p {-1-R)p+1-2R -l-2p (l-R)p + 2 + R j 



J2 



(3.10.2) 



We remark that transforming the period matrices above to the symmetric basis of Figure 
3.11 using the change of basis matrix (3.7.7), the period matrices obtained agree with those 
found in [BE06], eq. (7.36). 



3.11 The Ercolani-Sinha conditions 

As shown in [BE06], the curve X, with a = 0, satisfies the Ercolani-Sinha constraints for cer- 
tain values of /3 and 7, and hence describes indeed a monopole. 

In this section we express the Ercolani-Sinha constraints in the different cyclic basis, and 
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3.12. The vector of Riemann constants 



hence manage to reduce them as relations among periods on the genus 2 quotient curve X: 
this proves crucial for the solution, as will be shown in Chapter 5. 



The Ercolani-Sinha conditions here read 

n^A + ni^B = -2iy (0, 0, 0, 1). 

In solving the first and third of these conditions, using (3.9.8) one deduces that = rij and 
mj ~ rrij . Hence the remaining equations can be expressed as 

''^ +(TOo,3m) ' = -2^.(0,1), (3.11.1) 
xi yij \X, Y,j 

which are the Ercolani-Sinha equations on the reduced curve with n — tii, m — rm. Alterna- 
tively, in terms of the periods on the genus 2 curve X, one has 

(no, 3n)A + (3too, 3m)B = 6(^(0, 1), (3.11.2) 

Note that this also implies that 

U = Tr*U. (3.11.3) 
Also, for the Ercolani-Sinha vector of Theorem 1.7, one has 

c = 7r*(c), where c = noQo + 3rtOi + 3rno&o + 3mbi (3.11.4) 



One can see that condition (3.11.1) is satisfied at a = 0, 7 = ±5i/2 and 
2n — m 

R = , Hq = 5 72 — m, rriQ = —3n. 

m + n 

We remark that these are indeed the same values obtained in [BE06] , as one can easily check 
transforming the Ecolani-Sinha vector of Proposition 1.14 from one basis to the other using 
the matrix (3.7.7) (see Maple file change _bas i s _alphaO .mw, final part). 



3.12 The vector of Riemann constants 

Below we calculate the vector of Riemann constants for the genus 2 curve X 

= (a;3 + aa; + 7)2 +4/32. (3.12.1) 
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3.12. The vector of Riemann constants 



For applications we need the basepoint for the Abel map to be one of the points at infinity, 
say P+; but, in fact, it is easier for our calculations to choose one of the branchpoints, say Bi; 
then, to obtain back the vector of Riemann constants with respect to P+ we use the following 
relation 

/Cp+ =^p+(Bi)+/Cb,, (3.12.2) 
where Ap^ is the Abel map with basepoint P+. 



3.12.1 The vector /Cb, 

In the following we use the basis for the curve X as in Figure 2. The arc expansions for the 
cycles of eq. (3.8.3) can be equivalently expressed as follows 

ai =7i(2,l)+7i(l,6)+72(6,l)+72(l,2), 
bi = 71 (6, 5) + 7i(5, 4) + 72(4, 5) + 72(6, 6), 

ao - 71 (3, 4) + 72(4, 5) + 72(5, 6) + 71(6, 1) + 72(1, 6) + 71(6, 5) + 71(5, 4) + 72(4, 3), 
bo = 71 (3, 4) + 72(4, 5) + 71(5, 4) + 72(4, 3). 

To simplify the above expressions, we proceed as in Farkas and Kra ([FK80] VII. 1.2). Using 
that under the hyperelliptic involution J : {x, y) — > {x, —y) one has 

7fe+i(2j + 1, 2j) - - J7fe(2j, 2j + 1), 

which implies that 

Jfk{6,5) \J7fc(l,2) J-iki3A)J Jlk{lfi) \J7fc(2,3) Jik(4:,5}J 

Using these relations, one finds 
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3.12. The vector of Riemann constants 



From these one can solve to obtain the integrals between branchpoints as follows: 



7l(l,2) 



7l(2,3) 



71 (3,4) 



71 (4,5) 



i(e(0)+e(i)+^(0)) 



71 (5,6) 2 

Hence one can easily deduce the image under the Abel map, with basepoint Bi , of each 
branchpoints (their characteristics are also given) 



^B,(i?5)-^(e(°)+7r(0)) 







1 1 





1 

1 

1 

1 
1 

1 1 

1 



Note that there are only 2 (= g) odd characteristics, corresponding to B2 and Bq, and that 
the rest are even. 

Following an argument of Farkas and Kra ([FK80] VII.1.2), we recall that [^/ ] is even/odd if 
the parity of the characteristic [ ] is even/odd. Hence, by definition of theta with character- 
istic (A.2.3), the theta function vanishes at odd half periods. Since only g branchpoints have 
odd characteristic, we can conclude that the vector of Riemann constants takes the form 



1 

1 1 



(3.12.3) 
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3.12. The vector of Riemann constants 



3.12.2 The vector /Cp+ 

Once we have calculated the vector of Riemann constants with Bi as basepoint, making use 
of equation (3.12.2) we can change its basepoint; to simplify future calculations, we choose 
one of the point at infinity, say P+, as basepoint. 

Numerically, using the command AbelMap in Maple (developed by Deconinck and Patter- 
son^ ^ see [DP08]), one can compute the Abel map of Bi based at P+, obtaining 



Ap^Bi) 



=(1) 



r(l) 



i 



(3.12.4) 



Hence, the vector of Riemann constants with basepoint P+ is 



e(2) + -TT 



(1) 



r(2) 



1 1 

2 2 

1 1 

6 2 



(3.12.5) 



3.12.3 Thecasea = 

We recall that according to the results by Fay and Accola presented in section 3.4.2 (see 
[Fay73, Acc71]), the following relation between the vectors of Riemann constants on the two 
curves holds (cf. (3.4.9)) 



(3.12.6) 



n=0 



In the case where a = 0, we have that 3 / e A for any branchpoint, and moreover they are 

Si 

sent one to the other by a. Hence, Fay- Accola relation becomes in this case 



(3.12.7) 



Using Fay's expression (3.4.8) for e, we get 

'3-1 



e = TT 



2-3 



,0 - (1,0,0,0) = 0. 



In [BE06] (section 5.3) the vector of Riemann constants K.^ for the curve Xq is computed to 
be a half period, namely 



1111 
1111 



In view of the Fay-Accola relation (3.12.6), the expression (3.12.5) is indeed consistent with 
this result. 



i^Thanks to Matt Patterson for providing additional documentation, and help in several occasions. 
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3.13 Theta functions 



For completeness we observe also that applying Fay-Accola Theorem 3.2 we obtain a sim- 
ple expression for the theta function on X associated to the cyclic basis of section 3.7. In 
particular we have 

^(3zo,zi,zi,zi;'r) 



nLo^ 




- 

3 " 



(3.13.1) 



(2^o,2i;r) 



where as in Fay-Accola theorem co(f ) is a constant. 



Hence, the genus 4 theta function Q can be expressed as a product of hyperelliptic 
ones. This is particularly important when attempting to solve Hitchin constraints 113 in the 
Ercolani-Sinha formulation (1.5.19): despite this semplification, this proves to be a nontrivial 
task, and will be examined elsewhere. Here we only remark that, provided that the Ercolani- 
Sinha constraints (1.5.20) and (1.5.19) are satisfied, eq. (3.13.1) means that the monopole 
solutions can be expressed in terms of hyperelliptic functions. 
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CHAPTER 

4 

An interlude: elliptic subcovers 



Quotienting the monopole curve X by the action of the symmetry group C3 
we have obtained a genus 2 curve X, and are hence able to express several 
objects on X in terms of objects on X, as shown in Chapter 3. This provides 
a remarkable simplification of our task of imposing the Ercolani-Sinha con- 
straints, allowing us to work in genus 2 rather than 4. Therefore, before tack- 
ling this task, we can investigate whether a further simplification is possible, 
namely, whether the curve X also covers elliptic curve(s): if this is the case, 
obviously the problem would simplify much more. This is examined in this 
short chapter, making use of invariant theory. 



4. 1 Invariants and elliptic subcovers 

We want to investigate whether the genus 2 curve X admits a maximal covering ip ■ X ^ E 
to an elliptic curve E, or, in other words, whether X admits (degree n) elliptic subcovers. 
Here maximal means that the map ip : X E does not factor over an unramified cover of 
E. 

We remark here that such degree n maximal elliptic subcovers always occur in pairs, [E, E'); 
there is an isogeny of degree between the Jacobian of X, J{X) and the product E x £": 
one says that X has (n, n) split Jacobian. 

Here we only consider the case of a degree 2 elliptic subcover, as this is indeed the case 
for a = 0, as shown in [BE06]; for continuity, one expects then that also in the case a 7^ 0, if 
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4.1. Invariants and elliptic subcovers 



there is an elliptic cover at all, this should be of degree 2. 

Before analysing this case, let us introduce the invariants in some generality. 
4. 1 . 1 Igusa invariants for a genus 2 curve 

The J -invariants were introduced first by Clebsh and Bolza [Cle72, Bol87], and generalised to 
characteristic different from 2 by Igusa [Igu60] ; many results used here are due to Shaska (see 
[Sha05] for a review). These invariants are used to study the properties of moduli of genus 2 
curves in terms of the associated sextic. 

Recall that every genus 2 curve can be represented as a sextic^ given in general form as: 

OqX^ + aix^ + a2X^ + a^^x^ + a^x'^ + a^^x + ag — 
aoix - ai){x - a2)ix - a3){x - a4){x - a5){x - ag) 

One can define an integral (or relative) invariant J, as a continuous function on the space of 
sextics if it is invariant under the natural action of the group SL{2). One can explicitly build 
these relative invariants ([Igu60], section 2); they are homogeneous polynomials of weights 
2, 4, 6 and 10, respectively given as 

i<j <k<l<7n<n 

k = (4) X! [oii- ajY{aj - akY{ak~ aiY{ai- a„^Y{a,rl- an)'^{an~ aif, 

2<j<A;</<m<n 

J6 = ao X! [ai- ajY{aj - akf{ak- aif{ai~ a„if{a,n- ctnYian- aif 

i<j <k<l<m<n 

with Qfj roots of the sextic. Here the indices run from 1 to 6, so that j2 is the sum of 15 ele- 
ments, ji of 10 and jg of 60. 

One can then introduce absolute invariants, i.e. continuous functions on the set of hy- 
perelliptic curves, assuming the same value on birationally equivalent curves (see [Igu60] , 
section 6). These absolute invariants can be expressed in terms of relative invariants as fol- 
lows 

zi-144^, i2 = -1728 ^^2^^^^^, Z3 = 486 ^. (4.1.1) 

.?2 J2 J2 

' This is due to the fact that every genus 2 curve has a representation as a double cover of with six ramification 
points Xi. This establishes a bijection between isomorphism classes of such curves and binary sextics, as every 
unordered 6-tuple of points can be described by such a sextic. 
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It is a necessary and sufficient condition for two genus 2 curves (with-^ j2 ^ 0) to be iso- 
morphic that their absolute invariants are the same. 

4. 1 .2 Elliptic subcovers of degree 2 

We are now in a position to give conditions on the parameters of the curve X under which 
it is a degree 2 cover of elliptic curves, making use of invariant theory. One can show (see 
[Sha05] , section 5) the following characterisation. 

Proposition 4.1. A (non-singular) genus 2 curve has a degree 2 elliptic subcover if and only if 
the automorphic form xw vanishes, where 

X30 = - 19245600 iio^i4«2^ + 507384000 nn^i4^i2 + 77436^014^*2'' - 81^2^*6'' - «2^j4^ + 384^4^*6 

- 78 ^2^*4^ + 31104 ig^ - 3499200 110^2*6^ + 104976000 iio^i2^i6 + 972 ^1012^*4^ + 8748 iioi2\^ 
+ 1332i2^«4''i6 - 8910i2^j4^«6^ - 1728i4^«2^j6 + 6048i4''i2j6^ + 29376 i2^i4^j6^ - 47952^2*4*6^ 
+ 41472^10^4^ - 236196^10^^2^ + 159i4^i2^ - 80i4^»2 - 125971200000 iio^ - 592272 Hoi4''«2^ 

- 870912 iioi2^i4^ie + 3090960 ^1012^^4*6^ - 5832 i 10*2^ ^4 + 4743360 Hoi4^*2*6 - 6912 ^42^6^ 
+ 12^2^14^*6 - 54^2^*4^*6^ + 108i2'^i4jj6^ - 2099520000 jio^*4*6 - 9331200 iio*4^*6^ 

Proof. If a non-singular genus 2 curve has a degree two elliptic subcover, then it can be put 
in the form 

F2 = _ s,X^ + 52^2 - 1, 

with 27 - 18siS2 - sjsj + 4sf + 4s^ 7^ (see [SV04] for a proof). Note that this determines 
X up to a coordinate change in the subgroup of the automorphism group of E generated by 
(ri, T2), where ti : x 1/x and T2 : x ^ p^x, where is a 6-th root of unity. Invariants for 
this action are: 

u — S1S2, V = sl + si- 
In terms of these new variables, the relative invariants can be written as 

j2 = 240 + 16m, 

j4 48w + + 1620 - 504u, 

je = -2066w + 96w - 424^^ + 24-1*^ + 160™ + 119880, 
jio = 64(27 - 18u -u^ + 4w)2. 

Using these expressions, the absolute invariants *i, 12, *3 can be expressed as rational func- 
tions of u and v, for j2 7^ 0. Eliminating these variables and using eq. (4.1.1) give the required 

2if j2 = one can define new absolute invariants. 
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expression for xao- 

For more details about these calculations we refer to the Maple file invariants4 . mw. □ 

4.2 Igusa invariants for the curve X 

Let us now apply the previous considerations to the curve X of equation (3.1.8), where we 
can rescale to set /3 = 1 (see section 5.3): 

y2 ^ (a;3 _|_ ax + 7)^ + 4. 

The Igusa invariants can be expressed in terms of the parameters a, 7 in the following way 
(see Maple file invariants4.mw): 

ia =- 32 a^- 216 7^- 960, 

=64a^ + 864 0^7^ - 2496 + 2916 7"^ + 18144 7^ + 25920, 
= - 157464 76 - 17496007^ - 6397056 7^ - 7672320 ~ 69984 0=^7'' 
- 285120 0^7^ - 10368 a^-f'^ + 648960 ~ 3840 - 512 a^ 
iio = - 47775744 - 23887872 7^ - 884736 a^j^ + 3538944 - 2985984 7" - 65536 a^. 

Substituting these in the expression above for xso. the condition of vanishing of Proposition 
4.1 reads 

7^0^ (-8000000 + 1146000 - 53088 + 784 - 6480000 7^ + 327240 a^-f'^ + 7128 -f'^a'^) 
-1749600 7^^ - 10935 7^a^ - 1574647^)^ = 0. 

From this expression it is evident that xso vanishes identically when a — 0, hence the curve X 
with a = always has an elliptic subcover; this is in agreement with the results from [BE06] 
(lemma 7.2, also summarised earlier in section 1.6.3). 

As for the case a ^ 0, analysing the above equation one sees that there are real solutions for 
some values of a and 7, as plotted in Figure 4.1, where 7 is given in function of a. In contrast 
with the case a = 0, then, the curve X does not admit an elliptic subcover for all values of 
the parameters, but just for some special ones. 

We come back to this graph in the last chapter, when we analyse solutions to the Ercolani- 
Sinha constraints (3.11.2). In particular, in section 5.3.1 we find a curve of solutions to 
(3.11.2) in the 0,7 plane. We can hence examine the intersection of this curve with that 
of Figure 4.1, in order to find those Riemann surfaces solutions of the Ercolani-Sinha that 
also admit an elliptic subcover: these are particularly interesting as they correspond to 
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y 




Figure 4.1: Solution for elliptic subcovers of X 



monopoles solutions expressible in terms of elliptic functions. 
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CHAPTER 

5 

The AGM method 



Gauss developed a rather efficient method for computing elliptic integrals, 
namely the Arithmetic-Geometric Mean, or AGM, first introduced by La- 
grange, which was then generalised to a class of genus 2 curves by Richelot 
and Humbert. Here we use this method for the curve under consideration, of 
eq. (3.1.8); the AGM method cannot be applied as it is to this class of curves, 
but needs to be modified. This is done in this chapter, after a short introduc- 
tion to the AGM method in general; the findings are then applied to solve the 
Ercolani-Sinha constraints numerically. 

We are now in a position to solve explicitly the Ercolani-Sinha constraints, in the form 
given in eq. (3.11.2), namely in terms of the genus 2 periods: 

(riQ, 3n)A + (3too: 3m)S = 6i/(0, 1), 

In order to do this, we make use of the so called Arithmetic Geometric Mean (AGM in short) 
method to calculate the periods. The AGM method permits one to calculate, in a rather effi- 
cient way, all the integrals between branchpoints, and hence, from eqs. (3.8.3), the periods: 
in particular, it is possible to check the Ercolani-Sinha constraints for given values of the pa- 
rameters^ a and 7. As we know that the Ercolani-Sinha constraints have a solution for a = 
and 7 = ±5\/2 (tetrahedral monopole, see [HMM95, BE06] and Theorem 1.18), it is reason- 
able to look for a solution in the (a, 7) plane in a neighbourhood of these points. This is done 
numerically: starting from a — 0,7 = 5a/2 we let a and 7 vary infinitesimally, calculate the 
Ercolani-Sinha constraints for every such pair and look for a 7 such that these constraints 
'We can rescale so as to have 13 = 1, see section 5.3. 
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are still satisfied, with a, 7 e M. 

The results obtained are described in section 5.3.1: in view of this, we give some details about 
the AGM method. 



5.1 The AGM method 

The origin of the AGM method dates back to Lagrange, but it was Gauss who truly initiated 
its investigation, and in fact a large part of what is known today seems to be due (or at least 
known) to him. In particular, Gauss examined the elliptic case in detail, while the generali- 
sation to the hyperelliptic case was due to Richelot, even though it seems that it was known 
to Gauss already (for historical notes see e.g. [Cox84]). 

Here we describe briefly the elliptic case, and then, in some more detail, the hyperelliptic 
case; for a more extensive treatment, we refer to [Cox84, BM88] and references therein. 



5.1.1 AGM: the elliptic case 

Let a, b be positive real numbers: their arithmetic-geometric mean, denoted A/(a, b), is the 
common limit of the sequences defined as follows: 

ao = a, &o = b, 

fln+l = " ^ " , bn+l ^ V anbn- (5.1.1) 

These two sequences satisfy the following property 

a> ai> . . . > an> a„+i > . . . > 6„+i > b„ > . . . > bi > bo, 

which can be proven by induction, and ensures the existence of a limit. From the observation 
that bn+i > 6„, we have 

O'n+l ^ ^n+l < On+l — &n = 2 ^ 

So, the sequences a„ and &„ converge to a common limit, their arithmetic geometric mean 

lim a„ = lim &„ = M{a,b). 

Moreover, one has that 

0<a„-bn <2""(a-6), 
which ensures a fast convergence (which is relevant in the present work). 
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Lagrange and Gauss introduced the Arithmetic-Geometric Mean in connection with el- 
liptic integrals, obtaining the following result: 

Theorem 5.1. AGM [Gau99] Leta, b e M+, andM(a, b) their arithmetic geometric mean; then: 



^Ja? cos2 ^ + 62 sin2 2M(a,6) 

For the proof see e.g. [BM88]. 
Using this result, all elliptic integrals of the form /^^ ^^^^ can be expressed in terms of an 



appropriate arithmetic geometric mean (see [BM88]), after various changes of variables. 



5. 1 .2 The genus 2 case: Richelot and Humbert 

The above theory of Gauss' Arithmetic Geometric Mean was extended to the hyperelliptic 
case by Richelot, in [Ric37, Ric36] . This section follows the modern exposition of the work by 
Richelot given by Bost and Mestre in [BM88] . These authors remark, how Richelot's "changes 
of coordinates" is in fact a correspondence (cf. definition 5.2 below); they follow the con- 
struction given by Humbert, who in [HumOl] re-interprets Richelot's findings in terms of 
duplications formulae of 2-variable theta functions, i.e. isogenies (of type (2,2) ) on Abelian 
surfaces, hence generalising and simplifying Richelot's algorithm. 



We remark that the Richelot-Humbert construction is only given for the case where the 
genus 2 curve, represented as a two-sheeted cover of with six branchpoints, has all real 
branchpoints: in fact, it does not apply in a straightforward fashion to the case of a genus 2 
curve with complex roots. In the next section we argue a possible generalisation in this di- 
rection. 



Consider the genus 2 curve C" 

+ P{x)Q{x)R{x)=Q, 

(5.1.3) 

P{x) = {x-a){x~a'), Q{x)^{x-b){x-b'), R{x) ^ {x - c){x - d), 
where the real roots a, a', 6, 6', c, c' are ordered as 

a < a' < 6 < 6' < c < c'. 



2' 



The proof is based on tlie following change of variables 

. 2a sin tfi' 

sin <p = = (5.1.2) 

a + b + {a~b) sin^ (f>' 

and makes use of the identities (a^ cos^ t/i + b^sin^ (f>)~^^^d4> = (a^ cos^ 0' + 6jSin^ (7!>')^^/^d(7!>', proven by Jacobi. 
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Associate then to the triple of (real) polynomials (P, Q,R), another triple (C/, V, W), defined 
by 

U{x)^[Q{x),R{x)], V{x) = [R{x),P{x)], Wix) ^ [P{x),Q{x)]. (5.1.4) 

df{x) dg{x) 
where [f,g] := -^5(2^) - "i^^fi^)- 

In general, given two non-proportional degree 2 polynomials Pi,P2, the roots r, r' of 
[Pi, P2] can be obtained geometrically from the roots pj,pj of P^. If C is a conic, and g : — > C 
a unicursal parametrization, then g{r),g{r') are those points where the tangent to C passing 
through the intersection of the lines g{pi)g{p[) and 5(^2)5(^2) meets the conic C itself. 
As for the roots of the polynomials U, V, W, they are all real; denoting u < u',v < v',w < w' 
the roots of U, V, W respectively, the following relation between all the roots holds 

a<v<w<a'<b<w'<u<b'<c<u'<v'<c'. (5.1.5) 

These inequalities can be proven algebraically (see section 5.1.3, where we give explicit 
expressions for the roots of U, V, W); in view of the above interpretation of the roots of the 
bracket [Pi, P2], a geometric description is given in Figure 5.1. 

Therefore, this is a situation similar to Gauss' AGM case: to each pair of branchpoints one 




Figure 5.1: Roots of P, Q, R and U, V, W 
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can associate another pair of points which are closer than the initial ones, and we expect a 
relation between the integrals of corresponding pairs. Iterating this process, we can repeat 
the considerations of section 5.1.1 for every pair to argue the existence of a limit, and hence 
obtain an expression for the integrals in terms of these limits. The relation between integrals 
suggested and proven by Humbert is 



=d.T = 2VA r ^^"'^ =dx, 

) Jy J-U(x)V(x)W(x) 



^-P{x)Q{x)R{x) Jv ^-Uix)V{x)W{x) 

and similarly for the integrals between the other pairs of branchpoints, where^ A is the de- 
terminant of the matrix whose entries are the coefficients of P, Q, i? in the basis {l,x,x'^), and 
S{x) is a polynomial of degree at most 1. 

There is, however, an important element of difference between the elliptic and the hy- 
perelliptic case: while in the standard Gauss construction the change of variables can be 
interpreted as a map from one curve to the other, in this case it is given not by a map, but by 
a correspondence. 

Definition 5.2. [GH78] A correspondence T : E — > S' of degree d associates to every point 
P a divisor T{P) of degreed in S', varying holomorphically with P. 

A correspondence can be presented both as a holomorphic map S S'^"*' from S to the 
d-th symmetric product of S', or equivalently by its "curve of correspondence": 

D = {(P,P') : P' e r(P)} c C X C". 

In the case of the Humbert construction, the two curves are C, with equation + 
P{x)Q{x)R{x) = 0, and & , with equation Ay'^ + U{x')V{x')W{x') = 0. The correspon- 
dence between C and C' , which Humbert considers in [HumOl], is of degree 2, and is given 
by the curve Z in C x C" of equations 

P{x)U{x')+Q{x)V{x')=Q, 

(5.1.7) 

yy' ^ P{x)U{x'){x-x'). 

In analogy with the pull-back of a map, one can introduce also for a correspondence a 
linear map 8z : VL^(C') il^ (C). Then eq. (5.1.6) can be interpreted as the relation between 
differentials 

^(-')d.'U«d., (5.1.8) 



y J y 

together with an analysis of the image of the path joining a, a' (resp. b, b' or c, c') under the 



^The reason for the presence of this constant is clarified in the proof given in the next subsection. 
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correspondence Z. 



Furthermore, eq. (5.1.6) can be interpreted as an identity in J{C'), the Jacobian of C": it 
states that a half period in C is sent by the correspondence Z to a period in C', or alterna- 
tively, the image of a half period is zero in the Jacobian. This is a slightly weaker formulation, 
which is explored further in Remark 5.1.2. 



The proof of the identity (5.1.6) is given in section 5.1.4: we state first a version of the 
Arithmetic- Geometric Mean for genus 2 curves, bases on the remarks above, and on eq. 
(5.1.6). 



5. 1 .3 The AGM method for genus 2 curves 

Consider a genus 2 curve as in eqs. (5.1.3), (5.1.4). Define 6 sequences (a„), (a^), (6„), (6^), 
(c„), (c^) recursively by the conditions: 

• ao = a, Oq = a', fop ^ b, 6q = 6', Cq ~ c, c'q — c'; 

• a„+i, a^,^!, 6„+i, c„+i, c;,+i are roots of t/„y„VF„, ordered as follows 

a«+i < a^i < bn+i < Ki+i < Cn+i < c^+i, (5.1.9) 

where, for every n, 

Pn{x) ^ (x - an){x - a'J, Qn{x) ^ {x - bn)(x - b'^), Rn{x) ^ (x - Cn){x - c'J, 

U„ix)^[Qnix),R„ix)], Vn{x) = [Rnix),P„ix)], P^„(x) = [P„(x), Q„(x)]. 



Bost and Mestre give in [BM88] an explicit expression for these sequences: 







— Bn 


Cn + 


c'n ^ In 


-< 


bnb'n 


- a„a^ 


+ Cn 


bn + 


b'n - a„ 


- < 




- bnb'n - 


'rAn 



bn + b' - an - a' ' 



n n' — h h' — A 

''"+1 = ^r-VT, — : IT' ^.+1 = „ +r' ~h -h' ' • 

_ _ / Cj7,c^, — djid,^ -\- Bji 



a^ 



where 



An = \^ (bn 


- c„)(6„ 


-dnWn 


- c„)(6;, - c'„), 


Bn — \J (c„ 




-<)« 





C« = (a„ - 6„)(a„ ~ &;^)(a;, - bn){a'^ ~ b'J. 
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These can be derived finding the roots for ?7„, 14, Wn as foilows'* 





~ bnb'j^ ^ An 


c„ + 


c' -b -b' ' 

^n 




— a„aj^ =F Bn 


Cn + 


c'n ^ ~ o!n 


bnK, 


— ClnO^n ~F ^n 


bn + 


b'n ~ ~ O-'n 



(5.1.11) 



and ordering them as in eq. (5.1.9). One can see directiy from the expressions above that 

Vn < Wn < w'n < Un < u'^ < v'^ (cf. also cq. (5.1.5)); hence we set 

fln+l = Vn, = Wn, bn+1 = w'^, = Wn, C„+i = u'^, c'^^^ = U^; (5.1.12) 

the expressions in (5.1.10) then follow. 



Theorem 5.3 [Richelot [Ric36], Bost and Mestre [BM88]]. With the above definitions, the 
sequences {an), (a^), (&„), {b'^), (c„), (c^) converge painvise to common limits 

lim a„ = lim = a = M{a,a'), 
lim bn = lim b'^^ P = M{b, b'), 

n— >-oo n— >-oo 

lim Cn — lim cj^ — 7 = Af(c, c'). 



Furthermore, for any polynomial S{x) of degree at most one, the following relations hold: 

r' s{x)dx 

i[a,a ) = / — - = irl 



nb, b') ^ f , "I'l^ ^, ^ = .T ,^ (5.1.13) 







-7)' 








-7)' 


, ^(7) 





^-P(x)Q{x)R{x) (7-a)(7-/3)' 



where 



T=ftt„, ^„ = =. (5.1.14) 

„=o V + b'^-an- a'^){cn + d^-bn- b'J{cn + d^-~an- 



The proof of convergence for the sequences (a„), {a'^), (bn), {b'n), {cn), is similar, for 
every pair, to that for the elliptic case given in section 5.1.1. 



'*To derive this we only use eqs. (5.1.6) for the correspondence. For Cn,c'„ eq. (5.1.16) is also needed. 
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It follows from a direct calculation that 

P^{x) = [Q{x), R{x)]{-a ~ a' + c + c') = U{x){a + a' - c- c'), 
Qj^{x) ^ [P{x), Q{x)]ia + a' -b-b') = Wix){a + a' - b - b'), 
Ri{x) = [R{x),P{x)\{b^b' -c-d)^V{x){b + b' ~c-d). 

Let us focus on the fist pair a, a': the other cases can be dealt with in the same fashion. 
Using (5.1.6), the following holds: 



^ ' --dx = 2VA„ / ^ ' dx 



a„ PnQnRn J a„ + i \l ^ \Pm Qn\ [Qrii ^n] [^n i Pn\ 



Six) 



(5.1.15) 



Ja^^i \/ — Pn+lQn+lRn+l 



It follows that 



'a \/—PQR \k=0 J •^'^"+1 y~Pn+lQn+lRn+l 

Taking the limit for n — >^ oo 

- (n '0 r , , t' n - ^ (n 5 / , . . " 

S'fa;) 

= TtTZ ReSa;=a ^= , 

\/ ^ Pn+lQ7i+lRn+l 

withTasineq. (5.1.14), yields eqs. (5.1.13). 



Remark 5.1.1. The expressions for the integrals between branchpoints we have given are 
slightly different from those in [BM88] : in particular, the integral between b and b' has op- 
posite sign. This depends on a different choice of conventions. Indeed, Bost and Mestre, 
in the note 2, p. 51 of [BM88], claim that they want to recover the "classical identity" 
la - h + Ic = 0: with our choice of convention for sheets, the relation between integrals 
become la + lb + Ic = 0, which follows from the fact that the integral around a cycle encir- 
cling all the cuts, oriented so that the upper arc goes from negative to positive real values, is 
zero. 

5. 1 .4 The proof of a fundamental identity 

For the proof of eq. (5.1.6) we follow [BM88]: in particular, the authors foUow the lines sug- 
gested within a problem presented at the entrance examination to the faculty of Physics in 
the Ecole Normale Superieure, Paris, in 1988. The object of the problem was indeed to prove 
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Richelot theorem: a number of hints were given and the proof was divided into several steps; 
here we follow these steps, providing some additional details. 

For every pair (x, z) the following fundamental equality is valid 



Pix)U{z) + Qix)Viz) + Rix)Wiz) + {x- z)^A{P,Q,R) = 0; 



(5.1.16) 



this follows from evaluating the determinant 



det 



/ P{x) P'{x) P{z) ^ 

Q{x) Q'{x) Q{z) 
y R{x) R'{x) R{z) y 



Recalling eq. (5.1.7), we introduce the following polynomial of 2 variables (x, z), 



F{x, z) = P{x)U{z) + Q[x)V{z). 



(5.1.17) 



Its zeros define two functions, zi(x) and 2:2 (2^)- which are the (z-coordinates of the) images of 
the correspondence Z of eq. (5.1.7) in C 



F{x,z^{x)) =0, i = l,2 y X. 



(5.1.18) 



Following^ [BM88], define, for x e (a, a'), 



2/2(3;) 

One can compute 2/1,2(2;)^. to get^ 

2/1,2(2;) 



y{x) = ^-P{x)Q{x)R{x), 

P{x)U{z^{x))(x-zi{x)) 



y{x) 

P{X)U{Z2{X)){X-Z2{X)) 



2 _ 1 P{x)U{zi,2ix)fW{zi^2{x)) 



A 



Qix) 



A preliminary lemma can be proven. 



(5.1.19) 



(5.1.20) 



^We point out that Bost and Mestre use an absolute value under the square root for y{x). 
^Using the above expression for y{x) and (5.1.16) 



2 _ P(x)C/2(2i(x))(x-2i(x))2 



Q{x)R{x) 

1 P(x)U^(zi{x))[P(x)U{ziix)) + Qix)V{zi{x)) + R(x)W{zi(x)))] 



A Q{x)R{x) 
Recalling that, from (5.1.18), P(x)U{zi{x)) + Q{x)V{zi{x)) = 0, the result followfs. 
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Lemma 5.4. The following relations are valid 

yi{x) y2{x) y{xy 

Zi{x)Zi{x) Z2{x)z'2{x) X 

yi{x) y2{x) y{x)' 



(5.1.21) 
(5.1.22) 



Proof. Firstly, using 



U{z) = xi(z - ci)(z - c'l), 



with xi = 6 + 5' - c - c', one has 

yi{x) y2{x) 

y{x) ( z[{x) 



P{x) \ U{z,{x){x - zi(.t))) C/(z2(.t)(.t - Z2{X))) 
y{x) ( z[{x) 



z'^ix) 



Pix)xi { {zi{x) - ci){zi{x) - c[){x - zx(x)) {z2{x) - ci){z2{x) - c!^)(x - Z2(.x)) 
This expression can be further reduced to symmetric functions of zi {x) and zi (x), 

Z[{x) ^ z'2{x) 

z[{x) 



yi{x) 2/2(2;) 

_ y{x) 1 



1 



P(x)xi c — c' [ X — cJ 

y{x) 1 

P{x)xi {x — d){x — c) 



^ z'2{x) 

Zi{x) — d Z2{x) ^ d 

Z[{x) Z'2{x) 



Zl{x) — X Z2{x) — X 



y{x) [ i_ 

P{x)xi I X — d [dx 

_ yjx) 1 

P{x)xi {x — d){x — c) 
The symmetric functions 



ln[(zi(x)-c')(z2(x)-c')] 



1 



1 

X — c 



Z'ljx) ^ Z'2{x) 
Zi{x) — C Z2{x) — C 



dx 



ln[(zi(a;) — c){z2{x) — d)] 



—A + ^4 + H{zi{x) - x){z2{x) - x)] 

zi{x) — X Z2{x) ~ X dx 



Si{x) = Zi{x) + Z2{x), S2{x) = Zi{x)z2{x) 



are computed from the equality 



F{x, z) = 4)o{x){z - zi{x)){z - Z2{x)) 



via Maple. We do not report the expressions for si and S2 here as they are not very enlighten- 
ing, but using them the direct derivation of (5.1.21) is immediate. 
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Once (5.1.21) is proven, (5.1.22) follows. Indeed, we have 

Zi{x)z[{x) ^ Z2{x)z'2{x) 

yi{x) y2{x) 

y(x) ( z[{x) 4 (a;) 1 y{x) ( z[{x) z'^ix) 



Pix) {U{z,{x){x-z,{x))) U{z2{x){x~Z2{x)))l P{x) {U{z,{x)) U{Z2{X)) 

Using eq. (5.1.21) the first term in the right hand side simplifies to x/y{x). The second term 
vanishes, as one has 

z[{x) ^ z'^ix) d / (zi(a;) - ci)(z2(a:) - ci) ^. (5 123) 



U{zi{x)) U{z2{x)) dx \izi{x) - C2){Z2{X) - C2) J ' 

but the argument of the logarithm is a constant, namely 



b'cJ - c2 - c'2 - 2bb' + bc + be' + 2y/{b' - c'){b' - c){b - c'){b - c) 
-b'd + c2 + c'2 + 266' - 6c - bd + 2^(6' - d){b' - c){b - c'){b - c) 

Therefore eq. (5.1.22) holds. □ 

Adding (5.1.21) and (5.1.22), multiplied by arbitrary coefficients, one obtains the equality 

S{zi{x))z[{x) ^ S{z2{x))z'2{x) ^ S(x)_ ^ 
yi{x) y2{x) y{x)' 

where we remark that 5 is a polynomial of degree at most one''. 

Now we are in a position to prove (5.1.8). To do this we note that the first term on the left 
hand side of (5.1.24) can be transformed, by using (5.1.20) and then (5.1.18), as follows: 

o, I ^^^'lix) ^, , z[{x)^/K./Q{x) 
S{zi{x)) — -- = S[zi[x)) 



yi(^) "U{z,{x)),/P(x)^W{zi{x)) 

(5.1.25) 

Q( I \\ 2i(a:)\/A 

" ^-U{,z^(x))V{zr(x))W{z^{x))' 

and similarly for the second term. Thus, combining eqs. (5.1.24) and (5.1.25), we obtain eq. 
(5.1.8). 



5.1.5 Concluding remarks 

Remark 5.1.2. We point out that thus far we have only proven eq. (5.1.8); in order to obtain 
eq. (5.1.6), we need to integrate eq. (5.1.8) on the appropriate contour: on the left hand side, 
this is just a straight line connecting the two branchpoints, on the right hand side one needs 
to examine the image of the previous contour via the correspondence (5.1.7). 

'^In fact, one cannot find an analogue of eqs. (5.1.21) and (5.1.22) for higher degrees. 
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Note that this is not examined in detail by Bost and Mestre in [BM88] , although they do men- 
tion this issue. We are not able to give an abstract description of this in general, nevertheless 
we propose some progress in this direction. 

Consider the lift of the correspondence Z in eq. (5.1.7) as a map between Jacobians, 
namely C : J[G) J{C'), and denote Ba the branchpoint with z-coordinate a, so Ba = 
(a, 0), and similarly for the others. Then the kernel of this map is 

kerC - {Ba. - Ba, Bf,, - Bt, B,, - B, ,0}. 

This can be understood analysing the images of these branchpoints under the correspon- 
dence. 

For the first and second pair, i.e. a, a' and 6, b', using the notation of eq. (5.1.18) for the z- 
coordinates of the images via the correspondence, one has 

zi(a) = zi(a'), Z2{a) = Z2{a'), (5.1.26) 

and similarly for the second pair. It can be seen immediately from the first of (5.1.7) that 
roots of P correspond to roots of V; writing explicitly the 2 solutions'* of (5.1.7) for X shows 
that we are indeed in the case of eq. (5.1.26). 

Viewing the Jacobian as the lattice of periods, the image of the divisor Ba'-Ba (resp. B^' - Bt) 
is a point in the lattice. 

As for the third pair of branchpoints, Bc,Bc', recalling relation (5.1.16), then the first of 
eqs. (5.1.7) corresponds to R{x)W{z) + {x - z)^A(F, Q, R) = 0, hence, when Xc (resp. c'), R 
vanishes, and ^ 

zi,2(c) = (c,ri,2) =: Pl.2, ^i,2(c') = (c',fi,2) =: Pl,2. 
The fact that Be' - Be belongs to the kernel follows from the fact that 

C{Be' - Be) =Pl+P2-Pl-P2= div ( = Q 

\w — c J 

Note that the above considerations do not provide any kind of insight about the periods 
corresponding to the images of the above divisors: indeed, they involve objects in the Jaco- 

^Note that we cannot get Y from the second line of eq. (5.1.7) in this case, as we have 0, the i/-coordinates of 
the branchpoints, multiplying Y, the unknown that we need to determine: in this case, though, the two images 
of a branchpoint in the first two pairs are also branchpoints for the curve C' , so their Y-coordinates are uniquely 
determined. 

^Note that we still cannot determine the y-coordinates of the images, as for the other 2 pairs, but this time these 
images are not branchpoints of C", so solving in the equation of C we get 2 solutions, which we can take to be Yi 
and Y2 (resp. Y\ and 12). 
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bians, so by definition modulo periods. However, as the kernel has a group structure, there 
is a relation among those images (for instance, it suffices to find the images of Ba' - Ba and 
i?b' - Bb to obtain immediately that of Be' - Be). 

Remark 5.1.3. The integrals between other pairs of branchpoints, e.g. a' and b, can still be 
calculated using the extended AGM algorithm. This is achieved using an appropriate trans- 
formation: for example, the map 

f:x^ ^ J (5.1.27) 

2x — a ~ a' 

sends the roots a, a', b, b' , c, c' to 6 real numbers /(a) < /(c') < /(c) < f{b') < f{b) < f{a'), 
for which the above algorithm can be used. 

5.2 Generalisation to the genus 2 case with complex conju- 
gate roots 

We study here the possibility to generalise the AGM method to the case where the branch- 
points do not all lie on the real axis, as was the case in AGM analysis, but the polynomials 
P,Q,R axe still real: this corresponds to the three pairs of branchpoints being 

a' = a, b' — 6, c' = c, 

ordered as follows 

3?(a) = 5R(a') < 3?(6) = 5R(6') < 3?(c) = 3?(c'); 

in the following we also take 3(a) < 0, 3(5) < 0, 3(c) < 0, for definiteness. 
This is the case for the quotient monopole curve of eq. (3.1.8): recalling eq. (3.1.9) for the 
branchpoints, we see that indeed they can be split in complex conjugate pairs as above (cf. 
eq. (3.1.10)). 

We remark that all the polynomial relations found in the previous sections can still be re- 
covered in the case of complex conjugate branchpoints: hence, the relation (5.1.8) between 
the differentials on C and C" still holds true for this case. The difference lies, in fact, in the 
images, under the correspondence Z of eq. (5.1.7), of the branchpoints and the paths con- 
necting the branchpoints. 

In this case, since the initial branchpoints are complex, a relation analogous to (5.1.5) 
cannot be written; nevertheless, the roots of Uq.Vq, Wo are real, as can be seen by consider- 
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ing their explicit expressions in eqs. (5.1.11), and can hence be ordered. In contrast with the 
real case, though, this ordering is not unique: in the real case the ordering of u, u' , v, v', w, w' 
depends only on the relative ordering of a, a! , b, b', c, c' on the real line; in the complex case, 
instead, it also depends on their imaginary parts, as can be seen again from eqs. (5.1.11). 
Depending on the ordering of u, u' , v, v', w, w', equation (5.1.15) relating the integrals be- 
tween the three pairs of branchpoints on C and C" needs to be modified appropriately. We 
have considered in detail the case of the monopole quotient curve X: other cases can be 
studied using the same method. 

5.2. 1 The AGM method for the quotient monopole curve 

Consider the quotient monopole curve X of eq. (3.1.8), namely 

{x-^ + ax + "if + 4/3^ 
If we order the branchpoints as in section 3.1.1, we have: 

a^Bi,a'^Bi- b ^ B5, b' ^ B2; c^Sg, c' = Bi (5.2.1) 

We calculate u, u', v, v', w, w' via eqs. (5.1.11), and examining their relative ordering we find 
the following two cases: 

case 1 : a> 0, v < w < w' < u < u' < v'; (5.2.2) 
case2 : a < 0, u < v < w < u' < v' < w' . (5.2.3) 

Case 1 is exactly the same as the situation considered in [BM88]: in particular, eq. (5.1.12) 
still holds, and hence the sequences a„, a^, 6„, b'„, Cn, c'^ are still given by eqs. (5.1.10). There- 
fore, the first equality in eq. (5.1.15) holds, in view of eq. (5.2.2), so the integrals between 
complex conjugate pairs of branchpoints are still expressed by eqs. (5.1.13). 
Note that here, though, it is not possible to find the other integrals between other pairs of 
branchpoints with a change of coordinates (5.1.27) as in Remark 5.1.3: the images of the 
branchpoints after such a change of coordinates would not be complex conjugate pairs any 
longer. Later in this section we propose a solution to this. 

As for case 2, eq. (5.1.12) does not hold for the first step of the recurrence, due to the 
different ordering of eq. (5.2.3); it needs to be modified as follows, for n = 1, in view of eq. 
(5.2.3): 

fli ~ u, a'l — V, bi = w b'l — u' , Ci — v' , c'l — w' , 
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Note that this change only occurs in the first step: after that, the curve C has all real branch- 
points, and hence the Richelot-Humbert iteration can be applied as in [BM88] . In particular, 
let us denote by I{p, q) the integral on C between p, q on the first sheet, and attach a super- 
script to denote integrals on the curve C^^^ of equation y'^ + Pi{x)Qi{x)Ri{x) = 0. Then 
the integrals /(*^ can be expressed by equations (5.1.13), using the AGM method for the curve 
C. 

We obtain an expression for the integrals /(p, q) on C as follows. Equation (5.1.6) suggests 
that eq. (5.1.15) for n = is modified as follows (for the three pairs): 

I{a,a') = IqI (a^,5i), 

I{b,h')^m\c\,ai)^ to(-/'(a'i,6i)+/'(6'i,ci)), (5.2.4) 
I{cc!)=-tol'{b'^,ci). 

These relations are suggested by analogy with the Bost and Mestre case in view of the differ- 
ent ordering of eq. (5.2.3), and have been checked numerically in a variety of examples. In 
particular, in the Maple file iniages_paths_correspondence .mw we have considered^" the im- 
ages on C" of the contours of integration on C, i.e. straight lines between the branchpoints, 
via the correspondence (5.1.7), to obtain the contours of integration on the right hand side 
of eqs. (5.2.4). 

Moreover, from a direct analysis of the images of the paths between every branchpoints, 
again performed numerically, relations analogous to those above for every pair of branch- 
point emerge: 

/(a,fe)=ito(/'(ai,a'i)+/'(&'i,Ci)), 

I(a',h')=\t^ {l'{a,,a[)-l'{b[,ci)), 

\ (5.2.5) 

c) = -to (-/'(ai, a'l) + I' {al^ " ^'i)), 
nH, c') = ^to (-/'(ai, a'l) - foi) - fo'i)). 

Recalling that we are able to express the integrals /'(p, q) applying the AGM method to the 
curve C" (with real branchpoints) as in Theorem 5.3 (combined with remark 5.1.3), eqs. 
(5.2.4) and (5.2.5) can be considered an extension to the AGM method, which gives all the 
integrals between branchpoints on C. 



'"More specifically, we have studied the images on the curve C" of equation Aj/^ + U {x)V{x)W(x) = 0, in order 
to understand the first equality of (5.1.15) (as the second foUows immediately from the first). We hence obtain, for 
instance, that the straight line between a and a' on C, call it 71 (a, a'), is sent to a closed cycle encircling a'j and 61 
on C': noticing that in the second equality of (5.1.15) there is a factor of 1/2, absorbed in the definition of ti, we 
obtain that the image of 7(0, a') on the curve + Pi {x)Qi {x)Ri (x) = 0, i.e. C , is precisely the path from a'^ and 
bi. 
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Finally, we can apply considerations similar to the above also in case 1 of eq. (5.2.2), to 
find the integrals between non complex conjugate pairs of branchpoints. As before, we can 
express integrals between branchpoints on C in terms of those on C": 

/(a, a') = t^I {ai,a\), 
I{b,b')=tol'{h,b[), 
I{c,c')^tol'{ci,c[), 

I{a, b) =i to ( I'{a[,bi) ~ /'(ci, c'l)), (5 2.6) 

Iia,b') =^to il'iai,bi) + l'ici,c[)), 

I{b, c) =i to ( - /'(fo'i, ci) - /'(ci, c'l)), 

/(6',c')=^to (/'(5i,6'i) + /'(6'i,ci) + /'(ci,c'i))). 

5.2.2 General remarks 

We point out that the results obtained above hold true not simply for the monopole quotient 
curve, but for every curve satisfying either (5.2.2) or (5.2.3). We summarise the findings of 
this section as follows: 

Extended AGM method 1. The AGM method can be extended to the case of a curve C with 
complex conjugate branchpoints, as follows: 

1. perform the first Richelot transformation to obtain a curve with real branchpoints, C; 

2. use the AGM method ofeqs. 5.1.13 to find the integrals between branchpoints on C; 

3. use eqs. (5.1.11) to calculate u, u' , v, v' , w, w' , and examine their relative ordering: there 
are several possibilities 

(a) easel v < w < w' < u < u' < v' (same ordering as in Bost-Mestre case [BM88]): all 
the integrals between branchpoints on C are given by eq. (5.2.6); 

(b) case2u < v < w < u' < v' < w' : the integrals between branchpoints on C are given 
by eqs. (5.2.4) and (5.2.5). 

The other cases can be obtained in a similar fashion. 

We remark that here we have only examined two possible orderings of u, u' , v, v' , w, w'. 
To obtain the integrals between branchpoints on C also for other orderings, one can per- 
form a similar analysis of the images of the contours of integration, i.e. straight lines between 
branchpoints, via the correspondence (5.1.7), to obtain the contours of integration on C. 
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5.3 The AGM method for the quotient monopole curve: solv- 
ing the Ercolani-Sinha constraints 

We apply the results of the previous section to the quotient monopole curve X of eq. (3.1.8), 
where we scale /3 to be 1, to get 

^ {X^ + aX + gf +4. (5.3.1) 

If we split the branchpoints In complex conjugate pairs as In (5.2.1), to find the integrals 
between branchpoints on X we work in steps, as explained in the previous section 5.2.2: we 
perform the first Richelot transformation to go back to a curve with real branchpoints, X'; 
we can use eqs. 5.1.13 to find the integrals between branchpoints on X'; finally, to get all the 
integrals between branchpoints on X, substitute these in eqs. (5.2.4), (5.2.5) for a > 0, or in 
eqs. (5.2.6) for a < 0. 

Hence, for both cases we are able to use the AGM method to find all the integrals between 
branchpoints, and therefore, upon using the arc expansion (3.8.3), to find the periods on X. 
This has been implemented in Maple, for this specific curve: we developed an algorithm 
that, for fixed values of a and 7, performs the AGM method as described above for case 1 and 
2, giving as final output the periods on the basis of section 3.8. These have been checked 
against the exact results obtained using Northover's code (see Appendix C) for a range of 
values of the parameters, and agree within 10"''. The advantage of the AGM method is that 
is much faster, as it deals only with polynomial manipulations; moreover, the convergence 
of the sequences (a„), {a'„), (6„), (c„) only needs very few steps, usually 6 or 7, for the 
precision we require. These considerations allow us to use this method successfully in the 
iteration, described in the next section, to solve the Ercolani-Sinha constraints. 

Change of coordinates 

For convenience, we rescale /3 to be 1. This is obtained using the following change of coordi- 
nates for the curve X 

w^W ^ a^a = P^^^^a, 1^9 = r^l, (5.3.2) 

and for the curve X 

X X = p-^'^x, y->Y ^ p-^y. (5.3.3) 

For our purposes the only effect of this is on the constant v in the Ercolani-Sinha constraints. 
Indeed, recall that Theorem 1.7 requires that -25ik = iVi a+m S differentials 
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n^'-' = — gp — , fl^''' = — gp — , Using the parametrisation (5.3.2) we then have 

dw 



^ 3w2 + az^ L + aZ^ 



7 



Imposing -2 = /^.g+^.j, ^l^'^K we have that 

This means that the genus 2 Ercolani-Sinha constraints (3.11.2) read 



(5.3.4) 



with c as in eq. (3.11.4). The first of these does indeed impose a condition on a and g, which 
we will solve in the next section; the second will be used to find /3, and hence rescale back to 
a and 7. 

5.3.1 Solving the Ercolani-Sinha constraints 

We are now in a position to use the AGM method to solve the Ercolani-Sinha constraints, 
namely to find those values of the parameters a and 7 in the equation 5.3. 1 for X that satisfy 
the constraints in the form (5.3.4). We focus on the first of these, namely 



^dX = ( no / +3n / +3mo / +3m / ) Ui = 0, (5.3.5) 

ao J J ho <J bi, 



Y 



and proceed as follows. 

We know already from the early work in [HMM95] that there is a particularly symmetric 
monopole whose spectral curve belongs to the class 5.3.1: this is the tetrahedral monopole 
(1.4.5), corresponding to a = 0, g = ±5\/2 (in rescaled quantities). We then seek a curve of 
solutions containing this monopole, i.e. passing through the point a = 0,g — 5a/2, in the 
space of parameters a, g. The work of Braden and Enolskii in [BE06] (see also section 1.6.2) 
gives a solution to the Ercolani-Sinha constraints for the case of the tetrahedral monopole, 
corresponding to the following values for the integers in (5.3.5), expressed in the cyclic basis 
of Figure 3.14 (see section 3.11) 

m = 1. n = I, no = 5n ~ mg = —3 n. (5.3.6) 

As we seek a curve of solutions passing through the point a = 0,g = 5^/2, by continuity, we 
expect that on this curve the constraints (5.3.5) are satisfied with the same set of integers 
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Figure 5.2: Solutions to the Ercolani-Sinha constraints I 



(5.3.6). 

Note that for the tetrahedral monopole with a = and g = — 5\/2 the Ercolani-Sinha 
constraints are satisfied for different values of integers, namely 

m = 0, n = 1, no = 5n — m, nio ~ —"in. (5.3.7) 

We hence expect another branch of the solution curve passing through the point a = 0, g = 
— 5\/2, with constraints satisfied by the integers (5.3.7). 

To begin, we wish to find solutions to (5.3.5) starting from the point a = and g ~ 5\/2, 
with integers (5.3.6) : we proceed iteratively as follows. We start varying a by a small e, namely 
flj = i • e ; we then vary g in smaller steps, gi^ = 5a/2 -|- fc • e^. For every such pair (oj, g^jc) 
we calculate the periods using the AGM method, and hence compute the first constraint in 
(5.3.5) , with integers given by eq. (5.3.6): for every ai we take the g^.^ for which this con- 
straint vanishes. We repeat this a sufficiently large number of times, obtaining the curve in 
Figure 5.2. We now comment on this graph, referring to the Maple file AGM_ES_monopole .mw 
for more detail. 

We have used a step of 10~^ for a, and of 10^^ for a e (2.8, 3.0) to obtain a greater detail in 
this interval. 
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For values of a e {0, 3}, the outcome is that we have a curve of solutions in the space of pa- 
rameters passing through the points (0, 5\/2) up to the point (3, 0), which, however, does not 
belong to this solution curve. 

We remark that the point (3, 0) is in fact a singular point, as 4 of the branchpoints collide 
pairwise, giving two singular points at ±i: this results in a rational curve, with equation 

y2 ^ (2,3 32;)2 + 4 = (2:^ + 4){x^ + if. (5.3.8) 

Note that the curves of equation (5.3.1) with 5 = and a ^ 0, namely = (x^ + ax)'^ + 4, are 
all hyperelliptic: the only exception is precisely the case a = 3 above. 

As all the cycles on the curve (5.3.8) are homologous to zero, the first of eqs. (5.3.4) is trivially 
satisfied, but the second is zero as well: hence the Ercolani-Sinha constraints are not satis- 
fied for this curve. Thus the values g = and a = 3 do not correspond to a monopole. 

We are also able to extend this curve for negative values of a, on the left of the point 
(0, 5\/2): in Figure 5.2 we plot 100 points corresponding to a < 0, g < 0. 

When trying to extend the solution curve to g < 0, we do not observe any values of the 
parameters satisfying the first of eqs. (5.3.5) with integers (5.3.6). But we notice that, since 
the point (3, 0) does not belong to the solution curve, continuity arguments do not prevent us 
to use the different set of integers of eq. (5.3.6): with these, we manage to extend the solution 
curve through the point (0, — 5\/2), which corresponds again to a tetrahedral monopole. We 
point out that this arc of curve for g < is precisely the reflection with respect to the a-axis 
of the arc obtained for g > 0. 

Finally, using the second equation in (5.3.4), we can obtain f3 for each pair (a, g), again 
using the AGM algorithm; we prsent /3 as a function of a in Figure 5.3. This allows us to find 
the original quantities a and 7 in view of (??). 

5.3.2 Elliptic monopoles 

In section 4.1 we have examined the conditions under which the quotient monopole curve 
X admits a degree 2 elliptic subcover: the vanishing of the invariant xso gives a curve of real 
solution in the parameter space a, g, plotted in Figure 4.1. We point out here that the solu- 
tions to the Ercolani-Sinha constraints found above do intersect this curve. This is plotted in 
Figure 5.4, where we see clearly that there are four points that satisfy both the Ercolani-Sinha 
constraints and the xso vanishing conditions; their coordinates have been found numerically 
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-2-, 

-3- 



Figure 5.3: /3 as a function of a 

to be 

01 = 2.584590,31 = 0.795087, 

02 = 2.884209,32 = 0.209470, 

03 = 2.884478,33 = -0.209525, 

04 = 2.584985,34 = -0.796463. 

Tliis means that for the above values of the parameters the curve X covers 2 elliptic curves, 
and hence the explicit solutions corresponding to these two monopoles can be expressed in 
terms of elliptic functions. This will be object of further study in the near future. 

5.4 Hitchin constraint H3 

Until now we have dealt only with one Ercolani-Sinha constraint, namely that of Theorem 
1.6 (or equivalently Theorem 1.7), corresponding to Hitchin constraint H2, which ensures 
that the boundary conditions (1.1.6) are satisfied. Nothing has been said about the con- 
straint H3, which is equivalent to the nonsingularity of the monopole: we discuss this in the 
present section. 

Recall that Braden and Enolskii showed in [BE06] that H3 is equivalent to determine 
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Figure 5.4: Solutions to the Ercolani-Sinha constraints II 



whether the line sU - JC intersects the theta divisor 9 (see Theorem 1.9). This condition 
has indeed been checked in the case of the tetrahedral monopole in [BE09] : they use the 
same methods employed in this work of thesis to express the theta functions in terms of 
lower genus ones, which allows to check the above condition. Their conclusion is that the 
only monopole in the class (1.6.1) (corresponding to a = in our case) is the tetrahedral 
monopole (cf. Theorem 1.18). 

We use this result in conjunction with an analysis of the moduli space of cyclically sym- 
metric monopoles performed by Hitchin, Manton and Murray in [HMM95]. They show 
that cyclically symmetric monopoles form 4-dimensional totally geodesic submanifolds M^^ 
(where < Z < n) of the full moduli space of charge n monopoles. In our case, n = 3 and 
ignoring the rotational degrees of freedom we have a one dimensional submanifold. 
In the previous section we have established the existence of a one parameter family of curves 
satisfying the constraints HI and H2, containing the tetrahedral monopole curves. In par- 
ticular, the curve of Figure 5.2 were built in such a way that they describe the only possible 
curves satisfying the Ercolani-Sinha constrains passing through tetrahedral monopoles. 
Combining these two remarks, we expect that this conditions is satisfied on the whole curve 
of Figure 5.2. 

We point out that the constraint H3, namely the nonvanishing of 9{sIJ - K) for s e (0, 2) 
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translate into the vanishing of genus 2 theta functions. Once can show (see [BralO]) that 



it = 7r*(/Coo+ - e) 



Combining this with (3.11.3) we obtain that 



sU — /C = 7r*{sU — /Coo+ ^ s). 



Using Theorem 3.2 (see also eq. (3.13.1)) we have that the nonvanishing condition above 





can be expressed as the nonvanishing of the 3 genus 2 theta functions 



- 

3 " 



{z;t). This 



has indeed been checked numerically for several values of s, see [BDEIO], confirming the 
prediction above. 



5.5 Conclusions and discussion 

Before commenting in more detail on the results obtained so far, let us summarise the argu- 
ments presented in this work of thesis. 

The aim was to investigate the existence of charge 3 cyclically symmetric monopoles; this 
translates into verifying whether certain constraints on the coefficients of the spectral curve 
X are satisfied. The first step consists in using the C3 symmetry of the system to quotient 
the genus 4 curve X to get a genus 2 curve X. This quotient is an unb ranched cover, which 
allows us to use Fay-Accola theory to express several quantities on X in terms of analogous 
quantities on X: to do this, the use of a particularly symmetric homology basis is crucial. 
A fundamental upshot of this is that the Ercolani-Sinha constraints reduce to constraints 
among the periods of the genus 2 curve X. We focus in particular on the constraint corre- 
sponding to Hitchin condition H2: we show that this basically reduces to check the vanishing 
of § Solving this entails solving explicitly several hyperelliptic integrals: this is done by 

c 

adapting the so called AGM method. This method is well known in the case of real branch- 
points, we manage to extend it to the case of complex conjugate branchpoints. 
Using the AGM method iteratively, we are able to find a curve of solutions for H2 in the mod- 
uli space, and we argue that for these H3 is satisfied as well. Hence we have the main result 
of this work of thesis, namely we establish the existence of a one parameter family of cycli- 
cally symmetric charge 3 monopoles; we do so by finding explicit values for the parameters 
a, g, given in Figure 5.2 for which the curve (3.1.1) does describe a monopole. 

We conclude comparing our results with those of other authors. 
We remark that the behaviour of the solution curve in Figure 5.2 is consistent with the results 
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Figure 5.5: Log-log plot of the asymptotic behaviour of a versus 7 according to Hitchin, Man- 
ton and Murray (solid), Sutcliffe (dash) and here (dots). 

of [Sut97] , where Sutcliffe predicts that the curve (5.3.8) , describing a configurations of three 
unit-charge monopoles with dihedral symmetry, constitutes an asymptotic state for a 
3-monopole configuration (cf. eq. (4.16) in [Sut97]). Moreover, both Sutcliffe [Sut97] and 
Hitchin, Manton and Murray [HMM95] give the following picture of the scattering of three 
monopoles, corresponding to geodesic motion along one of these loci. Three unit charge 
monopoles come in at the vertices of an equilateral triangle, moving towards its centre, in the 
X1 — X2 plane. They then coalesce instantaneously into a tetrahedron. Finally the tetrahedron 
breaks up in a unit charge monopole moving along the positive x3-axis and an axisymmetric 
charge 2 monopole, moving along the negative xs-axis. The asymptotic behaviour at this end 
of the scattering is the following: 

a- (7rV4-362), ^ ^ 2b{b'^ + ir^ / A) , /3 - 0. (5.5.1) 

We plot this relation together with the set of points we obtain from solving the Ercolani-Sinha 
constraints, on a log-log plot; note that here we use the values of a and g of Figure 5.2, and the 
corresponding values of /3 to revert to a and 7. Figure 5.5 suggests that our results are indeed 
consistent with the asymptotic behaviour (5.5.1). Our results then confirm quantitatively the 
approximate predictions of [HMM95, Sut97]. 
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APPENDIX 

A 

Riemann Surfaces 



In this appendix we recall some basics about Riemann Surfaces, essentially to set notation 
and fix out conventions. For proofs of the theorems stated here, and a thorough exposition 
of these topics, we refer to [FK80, Fay73, Mir91], amongst the vast literature on the subject. 

A. 1 Cycles and period matrix 

Given a Riemann surface E of genus g, denote by ai, . . . , , bi, . . . , bg a canonical basis for 
Hi (S, Z), namely a basis such that its intersection matrbc is ± J, where 



Throughout this thesis we choose the minus sign, to be consistent with the conventions of 
[BE06] , as we aim to generalise the results therein. 

Denote by (ui, . . . , Ug) a basis of holomorphic differentials; we then define the matrices 
of periods A and B to be 



Given A and B we now construct the Riemann period matrix. We work with canonically 
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a-normalised differentials, the period matrix is 

r„ = BA-'; 

while for canonically b -normalised differentials it is rj, = AB~^; clearly T(, = r~^. We some- 
times denote the period matrix by r when neither normalisation is necessary. 



A.2 Theta Functions 

The canonical Riemann ^-function is given by 

d{z;T) = exp{i'Kri^Tn + 2i'jrz^n), (A.2.1) 

where r e N. The Riemann ^-function is holomorphic on C x and satisfies the following 
periodicity conditions 

e{z+p;T) = 6{z; r), 6{z + pr; r) = exp{-i7r(p'^Tp + 2z^p)} 6»(z; r), (A.2.2) 

where p e Z'". 

The Riemann ^-function 6a,b{z; r) with characteristics a, 6 e Q is defined by 

^'0,6(2:; r) = exp {in{a^Ta + 2aF{z + b)))] e{z + Ta + b; r) (A.2.3) 
= exp {i7r(n -|- a)^T(n -|- a) -|- 2z7r(n -|- a)^(z -|- 6)} , 



where a, 6 e Q'". An equivalent notation is 



(z;t). 



The function Oa,b{T) = Oa,b{0; r) is called the ^-constant with characteristic a, b. 



TheoremA.1 {Riemann theorem). If A is the Abel map based at a point Qo, andd isageneric 
point in J{T,), the function 9{A{P) + d) either vanishes identically, or has g zeros Pi,...,Pg, 
such that 

A{Pi) + . . . A{Pg) = d - K. (A.2.4) 
Here K is the vector of Riemann constants, independent of d but dependent on Qq and 
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on the homology basis. With this choice of sign convention for JC one finds 




A.2.1 Weiers trass reduction 

When a curve covers one of lower genus, the period matrix admits a reduction, i.e. can be ex- 
pressed in terms of a lower genus period matrix, and similarly, the associated theta functions 
can be expressed in terms of lower dimensional theta functions. This is known as Weierstrass 
reduction; a description purely in terms of the matrix of periods is given in [Mar92b] and 
more recently [BE01],[BE02]. 

Ag X 2g period matrbc n is said to admit reduction if there exist a maximal rank gi x g 
matrix A of complex numbers, a 51 x 51 matrbc of complex numbers IIi and a maximal rank 
2gi X 2g matrbc of integers M such that: 



where 1 < gi < g. 

When a period matrix admits reduction, it can be shown (Weiestrass' theorem) that there 
exists an element a e Sp{2g,Z) such that 



where r' and n have the properties of a Riemann period matrix and Q is a gi x (g — gi ) matrix 
with rational entries. 

Because Q here has rational entries, there exists a diagonal (g — gi) x [g — gi) matrix 
D = Diag((ii, . . . , dg^g^) with positive integer entries for which {QD)jk € Z. With (z, w) — 
{zi, . . . , Zg^,wi, . . . , Wg-g^ ) the thcta function associated with r may then be expressed in 
terms of lower dimensional theta functions as 



An = UiM. 



(A.2.5) 




(A.2.6) 



E 



m 



{Dw;Dt*D). 



9{z + Qm;Ti)e 



(A.2.7) 



m=(mi,...,mg_gj 

Q<mi<di-1 
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APPENDIX 

B 

A sketch of the Tretkoff and Tretkoff 

algorithm 



In this appendix we describe the construction by Tretkoff and Tretkoff, as in [TT84] , which 
is an algorithmic way to find a canonical basis for the first homology group of a Riemann 
surface. 

Its key elements are that the Riemann surface X is seen a branched cover of the sphere, and 
its Hurwitz system is used to find a cellular decomposition of X (in terms of CW complexes) . 
From this cellular decomposition one can find a system of generators for Hi {X, Z), and, once 
calculated their intersection matrix, obtain out of it a canonical basis. 
We now examine these elements in detail. 

B. 1 Details of the algorithm 

Let X be a Riemann surface of genus g, and let p : X — > be its branched cover of the 
Riemann sphere. The Hurwitz system consists of the following data: 

- n sheets; 

- t branchpoints &i, . . . , 5* in P^; 

- for every branchpoint, the monodromy group associated to it (cf. section 3.3). 
Also, we will use the following notation 

Z:=pi-{6i,...,6a, Z:=p-\Z). 
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Decomposition in CW-complexes 

One can choose on an arbitrary point a not coincident with any of the branchpoints 
61, . . . , fet. Joining a with all the bi one obtains t segments k. The 1- dimensional CW com- 
plex, or 1 skeleton, F of X is thus 

r = \Jp-\k). (B.i.i) 

We observe that the Riemann sphere minus the union of the k, i.e. D := — Uk is home- 
omorphic to the interior of a disc, namely a 2-cell: this is the prototype of the cells to be 
attached to the 1 -skeleton F to obtain a cellular decomposition of X. 
Recall that 

Hi{X,Z) = Hi{T,Z)/{dDi,...,dDn), (B.1.2) 

namely the quotient of Hi (F, Z) with the group generated by the 9Z)„ with the relation dDi + 
... + dDn - 0. 

Using property (B. 1.2) one can find an homology basis for X out of one for F: this is useful as 
it is indeed quite immediate to find a basis for Hi (F, Z), as we discuss in the following. 

A basis for i?i(F,Z) 

Let r e F be a maximal tree, i.e. a contractible subgroup containing all the vertices in F, and 
denote ei, . . . , the edges of F not in T. 

One can build a nontrivial cycle in the following way. For every two vertices ofT, Pi, Pj , there 
is a unique path in T connecting them, denote it so a cycle in F passing through Pi,Pj is 
the union of ^ and the edge e; connecting P^, Pj. In this way we obtain r cycles ci, . . . , c^, 
one for every edge e; in F not in T. 

It is immediate to see that the above is a basis for Hi{V, Z), since, contracting T to a single 
point, one gets a bouquet of R circles. 

Intersection numbers 

We briefly discuss the way the authors propose to calculate the intersection numbers for the 
cycles in Hi (F, Z), which is much simpler than counting "by hand", and has the advantage of 
being algorithmic. 

Consider an open set [/ e X such that: 

A\. T eU 

A2. it can be mapped homeomorphically and preserving orientation onto the unit disc, 
A3. dU meets the in 2 distinct points, Pi,Qi. 
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The idea is to replace the subarc of c, in the more complicated tree F with the chord con- 
necting Pj, Qi in this disc. After this replacement, the element in Hi {T, Z) represented by q 
is unchanged. Therefore the intersection between two cycles q, cj is the same as the inter- 
section of the 2 chords P^Qi and PjQj- The latter is indeed much easier to calculate, as it is 
shown in section B.2 for the case of the curve X of equation (3.1.1): the outcome of this is an 
r X r intersection matrix. 



A canonical basis for Hi {X, Z) 

In order to get a basis for Hi {X, Z) out of that for Hi (F, Z) found thus far, one would have to 
quotient by the group generated by dDj, and then consider the issue of canonicity. Tretkoff 
and Tretkoff remark that this can be achieved in a different, more algorithmic, easier way. 
Viewing the cycles q as elements of Hi {X, Z) , their intersection matrbc is the same as the one 
calculated above. But an intersection matrbc K for a genus g surface between an arbitrary 
number of cycles can always been put in standard form 



A.K.A' = 



V 



(B.1.3) 



Hence the matrix A selects 2g cycles among the set of generators above, such that their in- 
tersection matrix is that required; in detail: 



Cj, 



(B.1.4) 



We remark that the fact that the intersection matrix among these cycles is non-degenerate 
also ensures that they are independent. 



The algorithm 

The previous method can also be implemented algorithmically: we describe this, and apply 
in detail this algorithm to our case. 

The idea is to give a way to construct, algorithmically, a "planar version" A of the 1 -skeleton 
F, i.e. a tree from which one will obtain F after a quotient; then, since A is planar, it is easier 
to find the intersection matrix. 
The algorithm works as follows. 

1. Consider as basepoint a point on the Riemann sphere, call it a, and call its preimages 
under p on sheet fc; choose one of them, say ai on sheet 1. Then draw a family of circles 
d with center in ai , oriented clockwise. To construct the planar version of F, start from 
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ax, select points on d and join them to previously constructed vertex on Ci_ i . Starting 
from Ci, draw t points on this circle, and label them 61, . . . , 64: the labeling suggests 
that these will be identified with the branchpoints. Then join the point ai to them. 

2. On C2, for every bi on Ci put d vertices, as many as the elements of the monodromy 
group at that branchpoint; label them with the number of the sheet, in the same order 
as they are reached applying a permutation, i.e.ki tti (1), A:2 + 7r^(l ),..., fc^ = 7rf (1) 
(see example in section B.2). This means that for every ramification point on sheet 1 
one is applying a permutation to go to a different sheet, that is to a different a^. Then 
connect each hi with its corresponding vertices obtained in this way. 

3. On C3, for every point k on C2 coming from a 6j build vertices labeled again 

hi+i, . . . , 6f , 61, . . . , in this order, and connect the vertex k to these points. This 
means that one is following a path from each to a ramification point that is different 
from that one is coming from. 

4. Repeat this construction to get a sequence of base points on each C2k, a sequence of 
ramification points on C2k+i, connecting them as explained above, and stopping at a 
point with a label that has already appeared. 

The tree A built in this way is finite, and the 1 -skeleton V can be obtained as a quotient 
of A, identifying the edges between vertices with the same label. This is done in detail in 
the next section, where an explicit construction is given for a special case, and also the exis- 
tence of an open neighborhood U suitable for the computation of the intersection numbers 
is apparent, as explained above. 

B.2 The algorithm applied to the curve X 

We apply in detail the previous algorithm to our curve X of eq. (3.1.1): the result obtained 
following the steps above is given in Figure B.l. 
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B.2. The algorithm applied to the curve X 



The cycles obtained after quotienting are clearly labeled in the figure. Their intersection 
numbers are computed as explained previously: we consider an open neighborhood U' , in- 
dicated in the previous figure, such that its image under the quotient satisfies B.l; we map 
U' onto a disc, consider the points Pi , Qi on its boundary, image of the intersection of each 
cycle Ci with dU' and compute the intersections of the chords PiQi, PjQj, which is the same 
as the intersection of CiWith cj, as explained above. This is all represented in Figure B.l. 
From the figure, the intersection matrix is 
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(B.2.1) 



To have a canonical basis, this matrix needs to be put in the canonical skew- symmetric form. 
There is a well known algorithm that does so, and also gives the matrix of change of basis 
(see e.g. [Sie89], p.65): this has been implemented with Maple, and the result is that the 
intersection matrix above can be put in the form 



M.I.M^ = 



I 



^2,4 



^^4,2 
D4,2 
3)2,2 / 



(B.2.2) 
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where M is found to be 



M = 



B.2. The algorithm applied to the curve X 



1000000000 

1 00000000 

1000 -1 00 

0100-2 1 00 

000-10 10 

-1 1 -1 1 

-1 2 -1 -1 

1-1000 10 

1 -1 001 -1 00 



1 0-1100 



-1 1 



-1 1 



Thus, in term of the cycles defined earlier, the required basis is 



(B.2.3) 



(ci,C2, Cs - Cr, C4 - 2C7 + Cg, -C6 + Cg, -Ci + Cq - C7 + Cg, -C2 + 2 C7 - Cg - Cg, C2 - C3 + Cg). 
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APPENDIX 

c 

Maple tools: algcurves and extcurves 



Various results presented in this thesis have been obtained with, or rely on work done with 
Maple. In particular, the algcurves package for Maple provides a way to compute vari- 
ous objects on a Riemann surface; the ext curve code by Northover, together with its vi- 
sual aid cyclepainter, relies on this for implementing new tools which have been used 
throughout this work of thesis. Here we explain in short, referring to the relevant doc- 
umentation for more details, algcurves, developed in collaboration by Deconinck, von 
Hoeij and Patterson, has been part of Maple since Maple 6: an explanation of all the com- 
mands is included in Maple's help, and a more thorough description is given in [DvHOl]. 
Extcurves and cyclepainter, developed by Timothy Northover, can be downloaded at 
http://gitorious.org/riemanncycles, where the relevant documentation is also to be found. 



C.l algcurves 

The Maple package algcurves allows to calculate several objects on a Riemann surface. 
The Riemann surface is input as a polynomial in two complex variables, and is viewed as 
a branched cover of the Riemann sphere; singularities are admitted but the program works 
with the desingularised curve. From this, using standard results from Riemann surface 
theory, various quantities can be computed: the genus, a basis for the differentials, the j- 
invariant, etc. . Here we focus in particular on the basis for the first homology group, and the 
corresponding period matrix, as these are more relevant for this thesis. 
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Figure C.l: Paths chosen for the analytic continuation for X 



Monodromy 

Here a Riemann surface is viewed as a branched cover of the Riemann sphere, thus the 
first ingredient needed is the monodromy around each branchpoint. Recalling section 3.3, 
one starts by choosing a basepoint on the base P^, and ordering its preimages on S, 
P = (x^,j/f),...,P = {x^,y^). The basepoint is not allowed to be a branchpoint or a singular 
point on E, and the program chooses it so that it is far away from these "problem points" (see 
section 3 in [DvHOl] for a precise characterisation), taking a point with a sufficiently nega- 
tive real part. 

The next problem is the choice of the cut for each branchpoint: in this algorithmic approach, 
using simply straight lines is not feasible, as one needs a way to avoid it passing through 
"problem points". This is done by choosing a combination of segments and arcs as in Figure 
C.l (again, for a precise explanation see [DvHOl]). 

Now, for every branchpoint Bi one can build a path 7^ on that starts and end at x^, 
going along the cuts previously constructed and encircling the branchpoint anticlockwise. 
The n-tuple (yf , , . . . , y,^) is analytically continued along 7^, returning to the same set of 
w-values, but permuted {y^.(^iyy^^(^2)' ■ ■ ■ ^Va.in))- ^^^^ defines a permutation ai, the mon- 
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odromy data for the branchpoint B,; the permutations ctj generate the monodromy group 
based at . The analytic continuation is performed using the command Acontinuation: 
we describe it shortly as it was used in this work of thesis. 

Consider two points on P\ sufficiently close, x^ and x^ , and a path 7 connecting 
them; corresponding to x^ there is an ordered n-tuple {yf,y2, ■ ■ ■ lUn) ^- y{x'^)< i-^- the 
y-coordinates of the n ordered points Qi on S. Following the path 7 on from x^ to x^ , the 
entries of y(x^) also follow the lifts 7^ of 7 on S. If x^ to x^ are close enough, these lifts can 
be approximated as straight lines: 

y(x«) = y(a;Q) + y'[x^){x^ - x^) + 0{\x" - x'^p), 

where the entries of y'{x^) are the tangents to S at Q, = {x^ ,yf), obtained by implicit dif- 
ferentiation. If the two points are close enough, and the path 7 does not deviate much from 
a straight line\ the entries of y{x^) above give an accurate approximation for the n-tuple 
above x^, appropriately permuted: 2/^(2), ■ • • , y^^r.)) y{x'^')- 

Given arbitrary points x^ and x^ and a parametrised path "f{t) , the command Acont inuat ion 
iterates the above procedure along small portions of the given path, until the necessary con- 
ditions above are satisfied. It gives as output the ordered lists y{x^) and y(a;^), which means, 
for example, that the lift 7 of 7 with initial point the first entry of y(a:'3), corresponding to Qi, 
has final point the first entry of y(a;^), corresponding to i?cr(i)- 

The above results in the command monodromy (f , x , y) , whose output is a list with three 
entries: the first is the basepoint x^; the second is the ordered n-tuple (a^^ (1)' (2)' ■ ■ ■> 

(n))' the third is a list of the branchpoints Bi with the corresponding monodromy, given 
as a permutation ■ We give an example of how this work for the case of the curve X; the 
results are given to 4 digits for readability, but setting Digits in Maple yields the desired pre- 
cision. 

> f:=w"3 + w*z~2 + z~6 + 5*sqrt (2) *z~3 - 1; # define the curve 

> m:=monodromy(f ,z,w,showpaths) : 

#computes the monodromy (the option showpaths produces Figure B.l) 

> m[l]; # basepoint 

-2.0584 

> m [2] ; # sheets 

[-1.7952, .8976-2.5795*1, .8976+2.5795*1] 

> m[3]; # branchpoints and their monodromy 
[[[.9370-1.6904*1, [[2, 3]]], [.9954-1.6567*1, [[1, 2]]], 

'This can be made precise, e.g. bounding the second derivative of y(a;'3) to bound 0(|x^ — x'^\'^), 
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[-1. 9324-0. 337e-l*I, [[1, 2]]], [-. 2508- .4525*1 , [[1, 3]]], 
[-.2665-. 4435*1, [[1, 2]]], [. 5173-0 . 90e-2*I , [[1, 3]]], 
[. 5173+0. 90e-2*I, [[1, 2]]], [-. 2665+ . 4435*1 , [[1, 3]]], 
[-.2508+. 4525*1, [[1, 2]]], [-1 . 9324+0 . 337e-l*I , [[1, 3]]], 
[.9954+1.6567*1, [[1, 3]]], [.9370+1.6904*1, [[2, 3]]]] 



A homology basis 

The construction described above gives a way to view the Riemann surface as n sheets cut 
and glued along the cuts 7^ passing through the branchpoints , each with monodromy a b, ■ 
these are exactly the ingredients needed for the Tretkoff and Tretkoff algorithm of [TT84] (see 
Appendix B). Given the above data, the command homology (x,y) follows the Tretkoff and 
Tretkoff algorithm step by step; it gives as output both the cycles ci, . . . , c^, i.e. the homology 
basis for the l-skeletonT, and the canonical basis ai, . . . , Og, bi, . . . , bg for Hi{Y.,Z), together 
with the matrix of change of basis (in Hi (F, Z)) between the two. 

The cycle is given as a list: the first element specifies the starting sheet i; the second ele- 
ment is the coordinate of a branch point Bi in P^, together with the disjoint cycle of a, which 
contains i; the third element is a sheet j, and so on. Essentially, this reads: "from sheet i go 
to sheet j , by encircling the branchpoint Bi "; the list is cyclical, giving a closed cycle. This is 
h [cycles] below. 

The canonical cycles are also given as lists, but in a slightly different way: the first element 
specifies the starting sheet numbers; the second element is a branchpoint, together with 
a number indicating how many times one has to circle around it in the counter-clockwise 
direction (clockwise direction if this number is negative); the third entry is again a sheet 
number, etc. . This is h[canoiiicalcycles] below. 

> h:=homology(f ,z,w) : 

> h [basepoint] ; # basepoint 

-2.0584 

> h [sheets]; # sheets 

[-1.7952, .8976-2.5795*1, .8976+2.5795*1] 

> eval (h [cycles] ) ; # cycles 

table([l = [1, [.9954-1.6567*1, [1, 2]], 2, [-1 . 9324-0 . 337e-l*l , [1, 2]]], 

2 = [1, [.9954-1.6567*1, [1, 2]], 2, [-. 2665- . 4435*1 , [1, 2]]], 

3 = [1, [-.2508-. 4525*1, [1, 3]], 3, [. 5173-0 . 90e-2*l , [1, 3]]], 

5 = [1, [-.2508-. 4525*1, [1, 3]], 3, [-. 2665+ . 4435*1 , [1, 3]]], 

4 = [1, [.9954-1.6567*1, [1, 2]], 2, [. 5173+0 . 90e-2*l , [1, 2]]], 

7 = [1, [-.2508-. 4525*1, [1, 3]], 3, [-1 . 9324+0 . 337e-l*l , [1, 3]]], 

6 = [1, [.9954-1.6567*1, [1, 2]], 2, [-. 2508+ . 4525*1 , [1, 2]]], 
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10= [1, [.9954-1.6567*1, [1, 2]], 2, [.9370-1.6904*1, [2, 3]], 

3, [-.2508-. 4525*1, [1, 3]]], 
8= [1, [-.2508-. 4525*1, [1, 3]], 3, [.9954+1.6567*1, [1, 3]]], 
9= [1, [.9954-1.6567*1, [1, 2]], 2, [.9370+1.6904*1, [2, 3]], 
3, [-.2508-. 4525*1, [1, 3]]]]) 



> eval (h [linecircombination] ) ; 
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> eval(h[caiionicalcycles] ) ; # canonical cycles 

table ([b [3] = [[1, [.995401-1.656656*1, 1] , 2, [-1. 932407-0. 33715e-l*I, -1], 
1, [-.266483-. 443510*1, 1], 2, [.995401-1.656656*1, -1], 
1, [-.250849-. 452536*1, 1], 3, [-.266483+. 443510*1, -1], 
1, [-.250849-. 452536*1, 1], 3, [.995401+1.656656*1, -1], 

1, [.995401-1.656656*1, 1] , 2, [.937006+1.690371*1, 1], 

3, [-.250849-. 452536*1, -1]]], a[3] = [[1, [-.266483+. 443510*1, 1], 

3, [. 517332-0. 9026e-2*I, -1]]], b[4] = [[1, [-1. 932407-0. 33715e-l*I, 1], 

2, [.995401-1.656656*1, -1], 1, [. 517332-0. 9026e-2*I, 1], 

3, [-.266483+. 443510*1, -1], 1, [.995401+1.656656*1, 1], 

3, [-.250849-. 452536*1, -1]]], a [2] = [[1, [-.250849+. 452536*1, 1], 
2, [-.266483-. 443510*1, -1]]], b[l] = [[1, [.995401-1.656656*1, 1], 
2, [.937006+1.690371*1, 1], 3, [.995401+1.656656*1, -1]]], 
b[2] = [[1, [-.250849-. 452536*1, 1], 3, [.937006+1.690371*1, -1], 

2, [-.250849+. 452536*1, -1]]], a[4] = [[1, [-.266483+. 443510*1, 1], 

3, [-.250849-. 452536*1, -1], 1, [-. 250849+ .452536*1 , 1], 

2, [. 517332+0. 9026e-2*I, -1]]], a[l] = [[1, [.995401-1.656656*1, 1], 
2, [-1.932407-0.337156-1*1, -1]]]]) 

131 



C.l. ALGCURVES 



The Riemann period matrix 

A basis (ui, Ug) of holomorphic differentials on the Riemann surface is given by^ 

^ ^ Pk{x,y) 
dyF{x,y) 

with Pk{x, y) = J2i+j<d-3 (^kijx'y^, where d is the degree of F{x, y) as a polynomial in x and 
y. This is implemented in Algcurves by the command dif f erentials(f ,x,y). 



Now, given bases for the holomorphic differentials and for Hi (S, Z), the period matrix for 
the Riemann surface E is 

u={A,B); ^ Uj, = / uj. 

As the paths of integration are just straight lines and arcs of circles, they are parametrised 
by a; = j{t) with < < < 1. The lift of a; = 7(t), denoted hy y = 7(2/0, t), is obtained by 
specifying a starting value of y, yo, and by analytically continuing this value j/o along 7(i). 
These integrals are evaluated numerically using Maple's numerical integration routine. 



> differentials(f ,z,w) ; 



dz 



# basis of holomorphic differentials 

zdz 



> periodmatrixCf ,z,w) ; 

0.0278 - 0.0143 i 0.301 + 0.5269 i 
-0.4624 + 0.8095 i 0.3181 + 0.5522 i 

-0.4518 + 0.8229 i -0.6281 + 0.0068 i 

-0.4518 + 0.7961 i -0.6281 + 1.099 1 



# 2g X g period matrix 



0.0115-0.0360? 
-0.0115 - 1.065i 
0.0041 - 1.119i 
0.4592 - 0.2580 i 



-0.3174 + 0.5274 i 
0.3174 + 0.5162 i 
0.3066 - 0.5726 i 
0.0041 - 1.056 i 



> periodmatrixCf ,z,w,Riemcinn) ; 



-0.863- 0.048 i 
0.011 + 0.0110i 

0.0250 + 0.0055 i 
0.025 + 0.0167 i 

-0.435 - 0.863H 
0.435 + 0.2279 i 
0.017- 0.6156 i 



0.1206- 0.279 i 
0.1475 -0.3014i 
-0.3266 + 0.0159 i 
-0.3266- 0.6187 i 
-0.1549 - 0.2676 i 
0.1549 + 0.4722 i 
-0.4814- 0.0461 i 



-0.4770 + 0.360H 0.6531 + 0.1608 i 

# Riemann period matrix 



0.5477- 0.2433 i 0.0131 + 0.7222 i 
0.1712 + 1.2201 i -0.9952 - 0. 1650 i 



0.5704 + 0.9751, i 
-0.4986 - 0.2320 i 

0.0439 + 0.5845 i -0.9591 - 0.1642 i -0.3311 + 1.5383 i 

0.4113 -0.3140i 0.1469- 0.5196 i 0.2518 - 0.79099 i 



0.5027- 0.4584 i 
0.1609- 0.4954 i 
0.2554- 0.7874 i 
-0.8857+ 1.01434 i 



^We refer to [DvHOl] for the case of a singular Riemann surface. 
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C.2 Extpath and Cyclepainter 

In the present work on symmetric monopoles a fundamental simplification of the problem 
was the choice and use of a specific basis for the first homology group, namely one that ex- 
ploits the symmetry of the problem. This cannot be done with Maple's algcurves, as this 
package relies on the homology basis found via the Tretkoff and Tretkoff algorithm. There- 
fore the package extcurve has been developed, by Timothy Northover, to take care of the 
cases where symmetry plays an important role, which happens quite often in integrable sys- 
tems. 

The extcurves routines rely on algcurves, so, for instance, once again the Riemann sur- 
face is specified as the zeros of a polynomial of two complex variables. The single most im- 
portant feature of extcurves is the possibility of specifying an extpath, i.e. describing an 
arbitrary path on a Riemann surface, as a sequence of straight lines on the a;-plane together 
with a sheet specification at the initial point. 
In this framework, a homology basis is given as a list of extpaths. 

Paths can be created as extpaths in several ways: the one more used in this thesis is via 
CyclePainter. CyclePainter is a Java program to visualise paths on a Riemann surfaces, that 
can then be used by extcurve in Maple. The Riemann surface is pictured as a cut plane, 
where, very importantly, both the basepoint for monodromy and the reference point for 
sheet ordering can be chosen; the cuts are then straight lines from this basepoint to the 
branchpoints. The paths are drawn as a sequence of straight segments by simply pointing 
and clicking with the mouse; one can choose the starting sheet, and the segments have a 
different colour depending on their sheet; crossing a cut, hence changing sheet, results in a 
change of colour. 

This is very useful as it allows a graphical representation of a cycle, while also producing a 
file which can be directly read into Maple. 

We point out that this basepoints and cuts representation is just a pictorial representa- 
tion, when the path in exported as an extpath no use is made of that. 

Using CyclePainter, then, one can effectively draw a given homology basis for a Riemann 
surface (for example, the symmetric basis of section 3.7), and export it as a sequence of 
extpaths in Maple. Using extcurve, several functions can be made to act on extpaths. 

An important one is isect (curve, pathl, path2) , which computes the intersection 
number of pathl and path2. This has been widely used in this thesis to check if a given ho- 
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mology basis is canonical. 

Using this, another function is defined, findJiomology .transform (curve , homl , hom2) : 
this returns the matrix that transforms the homology basis homl to hom2. This has been rel- 
evant when comparing two different bases, like the cyclic and symmetric bases in section 
3.7.3. This is particularly useful when combined with the command 

trEinsf orm_extpath(initCurve , path, trans, finCurve := initCurve), which returns the 
image of path under the transformation"* trans. This last function allowed us to find 
the image of a given homology basis under certain involutions of the given surface, and 
findJiomology .transform resulted in the matrix form of these involutions on the given ho- 
mology basis (see for example section 3.9). 

Finally, the function periodmatrixC curve, hom) calculates the period matrix by trans- 
forming the algcurves version with the appropriate matrix of change of basis between 
algcurves basis and hom. This has been used to check numerically several results regarding 
the period matrix, e.g. its symmetry in section 3.9, or the Ercolani-Sinha constraints check in 
section 3.11. 

> marchCopen, "extcurves .mla") ; 

> with(extcurves) : 

> curveA, homA, namesA := read_pic("alphaO_A_f inal.pic") : 

# reads the cyclic basis from the corresponding CyclePainter file 

> curveBE, homBE, namesBE := read_pic("BE_hwb.pic") : 

# reads the symmetric basis from the corresponding CyclePainter file 

> Matrix (8, (i,j) -> isect (curveA, homA[i], homA[j])); 

# computes the intersection matrix for homA 
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caveat here is that this function only works with linear or Mobius transformations. 
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> M:=f ind_homology_trcinsf ormCcurveA, homA, homBE) ; 

# matrix of chsinge of basis from homA to homBE 



1 
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> Transpose (per iodmatrix(curveA, homA)); 

-0.0003 + 0.0369 i -0.0197 + 0.0358 i -0.0321 + 0.0192 i -0.1153 + 0.0734 i 

0.0328- 0.0194 i -0.0197 + 0.0358 i 0.00098 - 0.0367 i -0.1153 + 0.0734 i 

-0.0314- 0.0196 i -0.0197 + 0.0358 i 0.0325 + 0.0194 i -0.1153 + 0.0734 i 

-0.0000- 0.0000 i 0.0404- 0.0018 i -0.0000 - 0.0000 i -0.0057 - 0.4148 i 

0.0330- 0.019 i 0.0405 - 0.0007 i 0.0325 + 0.0187 i -0.0050 - 0.1388 i 

-0.0337- 0.0189 i 0.0405 - 0.0007 i -0.0329 + 0.0185 i -0.0050 - 0.1388 i 

-0.0006 + 0.038H 0.0405 - 0.0007 i 0.0002 - 0.0380 i -0.0050 - 0.1388 i 

-0.0000- 0.0000 i 0.0619 + 0.0322 i 0.0000 + 0.0000 i -0.3660 - 0.6182 i 

Note that here we see explicitly the symmetries described in section 3.9. 
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